INGENIEUR-ARCHIV 


UNTER MITWIRKUNG DER 
GESELLSCHAFT FUR ANGEWANDTE MATHEMATIK UND MECHANIK 
ZUSAMMEN MIT 


. _ A. BETZ : K.KLOTTER - E. METTLER - K. v. SANDEN 
aa 

_ — BF, SCHLEICHER - E.SCHMIDT - E.SORENSEN 

: ‘ : 

ae | 


HERAUSGEGEBEN VON 


R. GRAMMEL 
VIERTES HEFT 1957 


k A 4 oh sata AMOR IS si SF RE OIE SRE TET TS CTE ERY, mest 
SPRINGER-VERLAG - BERLIN / GOTTINGEN / HEIDELBERG 


Abgeschlossen am 6, Juni 1957 
Preis DM 11.40 


Das Fis , 
INGENIEUR-ARCHIV 


erscheint nach MaBgabe des eingehenden Materials zwanglos i in einzeln berechneten Heften, die zu Banden | 
vereinigt werden. 
Die fiir das Ingenieur-Archiv bestimmten Manuskripte sind unmittelbar an den Herausgeber 

Herrn Professor Dr.-Ing. Dr. R. Grammel, Stuttgart N, Robert-Bosch-StraBe 101 
oder an die Herren ) es 

Professor Dr.-Ing. Dr. A. Betz, Gottingen, Herzberger LandstraBe 39 A 

Professor Dr.-Ing. K. Klotter, Stanford (Calif.), 621 Alvarado Row. 

Professor Dr. E. Mettler, Karlsruhe-Durlach, Geigersberger Str. 12 

Professor Dr.-Ing, K. y. Sanden, Karlsruhe-West, Hertzstr. 16, (I. H. West). 

" Professor Dr.-Ing. Dr.-Ing. E. h. F. Schleicher, Dortmund, Plauener Str. 44 
Professor Dr.-Ing. E. Schmidt, Technische Hochschule, Miinchen, Arcisstr. 21 
Professor Dr.-Ing. E. Sérensen, Augsburg, MAN 


einzusenden. ; 


Die zum Druck angenommenen Arbeiten werden, soweit dies drucktéchnisch méglich ist, nach der: 
Reihenfolge ihres Eingangs beim Herausgeber veroffentlicht. 


Die Mitarbeiter erhalten von ihrer Arbeit zusammen 75 Sonderdrucke unentgeltlich. 


Fiir die Abfassung der Arbeiten wird auf das vom Deutschen Normenausschuf herausgegebene Heft 
Gestaltung technisch-wissenschaftlicher Veroffentlichungen™ hingewiesen. Die Vorlagen fiir Abbildungen | 
sind auf besonderen Blattern erwiinscht und konnen entweder in Reinzeichnungen oder in’ klarverstandlichen | 
Handskizzen bestehen; die Beschriftung und nétigenfalls die Reinzeichnung nimmt der Verlag vor. 


Mit der Annahme des pieinaltetst erwirbt der Verlag das ausschlieBliche Verlagsrecht fiir alle; 
Sprachen und Lander. 


Im ,,Ingenieur-Archiv“ erscheinende Arbeiten diirfen vorher an anderer Stelle nicht veroffentlicht | 
sein und auch spater nicht anderweitig verdffentlicht werden. 


Photographische Vervielfaltigungen, Mikrofilme, Mikrophote von ganzen Heften, einzelnen Beitragen 
oder Teilen daraus sind ohne ausdriickliche Genehmigung des Verlages nicht gestattct. 


Die Wiedergabe von Gebrauchsnamen, Handelsnamen, Warenbezeichnungen usw. ‘in dieser Zeitachrift | 
berechtigt auch ohne besondere Kennzeichnung nicht zu der Annahme, daB solche Namen im Sinne der Waren- ; 
ees und Blac honey unta peste bane als frei 2 zu betrachten waren und daher von jedermann benutzt werden } 

iirften. j 


i 
i 


Lay NEMLV EEL AE 


| 
; ; 
Heidelberg Berlin. W 35 
Neuenheimer LandstraBe 28—30 / Fernsprecher 279 01 Reichpietschufer 20 / Fernsprecher Sammel-Nr, 24 92 51 


| 
| 
Inhalt: 


Seite 


Kneschke, A., Rollreibung anf sBurbildender Fahrbahn. Mit 9 Abbildungen ... 227 


Stickler, B., Kin weiterer Beitrag zur Theorie der Bewegung rollender Fahr- 


zeuge. Mit 2 Abbildungen......... RO Mag us Abiko ks BE Ne i ae 244 
Scherer, A., KinfluBflachen einer Dreiecksplatte mit Aufpunkt am freien Rand. oo | 

Mit10 Abbildungen 96", ie ae ae te ae ae ee | 
Helmbold "Hi. B:, Gveineh des Tragfliigelauftriebs. Mit i abpideee Biot AN Pye 


Schwarze, G., Allgemeine Stabilitatstheorie der Schalen. Mit 3 Abbildungen ... 278 


Schmieden, C., Nichtlineare Schwingungen bei zwei Freiheitsgraden. I. Mittei-. 


lung. Mit 5 Abbildungen........4...-¢..0000005 j : _ 292 


Slee ee eee eee ee te tele ee ee 8 


INGENTEUR-ARCHIV 


XXV. BAND VIERTES HEFT 1957 


Rollreibung auf spurbildender Fahrbahn 
. Von A. Kneschke* 


1. Einleitung. Die gleichférmig rollende Bewegung eines schweren Kreiszylinders (Rad) 
auf einer Ebene (Fahrbahn) kann auf zweifache Weise aufrecht erhalten werden: durch eine Zug- 
kraft, die dem Bewegungswiderstand (Rollwiderstand) das Gleichgewicht halt, oder durch ein 
dem Rad aufgepragtes Triebmoment, das mit dem zwischen Rad und Fahrbahn auftretenden 
Reibungsmoment (Rollmoment) im Gleichgewicht steht. Fiir diese beiden Falle soll die rollende 
Reibung eines starren Rades unter der Voraussetzung untersucht werden, daB auf der ebenen Fahr- 
bahn eine als sehr zahfliissig angenommene, spurbildende Schicht liegt, die sich verschiedenartig 
auf den Rollvorgang auswirkt, je nachdem, ob es sich um das gezogene oder getriebene Rad handelt. 
Die zu untersuchende Rollbewegung vollzieht sich demnach unter dem EinfluB dieser diinnen 
Zwischenschicht, so daB die Anwendung der Mechanik zaher Medien nahe liegt. Die im einzelnen 
angegebenen vereinfachenden Voraussetzungen erméglichen, wie sich zeigen wird, eine sehr weit- 
gehende Beschreibung des Rollvorganges. 

Fiir die Erscheinung der Rollreibung, die durch die unmittelbare Berihrung eines elastischen 
Rades und einer elastischen Schiene verursacht wird, hat unter Beriicksichtigung der auftretenden 
Normal- und Schubkrafte L. Féppl+ eine strenge Lésung angegeben. 

Im Falle der zwischen Rad und Fahrbahn auftretenden zihen Zwischenschicht, deren Dicke 
als klein gegeniiber dem Raddurchmesser vorausgesetzt werden soll, wird es darauf ankommen, 
die in dem Rollspalt zwischen Rad und Fahrbahn beim Uberrollen der Zwischenschicht auftretende 
Verteilung der Strémungsgeschwindigkeit und des Druckes zu ermitteln, aus denen sich dann der 
Rollwiderstand bzw. das Rollmoment ergeben wird. Es wird sich zeigen, da die hierfiir not- 
wendigen Betrachtungen in einer engen Beziehung zur hydrodynamischen Theorie des Walz- 
vorganges stehen. 

Die rollende Bewegung des Rades mit der Umfangs- 
geschwindigkeit U* iiber die Fahrbahn kann man als 
einen Walzvorgang durch eine um eine ortsfeste Achse 
drehende Walze (Rad) und eine mit der Geschwindigkeit 
U gleichsinnig bewegte Ebene (Fahrbahn) auffassen. Da- 
bei sollen die in der iiberrollten Schicht zwischen Rad und 
Fahrbahn auftretenden Schubkrafte an der oberen und 
unteren Begrenzung des Rollspaltes vom Radumfang und 
der starren Fahrbahn vollstandig aufgenommen werden. 
Mithin wollen wir ein vollstandiges Haften der zahflissigen 
Zwischenschicht an den Randern des Rollspaltes voraus- 
setzen und demnach fir die Strémungsgeschwindigkeit u 


die Randbedingungen 
Ee gai Ut ps U* 


annehmen. Dabei ist h die von der Koordinate x abhangige Spalthéhe (Abb. 1). Uber die Ge- 
schwindigkeit U* des Rades in ihrer Beziehung zur Geschwindigkeit U ist vorlaufig nichts aus- 
gesagt. Wir werden erst dann, wenn wir die allgemeinen Betrachtungen schlieBlich auf ve eingangs 
geschilderten Falle fir den Antrieb des Rades anwenden, Genaueres iiber die GriBe U im Ver- 
haltnis zu U erfahren. Die auf Grund des Verformungsvorganges auftretende Schichtenbreitung 
senkrecht zur Radebene, die die durch die Héhenabnahme verursachte Verschiebung der zahen 
Zwischenschicht teilweise aufnehmen wiirde, soll unberiicksichtigt bleiben. Es handelt sich dem- 


nach um ein ebenes Problem. 


Abb. 1. Rollspalt. 


* Vorgetragen auf dem 9. Internationalen Mechanik-Kongref in Briissel am 11. 9. 1956. 
1 L. Féppl, Die strenge Lisung fiir die rollende Reibung, Miimchen 1947. 
17 
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2. Lagebeziehungen. Unseren Betrachtungen liegt ein Rollspalt zugrunde, der durch die 
Belastung eines starren Rades vom Halbmesser r auf die als zahfliissig angenommene Schicht von 
der urspriinglichen Dicke hy entstanden ist. Die Radbreite sei b. Die im AufriB durch den Grenz- 
winkelbogen a, dargestellte obere Begrenzung des Rollspaltes ist durch die Hingrifftiefe ho —hy 
festgelegt, die im Spalt veranderliche Héhe h durch den Winkel «. Alle Abmessungen werden auf 
das in Abb. 1 eingezeichnete rechtwinklige (x, y)-System bezogen. Es ist zweckmabig, den Kreis- 
bogen, der den Rollspalt abgrenzt, niherungsweise durch einen Parabelbogen zu ersetzen, dessen 
Parameter gleich dem Raddurchmesser ist. Wir erhalten dann die in der x-Richtung veranderliche 
Spalthéhe 


x2 
h=h+ oe 
oder, wenn wir den Winkel « ~ x/r und die Abkiirzungen 
gaol ae eis ee 
reas 1 yon ee r 
einfiihren, 
R=, = +o. (a) 


Fiir die im Rollspalt verainderliche GroBe (h — h,)/h = 0, die fiir h = hy in die relative Eingriff- 
tiefe 09 = (hy — hy)/ho, die wir kurz als Ein griffziffer bezeichnen wollen, iibergeht, gelten die 
Beziehungen . 


a h = A hy — 
Sayi—Hafs, 4=/P=yl—s. (b) 
Wir werden spater durch die Beziehung 
Ate Pp =a, (c) 


eine Hilfsvariable ~ einfiihren. Fiir diese gilt nach (a) 


(d) 


h 
COS Py = 3 =y1—6, 
Aus (a) und (c) folgt 
= hy 
ee eo (e) 


3. Formulierung der Aufgabe und ihre Integration. i eo i afte i 
Rollspalt zu praca ne wir von der Dieeacnulacnee Seen ee Seam 
eu _ i1 dp 

By de (1) 
aus, die von A. Sommerfeld! als Grundlage fiir die hydrod i i i 
verwendet wurde. Sie stellt eine Shee Fi cheneee cotropen Dikeacei a Seen 
der x-Komponente u(y) der Geschwindigkeit bei einer Strémung in einem von waralleleet Ebe : 
begrenzten engen Spalt dar. Dabei ist 7 die dynamische Zahigkeit der Flissigkeit. Man canna 


Differentialgleichung (1), wie in einer Arbeit von G. Duffing? bestatigt wird, aber auch dann an- | 


wenden, wenn die Strémung in einem Spalt erfolgt, dessen Begrenzungsebenen schwach geneigt 


sind. Dann ist p(x) nicht mehr linear von x abhangi i indigkei | 
é gig, und die Geschwindigkeitsko t 
ist auBer von y auch von x abhangig. Die Geschwindigkeitskomponente in ie Pa es | 


dabei vernachlassigt, und 0u/dx? wird gegeniiber 0?u/dy? als klein betrachtet. Diese Voraussetzungen 


sind bei dem nach oben parabolisch begrenzten Rollspalt weitgehend erfiillt, so daB wir berechtigt | 


sind, die Differentialgleichung (1) auch auf die Strémung im Rollspalt anzuwenden. 


* A. Sommerfeld, Z. Math. Phys. 50 (1904), S. 97 
* G. Duffing, Z. angew. Math. Mech. 4 (1924), S. 296. 
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Der Rollspalt selbst soll sich in allen weiteren Betrachtungen entsprechend Abb. 1 tiber das 
Intervall h, = h<h, erstrecken. Auch darin liegt eine die Rechnung vereinfachende Voraus- 
setzung. 

Unter Beriicksichtigung der Randbedingungen u =—U und =— U* 
die Liésung von (1) par E mi ae E Bi ee 
U— U* 

FP aa abi ar 06) Ne it amare i Oe (2) 


worin h nunmehr als Funktion von x anzusehen ist. Die DurchfluBmenge Q in der Zeiteinheit 
auf der Radbreite 6 ermittelt man aus 


h 
Q =—b fudy 
i 


Zu 
poe UF h® dp 
0 o( 2 ht wea ae) 


Die Kontinuitat des Fliissigkeitsstromes im Rollspalt ist gewahrleistet, wenn die durch seinen 
Querschnitt strémende Menge an jeder Stelle x einen konstanten Wert hat. Daher mu$ die Druck- 
verteilung p(x) so gewahlt werden, daB Q von x unabhangig wird, daB also 


U+U*dh , at #\=0 


2, dx |! dx 12 » dx 


ist. Dies liefert, wenn h* als Integrationskonstante eingefiihrt wird, das Druckgefalle 


d kes 
© —-—6n(U+ ut. (3) 


Die Konstante h* stellt damit die Spalthéhe dar, fiir die der Druck extrem wird. Der konstante 
DurchfluB ergibt sich demnach zu 


U = U* 
O=b for: (4) 
Aus (2) folgt, daB im Querschnitt bei h* die lineare Verteilung der Strémungsgeschwindigkeit 
a UN 
Ug 


herrscht. Wir wollen diesen Querschnitt FlieBscheide nennen. Bezeichnet man namlich die 


mittlere Strémungsgeschwindigkeit iiber einen Querschnitt mit U, so gilt die Bezichung 


Q — 
ae 
die in Verbindung mit (4) zu 
oe UL -U h* 
pe aie ie 


fibrt. Man erkennt daraus, daB links der FlieBscheide die mittlere Geschwindigkeit U gréBer ist 
als rechts davon. An der FlieBscheide selbst ist 
_ — ** 
aiea 
Es bestehen demnach im Rollspalt zwei hinsichtlich der parabolischen Geschwindigkeitsverteilung 
verschiedenartige Gebiete: ein Gebiet der Voreilung links von h* und ein Gebiet der Nacheilung 
rechts von h*. Im Bereich der Voreilung ist das Druckgefille positiv, im Gebiet der Nacheilung 
negativ, an der FlieBscheide erreicht der Druck einen maximalen Wert. 

Um die FlieBscheide festzulegen, benutzen wir die Tatsache, dafs am Radeingriff h=h) und 
am Radauslauf h=h, jeweils der Druck p=0 herrscht. Wir setzen dabei voraus, daB die nicht 
durch den Rollspalt abstrémende Restmasse der Schicht von der Dicke hy am Radeingriff nach 
beiden Seiten der eingreifenden Radbreite abstrémt, ohne dafs cin nennenswerter Staudruck dabei 
entsteht. Die DurchfluBmenge, fiir die Q< bUh, gelten muh, wird links der Auslaufstelle bei 
x — 0 auf der Fahrbahn abgesetzt und wird im allgemeinen eine Spurtiefe hy) —h verursachen, 
die kleiner als die Eingrifftiefe h, — h, ist. Fiir die GréBe des Strémungsdruckes wird vornehmlich 
der von h, bis hy sich erstreckende Rollspalt maSgebend sein. 


Mae 
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Aus (3) erhalten wir mit x = ro unter Beriicksichtigung von (c) und (e) 


dp 24 (U + U*) 


2 ie oa 
Gen ee (ee cos ®) 


woraus durch Integration iiber den Spalt die Beziehung 


Po 


ae (Po + sin Y cos Yo) + sin Yo cos? os 


67 (U > U* 

J ap Po) Pe ate z= too Sin Pp COS Py) — 

6 

folgt. Wegen des Verschwindens von p an den Randstellen ergibt sich daraus 
5 Po + SIN Yo COS Po 
La 


Fiihren wir hier die Eingriffziffer 6) ein, so wird nach (d) 


3 5 ; 
zy Po F Sin Po COs Fo) + Sin Po COS* Yo 


hes are sin /5) + oy 1 — dp) 


Te 


159 (1 — 0) 


- (are sin Jd) + Voy A 5) 


Um die Lage x* der FlieBscheide im Rollspalt von der Lange 


59 


0,000 
0,025 
0,05 
0,10 
0,15 
0,20 
0,25 
0,30 
0,35 
0,40 
0,45 
0,50 
0,55 
0,60 
0,65 
0,70 
0,75 
0,80 
0,85 
0,90 
0,95 
1,00 


0 G2 OY 06 08 10 
* * * 
Abb, 2. ee ey Me und ae als Funktionen von 6). 
Ip 155 1 


der Randbedingungen p(0) = p(y) = 0. Man erhalt 
r A$ [ 
|— 


3 h* ; 
Pox U*) see (p + sin p cos @) + sin p cost pore? + sin ~ cos ~) 


Mh 


L=Too = 


h* =h, += 


Tabelle 1 


h* 


x* 


hy 


1,000 
1,008 
1,017 
1,034 
1,051 
1,069 
1,087 
1,106 
1,124. 
1,142 
1,161 
1,180 
1,199 
1,218 
1,236 
1,254. 
1,271 
1,288 
1,303 
1,316 
1,327 
1,333 


hy 


1,000 
0,983 
0,966 
0,931 
0,894 
0,855 
0,815 
0,774 
0,731 
0,686 
0,639 
0,590 
0,540 
0,487 
0,433 
0,376 
0,318 
0,258 
0,195 
0,132 
0,066 
0,000 


l 


0,577 
0,572 
0,566 
0,553 
0,540 
0,526 
0,512 
0,496 
0,480 
0,463 
0,444 
0,425 
0,404 
0,381 
0,357 
0,330 
0,301 
0,268 
0,231 
0,187 
0,128 
0,000 


V2 7 hy dy 


zu ermitteln, gehen wir von der Beziehung | 


x2 


28 


A 


1,000 
1,000 
0,999 
0,997 
0,993 
0,987 
0,980 
0,969 
0,956 
0,940 
0,921 
0,898 
0,870 
0,836 
0,797 
0,751 
0,696 
0,630 
0,549 
0,449 
0,313 
0,000 


(5) 


(6) 


[aes 


A are sin V% 


CO 
6,220 
4,325 
2,957 
2,334 
1,954 
1,688 
1,489 
1,332 
1,203 
1,095 
1,003 
0,923 
0,853 
0,791 
0,736 
0,686 
0,642 
0,601 
0,564 
0,531 
0,500 


g, aus und erhalten mit h, = hy (1 — 6,) 


se ESE ie 1 


Abb. 2 zeigt den Verlauf von h*/h,, h*/hy 
= (h*/h,) (1—6,) und x*/I als Funktionen von 6); die Werte sind der Tabelle 1 zu entnehmen. 
Zur Berechnung der Druckverteilung im Rollspalt integrieren wir (5) unter Beriicksichtigung — 
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und daraus nach (d) 


i 3 oe ey oe ae 
Pome + UF — [a aresinyé + OB) + JL—9P 5 — 2 (aresin 5 + 80-9). (7) 
Man sieht, daB fiir 6 = 0 und 6 = 6, der Druck tatsichlich verschwindet, wenn h*/h, nach (6) 
eingesetzt wird. 


4. Resultierende Druckkraft und Reibungsmoment. Die Wirkungslinie der Resultierenden P 


der radial am Radumfang wirkenden Druckkrafte bprda und der tangential wirksamen hydro- 
dynamischen Reibungskrafte 


—bqrdx =brrda =by = r dx 
oy Metel 
bildet mit der Vertikalen den Winkel « (Abb.3). Dabei betragt 
die am Radumfang auftretende Schubspannung 7 nach (2) 


__ U—U* | h dp 
Ane De eee a a: (8) 


Die resultierende Druckkraft hat dann die Komponenten 


& Xo 
Psino = br f{ psinadx + br {[tcosada, 
6 a 


&q X 
Pcosa = br f{ pcosadx—br ftsina da. 
6 6 


Abb. 3. Krafte am rollenden Rad. 


Unter Beriicksichtigung der Randbedingungen fiir p wird 


a 


%o 0 Xo & 
; d CO. 
J psina dx = P cos x dx , | peosadx =—{ ®sinada, 
6 sae? 6 go 


so daB 


Xo 


Psina= br | (FF + tJeosada, 


0 


f 
st. Mittels (8) ergibt sich hieraus 


Psing =br [ B(1 +z, eosada+ bry (U—U%) | FS da. 
0 


%o 


Peost=—br | 71 + gy) sina da — Brg (U—U*) | Eds. 
¢ 0 


Bei der Anwendung dieser Betrachtungen auf das rollende Rad sind wir berechtigt, 


2 ae 
SILK Sp COS On, Ls 
zu setzen und h,/2r gegen Eins zu vernachlassigen. 
Dann erhalten wir naherungsweise 


h=h +50? 


Xo a2 
i 


2 2 
pense { #fi—S)a 4 ornew—oy | —— da, 


Xo ae 
P cosa =—br { Padx—bry(U—U*) | = do 
6 6 
und fiir ‘ 
h = 9 a=A, tgp 


~~ cos? p 
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mit Beriicksichtigung der Randbedingungen fiir p | 


Po 
AC bri? dp * 1 
Psing =— "7 | ate p dp +72 (U— U*) > | 


Po 
Peosk =—briy | Prgp dp + 2b (U—U") In cos go. 
6 


Wenn man (5) einfiihrt, folgt daraus durch Berechnung der Integrale 


* h* q , 
Pang te Po — Sin Po COS Yo — Fz (Po — SID Yo CoS Po + 2 sin? qo COS Pp) 
L it 
2b (U — U* 
a a eee 
P cos & = PCE EY aint gy] —gp_ (1 + cos gy)| + 2.6.0 (U— U*) In cos gy . 
al 


Fiir das Reibungsmoment beziiglich der Radachse ergibt sich in ahnlicher Weise 


oe dp 2br7n(U — U*) 
Mabe [gdx=— ie cozy Ay Po » 
6b U + U* : 2brn(U— U*) 
yim dit a arg eer: aS Po « 
1 1 1 


Wir wollen die eben ermittelten Groen in ihrer Abhangigkeit von der Eingriffziffer 6, mittels (d) 
anschreiben. Es ist dann 


Psinx = aoe are sin Vy — 5o (=) —Fp (are sin Jy + (2 d9—1) 759 (1— 6, (D6) 


2b (U— U*) 
i 


are sin Vo : 


-- 


P cos X= tay lesa a saa (9) 


Ni 6 ae + U*) are sin JO) — 37, (are sin 150 a 150 (I—6)) 


1 


2brn(U— U*) 
The, arc sin iets 


Diese allgemeinen Ergebnisse werden wir nun auf die eingangs erwahnten beiden Falle des gezo- 
genen und getriebenen Rades anwenden. Es wird sich dabei zeigen, daB die bisher noch will- 
kirliche Umfangsgeschwindigkeit U* des Rades in ihrer Abhangigkeit von der Fahrgeschwindig- 

keit U und der Eingriffziffer 6) festgelegt wird. | 


5. Das gezogene Rad. Da das gezogene Rad kein Moment aufnimmt, mu8 die Bedingung 


erfiillt sein. Das ergibt nach (9) 


h* are sin 5) +6 Gees) 
3(U + U*)| 1 — J 0 aE = ee 
See ( 2hy are sin |/d9 


woraus, wenn 


2 Sei 3) a — 9s) 


are sin Voo 


A= 


(10) 


gesetzt wird, 

U* A Osea ae ULAUMe, We 

U ~~ 2—A”’ U  - 2—A’ Op eae Co 
Damit sind U* und der relative Schlupf (U—U*)/U zwischen Rad und Fahrbahn in Abhangig- 
keit von der Eingriffziffer 5) festgelegt. Die Schichthohe h auf der tiberrollten Fahrbahn (Abb. 4) 
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betragt nach (4) 


so daB wir als relative Spurtiefe (hy—h)/hy = 69, die wir kurz Spurziffer nennen wollen, 
enh eae h* U* 
ee eae a os Ue 


(12) 


erhalten. Mit den in Tabelle 1 ersichtlichen Werten ist in Abb. 5 der Verlauf von A gezeichnet; 
mit den Werten von U*/U und 6, in Tabelle 2 ergeben sich die Bildkurven fiirU* /U und fiir die 


Spurziffer 5). Es tritt demnach ein mit der Kingriffziffer 
wachsender Schlupf zwischen Rad und Fahrbahn auf. Die 


Spurtiefe hy—h bleibt fiir 6, <0,7 merklich unterhalb der 


p N=Pcosa& 


A ? a, bp 


0 g2 O4 06 ~O7 O8 10 


Abb. 4. Krafte am gezogenen Rad. Abb. 5. Verlauf yon 4, U*/U und 6, fir das gezogene Rad. 


Eingrifftiefe hy) —h,. Fiir 6)> 0.7, also besonders groBen Schlupf, iiberschreitet die Spurtiefe ge- 
ringfiigig die Kingrifftiefe. Der Anteil des Gleitens am Bewegungsvorgang des gezogenen Rades 
ist in diesem Falle wesentlich gréfer als der des Rollens. Bei einer Abtrennung der Schicht durch 
eine iiber alle Grenzen wachsende Belastung kann das Rad nur gleiten, nicht rollen, da eine hydro- 
dynamische Reibung dann nicht vorhanden ist. Bei Punktberiihrung ist U=U7*, also reines 
Rollen, festzustellen. Beide Grenzfalle treten jedoch praktisch nicht ein. 


Die Normalkomponente N= P cos der resultierenden Druckkraft P = —_ tebele 
betragt nach (9) f Ee 5 
N P cosa 
bnU  byU 0,000 | 1,000 | 0,000 


5 ar & ss | 0,025 | 1.000 | 0,017 
=e i= Ain (Ol Sao) Ses : (ro ie 13a 0,05 0,999 | 0,035 
3A ira 0 Ge) — bo Dh, | 0) | ee) 0,10 | 0,994 | 0,072 
0,15 | 0.987 | 0,112 


Es zeigt sich, daB fiir 0,20 0,975 0,155 
Dip 0,25 0,960 | 0,201 

i 0,30 | 0,941 | 0,249 

0,35 0,916 | 0,300 

der Betrag von 0,40 0,888 | 0,353 
| —A 0,45 0,853 0,408 

25— In (I — by) 0,50 | 0,814 | 0,465 

0.55 | 0,769 | 0,523 

in dem Bereiche von 0<6,<1 unterhalb 1% des Wertes von one bea aa 


6 OAT een 
2—A\h,}/ 1—0) 0,75 0,533 | 0,756 
0,80 0,459 | 0,812 


liegt. Deshalb kann der erste Summand unberiicksichtigt bleiben, sodaB 0,85 | 0,379 | 0,865 


h* 0,70 | 0,601 | 0,699 
1—s-@— oy) 


: 0,90 | 0,289 | 0,915 
gan (1e2) osu tat, fe 0,95 | 0,186 | 0,961 
= a ee 2—6 13b 1,00 | 0,000 | 1,000 

bnU 2r- 2—A = Dh, | | Gir) 


ergibt. Durch die Beschrankung auf hinreichend kleine Werte von h,/2r wird auBerdem den 
eingangs erwihnten Voraussetzungen fiir die Anwendung der Differentialgleichung (1) Rechnung 
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getragen. Mit wachsender Schichtdicke hy oder abnehmendem Raddurchmesser 2r wichst bei 
konstanter Belastung N des Rades und konstanter Fahrgeschwindigkeit U die Eingriffziffer 6, 
und damit auch die Spurziffer Oy: Die Belastbarkeit bei konstanter Spurtiefe und F’ ahrgeschwindig- 
keit ist proportional dem Raddurchmesser und umgekehrt proportional der Schichtdicke. Mit 
eréBer werdender Fahrgeschwindigkeit nimmt bei konstanter Belastung die Spurtiefe ab. Eine 
Erhéhung von U erméglicht bei gleichbleibender Spurtiefe eine gréBere Belastung. Abb. 6 gibt 
den Verlauf der Belastungsfunk- 
tion (13b) wieder. 


Die Horizontalkomponente 
Z=P sina, 
die den Fahrwiderstand darstellt, 


ist auf Grund einer langeren, aber 
elementaren Rechnung 


Zo) hee 1 —Aarcsin|/ 6) 
bnU \; P 2A Ge 
(14) 
und wird durch Abb.6_ ver- 
anschaulicht. Die Zahlenwerte 


Neh, Z ho 
ba ae OE ee 
kénnen der Tabelle 3 entnommen 
werden. 
Tabelle 3 
; Pcosx hy Psing ,/ hy 
¥ bynU 2r bnU 2r 
0,00 0,000 0,000 
0,05 0,003 0,001 
0,10 0,012 0,004 
0,15 0,029 0,012 
0,20 0,056 0,026 
0,25 0,095 0,048 
0,30 0,150 0,082 
0,35 0,225 0,131 
0,40 0,324 0,199 
0,45 0,453 0,291 
0,50 0,625 0,413 
0,55 0,848 0,975 
0,60 1,142 0,788 
0,65 pooo) 1,069 
0,70 2,072 1,444 
0,75 2,834 1,955 
0,80 3,981 2,677 
i | & 950 | 9'608 «| Stora 
—____|___| _j|_. i ,00 9,014 
" Ue un a oe @ 0,95 20,495 9,798 
Abb. 6. Fahrwiderstand und Normalbelastung fiir das gezogene Rad. 1,00 co [e-e) 


Durch die Verbindung von (13b) und (14) kann man den Rollvorgang vollstandig beschreiben. 
Aus (13b) ermittelt man fiir die den gegebenen Parametern N, hy, r, 7, U entsprechende Eingriff- 


ziffer Jj. Daraus ergibt sich aus (14) der Fahrwiderstand Z und aus (12) die Spurziffer 6,. 


An drei verschiedenartigen Beispielen, fiir die wir 1 = 0,01 kg s/em? und b = 10 cm zugrunde 
legen, soll der Rechengang gezeigt werden. 


. h 
1b. perce wird oe = 0,02; U=100 cm/s. Fir N = 300, 400, 500, 600 kg sind Z, 
Onn On ein und die Leistung L zu ermitteln. Man erhalt die Ergebnisse der Tabelle auf S. 235 


oben links. Bei konstanten Werten von h,/2r und U wachsen mit zunehmender Belastung 
Fahrwiderstand und Spurziffer. Auch der relative Schlupf nimmt dabei zu. 
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a Be crmuereselzt wird U = 100 cm/s, N = 500 kg. Fiir ze == = ; = ‘ - ; 7 sind Z, dy, dys 
—q7 > «und L zu ermitteln. Man erhalt folgende Ergebnisse: 

rad 300 400 500 600 = S = ~ fa 

do 0,50 | 0,54 | 0,58 | 0,61 dy 0,70 | 0,58 | 0,51 | 0,47 

do 0,46 | 0,51 | 0,56 | 0,60 a 0,70 ies | 0,43 

eee 0,18 | 0,22 | 0,26 | 0,29 U = ue ae 0,26 0,19 | 0,16 
| e 
‘sa 28,3 | 37,8 | 49,0 | 59,0 ‘ 70,6 | 49,0 | 38,2 | 32,0 
Fs | 0,38 | 0,50 | 0,65 | 0,78 | 0,94 | 0,65 | 0,50 | 0,43 


Bei konstanten Werten von U und N nehmen mit wachsendem Raddurchmesser Fahrwider- 
stand und Spurtiefe ab, insofern die Schichtdicke konstant ist. Bei gleichbleibendem Raddurch- 
messer nehmen mit der Schichtdicke Fahrwiderstand und Spurziffer ebenfalls ab. Der relative 
Schlupf wird dabei geringer. 


3. Vorausgesetzt wird N = 500 kg, =. = 0,02. Fir U =50, 100, 150, 200 cm/s sind Z, 6,, 


Oo 23S und L zu ermitteln. Man erhalt folgende Ergebnisse: 


N=P cos 0 


Abb. 7. Krafte am getriebenen Rad. 


Bei konstanten Werten von N und h,/2r nehmen mit wachsender Fahrgeschwindigkeit Fahr- 
_widerstand und Spurziffer ab. Auch der relative Schlupf wird dabei geringer. 


6. Das getriebene Rad. Wir wollen nun annehmen, daf die gleichférmige Rollbewegung 
durch ein auf das Rad wirkendes Drehmoment aufrechterhalten wird und gleichzeitig eine Zug- 
kraft Z an der Radachse wirksam ist (Abb. 7). 

Aus der Gleichgewichtsbedingung 

Psnax+Z=0 


ergibt sich mittels (9) durch Auflésen nach U*, wenn man wieder 


(2 a — 9) 1004) 


verse pee ae 
ent; are sin Veo = 
U* 2—A Vis, Z hy 
U A A arc sin |/ dy bn UY 2r’ 
Us Ue 2 Vi=6, Z es 
Uy fd A arc sin dy bn UY 2r’ >) 


rail 5 eee V1—6, Z ee 


U A Aare sin | 59 bn UF 2r° 


ee 
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Demnach betragt die Schichthéhe auf der iiberrollten Fahrbahn 


Vi=% 


Z 


= OU UP ht hae (2 
b= Gp ——p o  ed 
so daB wir fiir die Spurziffer 
ga1-E (24 eae 
— 2hy\A ' Aare sin dy bn UY 2r 


erhalten. 


Tabelle 4a. 

0,000 

0,025 —3,010 
0,05 = 1161 
0,10 —0,473 
0,15 —1,320 | —0,153 
0,20 —0,928 | +0,049 
0,25 —0,647 0,197 
0,30 —0,426 0,319 
0,35 —1,572 | —0,240 0,425 
0,40 —1,279 | —0,076 0,525 
0,45 —1,019 | +0,077 0,624 
0,50 —0,778 0,225 Onze 
0,55 —0,547 OSSii7 0,838 
0,60 —0,316 0,538 0,964 
0,65 —0,074 0,717 (iA) 
0,70 —2,016 | +0,192 0,928 1,296 
Ds) —1,557 0,502 1,189 1532 
0,80 —1,031 0,894 1,535 1,856 
0,85 —0,364 1,439 2,040 2,341 
0,90 +0,636 VS 2,893 3,175 
0,95 +2,731 4,323 4,854 5,119 
1,00 le) fore) fone) fexe) 


0,15 1,143 | 0,622 | 0,309 | 0,100 
0,20 0,969 | 0,552 | 0,301 | 0,134 
0,25 0,856 | 0,512 | 0,305 | 0,168 
0,30 0,778 | 0,490 | 0,317 | 0,202 
0,35 1,209 | 0,723 | 0,479 | 0,333 | 0,236 
0,40 1,096 | 0,683 | 0,477 | 0,353 | 0,271 
0,45 1,006 | 0,656 | 0,481 | 0,376 | 0,306 
0,50 0,934 | 0,638 | 0,490 | 0,402 | 0,342 
0,55 0,878 | 0,629 | 0,504 | 0,429 | 0,380 
0,60 0,833 | 0,626 | 0,522 | 0,459 | 0,418 
0,65 0,800 | 0,629 | 0,543 | 0,492 | 0,457 
0,70 1,191 | 0,776 | 0,637 | 0,568 | 0,527 | 0,499 
0,75 1,088 | 0,761 | 0,652 | 0,598 | 0,565 | 0,543 
0,80 1,004 | 0,756 | 0,674 | 0,632 | 0,607 | 0,591 
0,85 0,938 | 0,762 | 0,703 | 0,674 | 0,656 | 0,644 
0,90 0,892 | 0,781 | 0,744 | 0,725 | 0,714 | 0,707 
0,95 0,876 | 0,823 | 0,806 | 0,797 | 0,792 | 0,788 
1,00 1,000 | 1,000 | 1,000 | 1,000 | 1,000 | 1,000 


)=1 


—0,177 
—0,086 | 
—0,007 
+0,065 
0,131 
0,194 
0,255 
0,314 
0,372 
0,430 
0,489 
0,550 
0,615 
0,688 
OG) 


pe a 
—#(G+)) 


1,000 


—0,253 
—0,141 
—0,043 
-+0,046 
0,130 
0,210 
0,286 
0,360 
0,434 
0,508 
0,585 
0,670 
0,771 
1,000 


is 
AS Aare sin 6, 51 U Ta ee 


2,423 

2,330 

2,267 

2,231 | 3,234 
2,223 | 3,146 
2,244 | 3,098 

2,300 | 3,091 

2,400 | 3,136 | 5,344 
2,562 | 3,248 | 5,308 
2,819 | 3,460 | 5,385 
3,242 | 3,843 | 5,646 
4,021 | 4,585 | 6,278 
5,915 | 6,446 | 8,037 
co co 


—0,249 
—0,119 
+0,002 
0,115 
0,222 
0,325 
0,426 
0,527 
0,632 
0,753 
1,000 


—0,193 
—0,002 
+0,178 
0,351 
0,521 
0,700 
1,000 
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Bee 

Mit den in den Tabellen 4a und b oe 
gegebenen Werten sind in Abb.8die 60 

Bildkurven von U*/U und 6o fiir 


die Parameterwerte 


ey ae 
sult (0 0e CEOs 
Sires 0201's. 0 


im Intervall 0 <6) < 1 gezeichnet. 
Daraus geht hervor, daB zwei ver- 
schiedene Falle vorliegen, je nach- 
dem, ob die Kraft Z entgegengesetzt 
(Z> 0) oder gleichgerichtet (Z < 0) 
der Fahrgeschwindigkeit U ist. 

a) Fall Z=0. Im Gegensatz 
zum bereits behandelten gezogenen 
Rad ist fiir das getriebene Rad 
U* = U. Die Spurziffer 6, kann 
niemals negativ sein; denn h>h, 
wurde bedeuten, da durch das 
Uberrollen eine Schichterhéhung 
am Radauslauf auftreten wiirde. 
Durch die physikalisch notwendige 
Forderung 6, =0 wird der Eingriff- 
zifferbereich = demnach  einge- 
schrankt, und zwar um so mehr, 
je gréBer der Parameterwert 


Z i 
bnUY 2r 
ist. Man erhalt somit die Giiltig- 
keitsbereiche 


Mee 6-1 far .7 —0, 
0,225 <6,<1 


fe = 02h e uy, 
0,355 <d,<1 a 

= Ob UV 
0,475 <6, <1 

fir 7 = 1087 ole 


0,600 <6, 1 
fir Z=20b,U|/5", 
0 
0,750 <6, 1 
fir Ze 5,eb7 UV. 
ho 
Die Spurziffer 6, ist immer kleiner 
als die Eingriffziffer 6). Am linken 
Rand jedes Bereiches wird sie gleich 
Null, d.h., das getriebene Rad 


hinterlaBt keine Spur. Am rechten 
Rand 6, =1 wird sie in allen Fallen 


Abb. 8. U*/U und Oo fiir das getriebene Rad. 
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gleich Eins. Dieser Fall kann nur bei unendlich groBer Normalbelastung des Rades erreicht 
werden; er ist also als mechanisch nicht realisierbar auszuschlieBen. 


Weiter erkennen wir aus Abb. 8, daB bei konstanter Fahrgeschwindigkeit U das getriebene | 


Rad nur dann eine bestimmte Zugkraft Z aufbringen kann, wenn seine Umfangsgeschwindigkeit 


U* Werte annimmt, die den Hingriffziffern des zu (Z/b 7 U) Vho/2r gehorigen Giiltigkeitsintervalles _ 
entsprechen. Da die Eingriffziffer wesentlich von der Normalbelastung des Rades abhingt, ist | 


bereits jetzt ersichtlich, daB das getriebene Rad nur dann eine bestimmte Zugkraft ausiiben kann, | 


wenn seine Normalbelastung einen gewissen Kleinstbetrag tiberschreitet. Das dafiir aufzubrin- 
gende U* ist dann, wenn 6, festliegt, aus den (U*/U)-Kurven zu ermitteln. 


b) Fall Z <0. Wenn auf das Triebrad eine Kraft Z<0 in der Bewegungsrichtung wirkt, 
kann die Umfangsgeschwindigkeit U* des Rades bei konstanter Fahrgeschwindigkeit U positive 
und negative Werte annehmen. Da fiir die Spurziffer immer 6,<1 gelten muB, folgt aus (16) 

Uf 
AG: 
Daher wird auch hier der Eingriffzifferbereich fiir den Parameterwert 


LVR 


>—l. 


bnU 
abgegrenzt. Man erhalt die Teens 
Nl So= 1 fin OL u/%, 
0,0600=<6,=1 fir Z=—05byn oe 
0190 0, 1. | fur Ze 0b we 
0,455<6,<1 fir Z=—2,0b7U a 
0.805 <6, = Diet 7 0 be =. 


In diesen Bereichen kann die Spurziffer 6), wie Abb. 8 zeigt, gegebenenfalls gréBer als die Eingriff- 
ziffer dy sein. Dies tritt vor allem dann ein, wenn zur Aufrechterhaltung der Rollbewegung 
U = konst. die Umfangsgeschwindigkeit U* negativ, das Rad also wegen der stark schiebenden 
Kraft Z gebremst werden mu. 

Einen tieferen Einblick in den Rollvorgang des getriebenen Rades erhalten wir erst dann, 
wenn wir die Normalkomponente der resultierenden Druckkraft und das Rollreibungsmoment 
in ihrer Abhangigkeit von 6) und (Z/b1 U)/h)/2r genauer untersuchen. 

Fir die Normalkomponente N = P cos & gilt nach (9) die Beziehung 


N U+U*2r 6 h* U+— U 


Dabei sind 
U + U* U*—U 
Ti und i 


durch (15) zu ersetzen, so daB die linke Seite von der Eingriffziffer 6, und dem Parameter 


(Z/b7 U) Vho/2r abhangig ist. Wie beim gezogenen Rad kann fir hinreichend kleine Werte von 


h)/2r auch hier der zweite Summand unberiicksichtigt bleiben, so da& wir aus 


NOt oem, Use aU gh, h* 
by U2r (hie aes f ap, 2 — 90) 
* 
durch Ersetzen von a y nach (15) schlieBlich 
Naas he 3 dy h* \(3 re Lr ee 
pS eid ee eee | a (Oe pessoa [ae ei A gee | Ye 
bynU 2r  1—6y rh, | 0) rt qieterole) aa 


erhalten. 


ib es 
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Fir das Reibungsmoment (Rollmoment) ergibt sich aus (9) bei Beriicksichtigung von (15) 


M hy _ 9 U*—U aresin Jy Z j= 


brnUY 2r U Vis, ~ bnU V 2r 
und schlieBlich 
Me tn ele A arcan |. 2— A 7 ce (18) 
brn UY 2r A jase AN LpagQel |f Pie 


rae 
onl | ér 


y= 0,80 


-5 


,-G8 


Abb. 9b. 
Abb 9 a,b. Normalbelastung und Reibungsmoment fiir das getriebene Rad. 
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Mit den in den Tabellen 5a und b verzeichneten Werten von 


sind diese in Abb. 9a, b iiber 


Nh 


b 


on 


M 


iy 


bryU ar 


Ingenieur-Archiv 


0,000| 0,000 
—0,001| 0,000} 0,001 
0,000| 0,002) 0,003 
0,003| 0,008} 0,012| 0,015 
—0,005| 0,012| 0,022) 0,029) 0,036 
+.0,002| 0,029] 0,046] 0,057| 0,068 
0,017! 0,058| 0,083] 0,099| 0,116) 0,140 
0,045| 0,102| 0,137] 0,160| 0,183] 0,218 
—0,067| 0,089 0,167] 0,214| 0,245] 0,277] 0,324| 0,402 
—0,048| 0,159} 0,262| 0,324| 0,365] 0,406| 0,468| 0,572 
—0,005| 0,264] 0,398] 0,479| 0,532] 0,586| 0,667| 0,801 
+0,076| 0,422] 0,595} 0,698] 0,767] 0,836| 0,940] 1,113] 1,458 
0,217} 0,660| 0,881] 1,013/ 1,102) 1,190] 1,323] 1,544! 1,987 
0,455| 1,021| 1,305| 1,475| 1,588) 1,701] 1,871| 2,155| 2,722 
0,855| 1,585] 1,950| 2,169] 2,315| 2.461] 2,680| 3,045| 3,776 
1,542} 2,495] 2,971| 3,257| 3,448] 3,638| 3,924| 4,400] 5,353| 8,210 
2,.718| 4,047| 4,681] 5,062| 5,316/ 5,569| 5,950] 6,585| 7,854] 11,661 
5,166| 6,917| 7,792/ 8,317| 8,667| 9,017| 9,542| 10,418| 12,168| 17,420 
10,400} 12,963 | 14,244| 15,013| 15,525) 16,038| 16,807| 18,088] 20,651| 28,340 
24.805 | 29,009| 31,112| 32,373| 33,214| 34,055 | 35,316| 37,418| 41,623| 54,236 
92.010 | 101,181 [105,766 |108,518 |110,352 |112,186 |114,937 |119,523 |128,694 | 156,207 
CO co CO [e.e) [o,e) [o.@) [o.e) [o.2) [o@) [e.e) 
Tabelle 5b. _™@_7/ Fo 
brynUYy 2r 
— 0,2 0,0 0,2 0,5 1,0 2,0 5,0 
—0,200| 0,000 
—0,500 |—0,200| 0,000 
—0,500 |—0,200| 0,001 
—0,499 |—0,197| 0,004 | 0,205 
—1,002 |\—0,495 |—0,191| 0,012 | 0,214 
0,999 | 0,486 |—0,179 | 0,026 | 0,231 
—0,991 |—0,470 |—0,158| 0,050] 0,259 | 0,571 
—0,976 |—0,444 |—0,125| 0,088 | 0.300 | 0,619 
2,039 | 0,948 |_0,402 |—0,075 | 0,143 | 0,361 | 0,689 | 1,234 
2,029 | 0,903 | 0,339 |—0,001 | 0,224 0.449 | 0.787| 1.351 
2,003 | —0,831 |—0,245 |+0,106| 0,340] 0.575 | 0.926] 1.512 
1,949 | —0,721 |—0,107| 0,262| 0,507] 0,753| 1,121] 1.735 | 2,964 
—1,853 | —0,553|-+0,097| 0,487| 0,747 | 1,007] 1.397] 2.047] 3.347 
—1,686 | —0,295| 0,400] 0,817] 1,096 | 1,374] 1.791 | 2:487| 32878 
—1,405 | +0,104] 0,858] 1,311| 1,613 | 1,914] 2.367| 3.121| 4.630 
—0,925) 0,739] 1,571] 2,070] 2,403 | 2,736 | 3,235 | 4,067) 5.731 | 10,723 
—0,084} 1,791] 2,729] 3,292| 3,667| 4,042 | 4,604 | 5.542 | 7.417 | 13.043 
+1,474| 3,651] 4,739] 5,392) 5,828 | 6,263 | 6.916 | 8,005 | 10,182 | 16,713 
4,659| 7,301] 8,621) 9,414] 9,942 | 10,470 | 11,262 | 12,583 | 15,294 | 23.147 
12,494 15,951| 17,680 | 18,717] 19,408 | 20,099 | 21,136 | 22:865 | 26,322 | 36,692 
41,987 47,371) 50,064 | 51,679 | 52,756 | 53,833 | 55,833 | 58,140 | 63,525 | 79.678 
[o.2) CoO [e,2) [e.e) [e,e) [o.@) [e.@) [o.e) CO [o.e) 
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mit der Eingriffziffer 6) als Parameter durch eine Geradenschar dargestellt, die entsprechend der 
angegebenen Giltigkeitsbereiche jeweils durch eine Kurve abgegrenzt wird. Fir 5, = 0, Z> 0 
erhalten wir aus (16) fiir (17) und (18) mit 

Qhe2 io, Z pes 


ee A are sin // 5p byiU 


die Grenzkurven 
1 nee Oo h, 
bn U2r 8 aap? (i 1) +44), 
M Thy __ 2 are sin |/ dy Da (2 Ay 1 
brayU ¥-2r isan h* 


fir 65 =1, Z<0 


INS Jb M 

ec a BEN, fo 4 y/ho 

by U 2r bryaU ¥ 2r ee || Pap 
Ist ein bestimmter Wert von (Z/b 7 U) V(Ap/2 r) vorgegeben, so liefert die Grenzkurve einmal den 
Wert von (N/b7 U) (ho/2r) und zum anderen den Wert von (M/bry U) V(hg/2 r), oberhalb derer eine 
stationare Rollbewegung méglich ist. Fiir eine bei vorgegebener Fahrgeschwindigkeit geforderte 
_ Zugkraft Z wird dadurch ein Grenzwert fiir die Normalbelastung N und das aufzubringende Roll- 
moment M angegeben, dessen Unterschreitung eine stationére Rollbewegung ausschlieBt. Damit 
ist der Zusammenhang zwischen den GréBen Z, hy/2r, N, M, 69, 69 und U*/U ausreichend auf- 
geklart. 

Es ist ohne weiteres ersichtlich, da8 der Fahrwiderstandsparameter 


i [= 
bynU V 2r 
fiir jede Eingriffziffer 6) des gezogenen Rades (Z< 0), fiir das M=0 gilt, durch die Schnittpunkte 


der Geradenschar in der Momentendarstellung von Abb. 9b wiedergegeben wird. Fir M=0 
in (18) wird namlich in st oad mit (14) 


4i—4 arc sin |/ 9 


rs U 2 Se 45 =A WVi—5, * 
und mit diesem Wert folgt aus (15), (16) und (17) 
Us - A Se eo, 
es a a a? 
N hy 6 ie lye bE Hie =a9): 


Das sind aber gerade die Werte (11), (12) und (13b), die wir fiir das gezogene Rad erhalten haben. 
In Abb. 8 und 9 ist dieser Sonderfall durch gestrichelte Kurven gekennzeichnet. 

Braucht das getriebene Rad keine Zugkraft aufzubringen (freies Triebrad), ist also Z = 0, 
so wird 


Cee ae oe ee Ue ee 


vic ey ba: U As U As 
und iy 
Oo — 1 Sh Al ° 
Fir die Normalkomponente N= P cos « gilt dann die Bezienung 
No ah 6 5, h* 
bynU P=70-5|! an, 2 Oo)) » 


fiir das Reibungsmoment e. o 
M hy ii 1—A arcsinJ/ dy 


brnUVy2r Ao, 
Bei der numerischen Auswertung berechnet man zuerst 
Ties Pie 
by U 
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und ermittelt fiir eine gegebene, der erdrterten Beschrankung geniigende Belastung N aus 
N hy 
bn U 2r 
an Hand von Abb. 9a, b die Eingriffziffer 6). Darauf bestimmt man aus 


_M Ve 
brnU ¥ 2r 


das Reibungsmoment M. Die Gleichungen (15) und (16) liefern dann die Werte U*/U und dy: | 
Den Beispielen legen wir wieder die Werte 7 = 0,01 kg s/em? und 6 = 10 cm zugrunde und | 


schlieBen sie teilweise an die des gezogenen Rades an. 


1. Vorausgesetzt wird h,/2r = 0,02, U = 100 cm/s. Fiir die Belastung N = 300 kg des gezogenen 


Rades hatten wir die Zugkraft Z = 28,3 kg errechnet. Mit den Werten 
Z y/ho 


= Cn Te es bn U Peis 


erhalten wir folgende Ergebnisse: 
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ue 120 225 | 450 920 | 1300 $M und der Leistung L angegeben, oberhalb derer 
eine stationdare Rollbewegung méglich ist. Es zeigt | 
On 0,32 | 0,40} 0,50| 0,60| 0,65 sich also, daB eine Belastung des Triebrades von 
mindestens 120 kg erforderlich ist, wenn eine Zug- 
do 0,00 | 0,11) 0,22) 0,34) 0,39 kraft von 28,3kg aufgebracht werden soll. Fir 
diesen Grenzwert der Belastung wird die gesamte 
Diisess Uf 0,64 | 0,61 | 0,63 0,73 | 0,83 Schicht tiberrollt, ohne daB Teile von ihr seitlich 
U abgedrangt werden. Die Spurziffer ist fiir diesen 
M Grenzfall gleich Null. Die Leistung und die Spur- 
kgm OPE ERS iOS a2 Lo ae Clam itter eas mit der Belastung et Der solute 
T Schlupf (U*—U)/U hat bei der angenommenen 
PS 0,83 | 1,02) 1,54) 2,69) 3,79  seringen Zahigkeit der Schicht verhaltnismaBig hohe 


In der ersten Spalte sind die Minimalwerte von N, 


Werte. Wir werden im dritten Beispiel den EinfluB 
gréBerer Zahigkeitswerte 7 baw. gréBerer Fahrgeschwindigkeit U auf die Rollbewegung erlautern. 


ee 1,00 0,69 0,47 0,31 
2. Vorausgesetzt wird 5p 0,50 | 0,45) 0,40) 0,35 
U = 100 cm/s, N = 450 kg, Z = 28,3 kg und r = 25 cm. . 
Wir erhalten fiir die Schichthdhen h, = 1.00, 0.69, 0.47, ee Daren cs 
0.31 cm folgende Ergebnisse: VU 
ae 0,63 | 0,54 0,46 0,39 
U 
M 172689) 05-60 S13 75nd 20 
Bei konstantem Raddurchmesser nehmen mit der Schicht- kgm ‘ ; : ‘ 
dicke Schlupf und Leistung ab. 7 
PS 1,54 1,28 1,07 0,90 
U 100 150 225 340 535 
3. Vorausgesetzt wird ee ee 
N= A50lke, 226.3 koe hie seem 80] 705905) bn 0,450; A0N SSS eam 
Damit erhalten wir fir U = 100, 150, 225, 340 
9395 cm/s folgende Ergebnisse: °o sans waa eee fey tene 
U*—U 
7 OO) |) Oza I SYN | OhOls 0,18 
M | 17,68 | 17,3 
Mit zunehmender Fahrgeschwindigkeit U oder zu- kgm ; 27 nL 20g eee 
nehmender Zihigkeit 7 nehmen Eingriffziffer und L . 
Schlupf ab, wahrend die Leistung ansteigt. PS 1,54 | 2,04 | 2,72) 3,75 | 5,39 
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4. Vorausgesetzt wird 


» 300 480 735 

Z=0, hyo =1lem, r = 25cm, U = 600 cm/s. : 
Fir N = 300, 480, 735 kg erhalten wir folgende Ergebnisse: Oo 0,200 080 70,35 
do 0,17 | 0,20 | 0,24 
eo. 64 || 0,06. |. 0,09 

U 
Das freie Triebrad hinterlaBt bei gréBerer Belastung auch eine sa 4,45 | 6,68 | 9,65 

tiefere Spur. Der Schlupf ist gering. Schlupf und Leistung aed P 

nehmen mit der Belastung zu. ae 1,48 | 2,27 | 3,37 


Herr Dipl.-Math. W. Miiller hat die numerische und zeichnerische Auswertung durchgefiihrt. 
Ihm spreche ich an dieser Stelle meinen Dank aus. 


(Eingegangen am 15. Oktober 1956) 
Anschrift des Verfassers: Professor Dr.-Ing. A. Kneschke, Freiberg (Sachsen), Richard-Wagner-Str. 13. 
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Ein weiterer Beitrag zur Theorie der Bewegung rollender Fahrzeuge 


Von B. Stiickler 


1. Einleitung. Die vorliegende Arbeit stellt eine Erweiterung und Erganzung der Arbeiten [1 
und II? dar. Vor allem werden die bei einer beliebigen ebenen Bewegung eines rollenden Kraft- 
wagens an den Radern wirkenden Haftreibungen ermittelt, und zwar in Abhangigkeit von belie- 
bigen auBeren und inneren eingepragten Kraften und den Massenbeschleunigungen, so weit es 
unter den zugrundegelegten Voraussetzungen tiberhaupt méglich ist®. Dabei wird der Antrieb 
vom iibrigen Fahrzeug gesondert behandelt (wie in I). Dieses unterscheidet sich von dem in II 
behandelten Wagen (Fall A = 0) hinsichtlich der Kinematik nur im Vorderradsystem. Es brauchen 
also nur noch diejenigen Massenbeschleunigungsglieder bestimmt zu werden, die bei 2 + 0 zu 
den in II abgeleiteten Beziehungen hinzukommen. Obwohl diese Glieder keine beachtliche Rolle 
spielen, solange die entsprechenden Tragheitskrafte der Vorderrader verhaltnismaBig klein sind, 
wird die Rechnung zunachst allgemein ohne weitere Vernachlassigungen durchgefiihrt, um die 

Leistungsfahigkeit des in II entwickel- 
ten Verfahrens zur Bestimmung von 

Reaktionskraften bei komplizierten 

Systemen zu demonstrieren. In Ziffer 6 

werden dann aber alle Glieder ge- 

strichen, die bei einem normalen Kraft- 
wagen entweder klein im Verhaltnis 
zu immer gleichzeitig auftretenden 
andren, oder die auch im ungiinstig- 
sten Falle nur Tragheitskrafte von 
einigen Kilogramm ergeben wiirden. 

Die ganze Rechnung wird mit einer 
kleinen Modifikation gegeniiber I durch- 
gefiihrt. Es hat sich namlich als zweck- 

maBiger erwiesen, den Bezugspunkt B 

fiir die Geschwindigkeit v von der 

Lenkachse des linken Vorderrades* 

genau zwischen beide Lenkachsen zu 
az  verlegen (Abb. 1). Das scheint aller- 

Abb. 1. Bezeichnungen am Vorderradsystem. dings dén Nachteil zu haben, daB die 

in I angegebenen Rechnungen zur 

Aufstellung der Differentialgleichungen I (34) fiir die Bewegung der Fahrzeuge wiederholt werden 

miuften. Da aber aus ihnen unschwer zu ersehen ist, welche Glieder auf das linke und welche auf 

das rechte Vorderrad entfallen, lassen sie sich ohne Miihe auf die neu gewahlte Lage von B 

umschreiben. Lediglich zur Kontrolle wurde die Rechnung unabhangig von I mit Hilfe des in II 
angegebenen Verfahrens noch einmal durchgefiihrt. 


2. Zur Kinematik der Vorderrader. Da sich von den in II gewahlten Verschiebungen }s 
nur Ds; und 0s, auf Punkte der Vorderrader beziehen, kénnen nur diese im Falle 2 +: 0 in andere 
Weise von den wahren Koordinaten q, abhangen als in II. Genau so steht es mit den Geschwin- 
digkeitsgréBen «,, so daB fiir diese mit Ausnahme von ws und @, die Beziehungen II (21) bis 
II (24) auch bei A + 0 gelten. Infolgedessen kénnen sich in der Koeffizienten-Matrix (b,,) II (26) 
nur die dritte und vierte Zeile andern und in der reziproken Matrix (B,) ALT) Ane ahs fiir p 
baw. @. geltende fiinfte bzw. sechste Zeile. Diese Anderungen sind zundchst zu bestimmen. 

Bei einer ganz beliebigen, d. h. auch das Gleiten einschlieBenden Bewegung des Fahrzeuges 
auf der Fahrebene, gilt fiir den Geschwindigkeitsvektor vg der Lenkachse des rechten Vorder- 


1 B. Stiickler, Ing.-Archiv 20 (1952) S. 337. 

2 B. Stiickler, Ing.-Archiv 23 (1955) S. 279. 

3 Siehe I, Vorbemerkungen. 

4 Siehe I, Abb. 1. 

® In If (27) muB es (B,,) statt (B,,) heiBen. 
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rades, wie leicht aus Abb. 1 zu ersehen ist, 
a Ve = 0g + 8, xb’ (1) 
mit b’ als Vektor von B nach C. Daraus ergeben sich die Komponenten 
v, =x cos 0, + y sin J, + B’ dy C08 Yo ; | 
v, = —xsin 0, + ¥ cos 0, —b' d, sin y, .{ 
Nach I (8) ist aber 


1 : : 5 : 
oa, (ook Aang) g(A.xy=+yl +A sin 2y + A? sin? y , sin yy =~ sin 7, (3) 
so daB aus (2) wird 


1 : ‘ 
v, = | [(cos 8 + A cos #, sin x) x + (sin 0 +A sin J, sin x) y +12 (cos y +A sin x) i], 


4 
= = (sin } + A sin J, sin y) x + (cos # +A cos J, sin y) y —lA sin 0]. o 


Die Umrechnung auf die GeschwindigkeitsgréBen w, ergibt mit Hilfe der drei ersten Zeilen von 


II (27) 
1 =8'G)—* (cosy +Asing)or+eo,, sat, (5) 
ee g(r + Asin y Ws.) . (6) 
Nun gilt fiir die Geschwindigkeit @, des Beriihrungspunktes mit der Fahrebene (siehe Abb. 1) 
o,=,+¢8,—rg, 8, =0, +%; (7) 
also nach (4) und I (10f) 
wo, =< (cos 6 +A cos #, sin y) # + + (sin 4 Asin 9, sin x) ¥ 
se e+, (cosy +Asing)| by + ei —rbe- (8) 
Hieraus ergibt sich die entsprechende Geschwindigkeit w; am linken Rade einfach dadurch, daB 
4 durch — A ersetzt wird, ¢ durch — c, @, durch @ und g durch 
g=g(—A,y%) = +yl—Asin2 7 +7? sin? x. 


Somit sind in II (26) die dritte und vierte Zeile zu ersetzen durch 


1 


z (sin  —/ sin #, sin x); 


(b3,) = lz (cos } —/ cos J, sin x); 
ee ae — ri 020: o}, 

g g 
(sin } + Asin J, sin x); 


oe | 


(i ye ie (cos } + A cos J, sin x); 
c Sey + 4 sin x); a 0: ——=7 (0; 0}. 
3 s | J 
Dadurch ergeben sich fiir II (27) folgende Anderungen: 
c g’ c A : I ee Ce ; c 
(B; ,) =|— Fons + ny ad a sin ¥) 5 eng 05-0; 0; (én) : — ap 


/ 


c gs. 
(Bey) = jz 0082 + — 


Le 


A 2 1 1 a ‘ 
rye Lia serra +(8 + sing); ral 
1 
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Damit kénnen die f,,,, berechnet werden. Die Anderungen gegeniiber II (28) lauten: 


Lae, Ase i} 7k 
ame Pot = pg he Piss = a> Pass = — 7a 8D Y > 
1 (11) 
sy 


ne 
Be’ ET ae te 


Pia = rgsin g Boa = go %> Bras = 


Hinzu kommen fiir die Bewegungsgleichungen 


Psr7 = 1, Bug =—1 


3. Die kinetische Energie der Vorderrader. Auch hier braucht man die Rechnung nur fiir ein | 
Rad durchzufiihren, denn der Anteil des anderen erhalt man aus dem Ergebnis durch einfache | 


Vorzeichenanderungen. Fiir das rechte Rad gilt zunachst (siehe Abb. 1) 


2 Ta, = my vt + AB + CGH 42 mpig (I,Xx2), (12) 


mit e als Vektor von C zum Schwerpunkt S des Rades. Ferner gilt, wenn gleich quadratische 
Glieder mit zwei beim reinen Rollen verschwindenden Oy fortgelassen werden, 


ve sot + U2 + 25,(0, xB), 
85(0, xB!) = —W A, (@, sin y— 07, €08 7) (13) 
are en | 


Dieses und (5) in (12) eingesetzt liefert 


2 Tra = Maoh + 208 + CHE L(A + my) HE +2 dy [mg (ee! —¥ sin) 0 


A é , A 
Mp e-- (cosy +A sin YX) og + imp (eg + 6 cos x) a, +20, (14) 


+2mp a gf 0, —~ (00s x +A sin x) Os + B07] 09. 


Hierin sind noch b, gemaB der dritten Zeile von II (27) und @, nach (10) durch die w, auszudriicken. 
Mit der Abkiirzung 


(15) 
ergibt sich, wenn Glieder gestrichen werden, die zusammen mit den entsprechenden des linken 
Rades Null ergeben, 


Tre = 308m’ + (B +m! 22) sint x +m" sin 7 8| + 330% | 


2 


+ olor} a (js tm 2B) sin 2 +5 m’ ‘(sin z8))] +6 [Fa 


cosy +c mea 


a | 


+ Wy!@_ (e + as alae — ne (16) 


+ Ws le? s+ cho E (cos x4 +Asin | 


G Clam Cc A’ Ue 
4-3 (c i Bee ae 5 sin y 4 ‘= | 


4. Ermittlung der Reaktionsgréfen 4, und der Bewegungsgleichungen. Zunichst sind die 
Anteile des rechten Rades an den Impulskomponenten P,, durch Bildung von 0TR,/0w, zu be- 
stimmen, woraus sich dann die Anteile des linken Rades Geieder einfach durch Piseted von A, 
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und g durch — 7, — c und g ergeben. Die Rechnung liefert an Stelle von II (32) 


P= Oe [ra eos 24+ ie + 72 m') sin 2 y + - m’’ [sin y (g — ar ) 
A’ 1 1 ey ’ 
+097 0087 (z+ ge) tem oak 
P=", Ime cos ¥ + (mg + 22? m’) sin y 
SP ot & 2—Asn2y 2-+Asin2y) 
AS me 2a 2% a 2¢ a 


A’/1 GN md 1 . 
= 067 (gat gz) beam a (008 2 + A sin x) ayee 06 1 tem) 


| 1 ay 


C |/- : 

Pra S|é—fani)or—for) 
C 

Py = = (c : sin i) @, + a5 , 


AS .4 il I 1 1 
+o sin 7 |= eed = 
8 1 ye ao fa T g E 


ae fi, 98 ot 5 AS Ph We 
Paso | Tsing (gat ga) tem F z) 4'(a I ee } 
Durch Einsetzen von II (28) in II (14) bekommt man unter Beachtung von (11) anstatt IT (33) 
. 1m 1% : 1 IP 
Ay =— Ky + Py +9" Py 05 x Boe aa a + Pa)sing + (ptt Pies 
G WB P : : 
A, = — K, + P, a iP + P, cos y + ak sin? y — (P; + P,) sin | 
: P db 
—Asing (3B — 3) Ws > (18) 
is an aa K, =i Ps ? 
As = K, ae i ? 
As a K; ae P; ? 
ee tk oP, ) 


und fiir die Bewegungsgleichungen nach [I (11) 
K,= P, — P, 0, 


: oT, (19) 
Ko =.P; ay Py Ons 
Fiir die eckigen Klammern in 4, und /, ergibt die Rechnung jetzt 
Is : 1 1 ; : 1 1 
Ay :|mp = | 2 mp 4 map sing z))sinze +mnesinz (— ss) 2 | 
my Is* , +m age 1 at: ° 
=m sin 74 v R sin (3 a)ie 
[ pare ee ie { (20) 
Age my (17) e082 + My — sin X + Mp pp sin' x Et eye 
C 1 . I5* “i * i A 
— mp eA sin? 7 iat a =—m{(l —T) eos +=; sin x sin YU 
asint (a + 3)i ) 
mat LT IAG Ge 


Auch hier lat sich zeigen, daB die Glieder mit v nichts anderes als die Komponenten des Impulses 


m vs+ des ganzen Fahrzeuges in @,- bzw. w,-Richtung sind, 
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Mit (17), (18) und (20) bekommt man an Stelle von II (36) 


4,=—K, + {mee cos 24 Her + A? m' sin 2 ¥ +m’ [sin x (g —§)y'l 6 


+ mip sing & + [7 cos (a + alee” (gale 


: 1 1 ” ea U 7 
+2 mai sing (3s — ga) +2 : m cox —8)] vi 


A 1 NG 1 1 
+5 cosa (aa + Bl +m ¢(se— 35] 


| 


+2cm’? sin? (2 cos 2 y — A? sin? y) (as 2) — cos 22(35— alle > 


g g 
fy = — Ky — |p 08 7 + (mg + 272m’) sin x + pm" (g — 2) 
2—Asn2y 2+Asin2y]). m Is* se. : 
+ 2g — 2g et Fb —Teosa +7 sing] sing 


g e 
1 


— |(ms + 228m!) 008 y— Mp, sing + A sing (e— F 


Al i Tied ae dese CTs 
re | +Acm cosy (a +g) t Bem sing (a— gi] x 


j 


(21) 


5 if 1 a ; F 
Sir qa sintx (3 | a) PAT (2 — 3 sin? y + / sin 2 ¥) 


+ AT @—3sint yz —dsin2 | 04 


PAGAN Ths ae 1 1 3 1 
4 a r(ge A) + 2A2cem cos x (Ga — a) + ma ed sin (Ze + 


+ BR em sind x 008, (Ga — Fa) — 219 om’ sin? 4 (5+ 

Gh ie CRS. Came 

A, = R,, = — K, aa (a ; sinz)o — ii), 
Cd Car Cur 

A= Ry =— Koa ale + poing)o+ Si, 
C, d ° 

dg= By =—Ky— FG (008% —T sin), 

Ag Ry, =—Ky— FG (co +7 sing) 9. 


1 


e 


, 


J 


Im Gegensatz zu II (36) gilt aber im Falle 2 + 0 i. a. nicht mehr einfach A, = Ry und A, = Ry, 
denn nach dem in II iiber die ReaktionsgriBen /,, Gesagtem gilt zunachst fiir die seitlichen Haft- 


reibungskomponenten (siehe Abb. 1) 
Ay @, +A, @, = Ry, v, + Ryd, + Ry @,. 
Daraus wird nach (6) und wegen 0, = v,(— 2) 


A, ®, +1, 0, = Ge + AA. @, + Ru +A sin 7( 


& 


woraus wegen der Unabhangigkeit der w, folgt 


dy = AE 4 AY, 
Ay = Ry +A sing (“y"— a 


Ryr 


(22) 


(23) 
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Da8B trotzdem in der Praxis nahezu 


A, = Ry, + Ry, hy = Ry 
ist, wird sich in Ziff. 6 zeigen. 


SchlieBlich ergeben (17) und (19) die zum System (21) gehérenden Bewegungsgleichungen 
ko m, + My, sin2 7% + (m, + 2 1? m’) sin? y + i MO Nas 0 


Mee et Pee. VA oe 
+ [renee t altem'(G—F) i 


laces 2 Mis gong Sk 1 Fron ies a (24) 
\nre a+(3 +A m’) sin 2 x 4 ym’ [sin x (g RE. 

+ [tamale pl tem GE 9) | 
Ky = [Fsing(s + a) tem’(G—allet+4'(G ae 

+ffonx(b pen (Ss i 


pany It 1 1 ; 1 1 : 
+A A 2cos 27 (s4— ge) —Asin2x (is + go] 2 J 
5. Der Antrieb. Die vorstehenden Ergebnisse sind nun noch zu erganzen durch diejenigen 
Massenbeschleunigungen, die bei voll eingekuppeltem Motor infolge Drehung von Motor- und 
Kraftibertragungsteilen auftreten. Unter den in I, S. 344, angegebenen Voraussetzungen gilt 
fiir die Rotationsenergie T, des Antriebes zunachst 


2 T= DIV +d Ines +d vedo (26) 
v= v=1 Vi 
Daraus wird nach I (38) und I (43) 
27, = Jy (P+ 2)? + Ja Pi + 2)? + 2 Iu (P + 2) (Pr + Ps) « (27) 
Fiir die Impulskomponenten P,,, folgt 


On 5 é 5 Om. $ 
Lr I, = [Jy (PY + %2) + Iu(Gi + Pala ++ Pz) 
(28) 


2 - F 5 OD) pc : 
+ [Ju 1 + Ps) + Suc (P + $2) 1 5,,- (Pi + Ps) - 


On 


Da nach der Differentiation nach w, die fiir das reine Rollen geltenden Bedingungen w, = 0 
(u = 1, 2,...6) eingesetzt werden diirfen, folgt aus (10) 


r@ +p) = + Bote (a ze (29) 
und aus der siebenten_und achten Zeile von II (27) 
1 (ji +s) = cosy +rsiny)o,, t= >". (30) 
Damit erhalt man aus (28), wenn noch die Abkiirzungen 


; 31 
D=V(g+2)+M (2cosy+7sinx), (G)) 


E=M (g +8) + H(2 cos xy + 7 sin x) 


eingefiihrt werden, 
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AN, ; il 1 Pt 
Pi, =(D(g +8) +£EC cos —2sing)] 0, + ¢(2s — [Vig+eg)y } 
+ M (r cos y — 2 sin y)] wg , 


ioe =—|Pheosx(2— 4) + D® sing (+ z) ty Wy 


1 1 i iy : 1 1 
—e(a— a) (V4 coox(e— 3] +V2?? sing (_ +. ae MT| gz, 
P,, = P,,=— Do, + Vela g) oe (32) 
1 1 
Po3= P44. =— Be; 4 Mela 2) a» 


P,, = [D(g +8) + E(2 cosy +7 sin x)] o, 
alas a) +g) + M (2 cosy +7 sin x)] wg, 


1 1 ul 1 \2 
Py = De(zs— a) . aa 7) Wg- 


Diese Impulskomponenten treten bei voll eingekuppeltem Motor auf den rechten Seiten von (17) 
additiv hinzu. Es zeigt sich, daB der Antrieb keinen Beitrag zu den eckigen Klammern von (18) 
liefert, so daB er nur zu folgenden Zusatzgliedern auf den rechten Seiten von (18) und (19) fiihrt: 
r Pa Pa } 

hao rly at Ps) Os» 


has = has = Rae (33) 


kas = ee 


Hieraus erhalt man durch Einsetzen von (31) und (32) schlieBlich fiir den Fall t = 0, d.h. fiir 
a = b, was beim Kraftwagen wohl stets gilt, 


Ajy={— 2 Hsin2y + V(g +8) (g +8)’ + 2M [cos x (g + 8)’ — sing (g + g)]} 0 
+e(a— gM +8) —2 Maing) 
+ [4H sin? y + V (g + g)*—4M sin x (g + g)'] 07x 
+e(— pa) (Ve +8! —2 Min gl ft 


Aag= loos (¢ — g) t4sina(; + 3)|4ea> 
i (34) 
ha Aaa =—(V (6 +8) +2 Meosy]i—Vo(ss— 2] 
— @ +8) —2 Msinz] 0 —Ve(e— 3.) i? 
fag = hog = — [2 Hoos x + M(g + B)16—Me(— i.) 7 
—[-2Hsing +MQ+8)]94—Me(a— 3) 2, | 
br [A H cost + V (G+ 8% + 4M cosy (+ B)16-+o(2— Z)IV (+8) +2Meosy]j 
+ [V(g +8) (¢ +8) —2Hsin2 y — 2 Msin x (g + g) +2 Mcosy(g+4 g)'] vx (35) 


1 JEN = 5 
toh) tet D429 Mcoezis 
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1 i y 
bee = ¢( 3s — ave + g) +2 M cos x] v + Vela ali 
; . a (36) 
FIG +8) —2 Min hog + Ve(2— LY #1. 


Am einfachsten steht es beim Hinterrad-Antrieb, fiir den V und M durch Null zu ersetzen sind: 


Jay = —2Hsin2 yi + 4H sinty vy =—4Hsiny 5 (v cos 7), 
Age =4as =A,,=90, 
das Aye =—2 Hoos yi + 2H siny vy =— 2H“ (vcosy), (37) 


k,7 =4H cos? y¥v —2 Hsin2yvy = 4 H eos 7“ (v cos x), 
kg OF J 
Nun ist aber v cos y nichts anderes als die Geschwindigkeit v, des auf der Langsachse liegenden 
Punktes der Hinterradachse, so daB einfach gilt 
Agi = — 4H siny 0, , 
Ags = 458 = — 2H 0,, (38) 
kg =4H cos 70, . 


6. Naherungsformeln fiir Kraftwagen iiblicher Bauart. Bei solchen Wagen ist 4 nie viel 
gréBer als 1/4 und y,,,, ~ 30°. Daher ist 


g=ylt+Asin2y +7? sin? y ~ 1+ 4sin2z, 


Ps 7 nae 

g=g(—A) zx 1— , sn2y, (39) 
L7m 

Bae? 


wobei der gréBte Fehler etwa 0,2% betragt (bei A = 1/4 und y = 30°). Es kann also naherungs- 
weise gesetzt werden 


Es 1 1 ae 
81rgE=2, Fea ee (g +g) =0, 
§ & 
-) ee 1 1 asin 2 1 1\/ (40) 
cinema tii oner Far <) ie 


Damit ergeben sich schon erhebliche Voreinfachungen. Streicht man auBerdem alle Glieder, die 
entweder klein sind im Verhaltnis zu immer gleichzeitig auftretenden anderen oder die — von der 
Momentengleichung des Vorderradsystems (~ = 8) abgesehen — auch im ungiinstigsten Falle 
nur Tragheitskrafte von wenigen Kilogramm ergeben, so wird man (jedenfalls solange die Ver- 
haltnisse nicht wesentlich anders liegen als beim zugrundegelegten Zahlenbeispiel aus I, Seite 348) 
mit folgenden Naherungsformeln auskommen: 


} 


LL) 


dy —K, + |m T0082 + "2 sin 24 —4 (H cos x + M) sin z) 6+ mp sey 


| im Ms m7 sin 2 — (me — 4H) sin? 7} 7, 
se are m Ls eeSiee es ; 3 
fy Re —Ky—(m “F008 4, + msing)é + 7 |(1—F)eos x + rsing|sin 


* 


— (macosx —m*> sin) 0 


Ang = Ry & — Kz (7 eat F2 Moos z)i + (pheos2z +2 Msin yz) vj, (41) 


Awae=R ww —Ky— (a +2V + 2M cos x) +(—jrheos2y +2Msin| 07. | 
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hag Ry & — Ke — [2008 4 —F sing) +2 Heong +2 Mi 
1 
ae [P(e0s + sin x) + 2H sin 2] v7, 
ZL 
Ago Ry, & — Ke — | (008 x + “sin z) +2 Hoos z+ 2M)? 
ab 


as Fi (— cosa + sing) +2 Hsing] 07, 
1 
K, [my +m sin 2 bmg sind y +4 (V + Heosty +2 M eos )) 6 | 
(42) 


+ [m5 005 2 + “GPsin 2. — 4 sin z (H 00s x + M) OX | 


47 


Keon 2 Fsiny—Aedsin 2a (+ V+ Meosy)|b+ CA +4 VAM sin2 2) 


+ 2 cosy 2m eh cos2y +4 Medsiny sin2 y] v7 


+ 2Asin2 7 [A’A (6 cos2 ¥—1) +4 Vc? A cos 2 4] x. ; 
Diese Formeln gelten, wenn der Motor voll eingekuppelt ist; beim Frontantrieb sind H und M, | 
beim Hinterradantrieb V und M durch Null zu ersetzen. Bei ausgekuppeltem oder nicht ganz — 
eingekuppeltem Motor kénnen H, V und M niaherungsweise durch Null ersetzt werden’. Be-_ 
merkenswert an (41) ist, daB in A, die Massenbeschleunigungen des Antriebes nicht mehr vor-— 
kommen und in A, nur bei 4-Rad- oder Hinterrad-Antrieb. 


Setzt man ferner (39) in (23) ein, so ergibt sich zunachst 


Ae 
A, © Ry, + Ry, + sin2 x (Ry — Ry,) » 


eee (44) 
A, % Ry —A sin x (Ry;— Ry,) —Z sin sin 2 y (Ry, + R,,,) 
Nun ist aber fiir 2 = 1/4 und y,,,,, = 30° 
ATs : AP 2 
> sin2y < 0,11, Asin y < 0,125, > siny sn 2 7 < 0,014, 
so da meistens sogar mit 
A=Ryt+ Ry = Rp, A=Ry (45) 


gerechnet werden kann. Dann kénnen aber die Naherungsformeln fiir A, und A, in (41) auch als 
Naherungsformeln fiir Ry und Ry genommen werden. Diese unterscheiden sich von den fiir 
A = 0 genau geltenden Formeln fiir Ry und Ry in II (36) nur dadurch, da die y-Glieder beim 
Kraftwagen wegen Belanglosigkeit fehlen, in my, a = b gesetzt ist und in Ry bei voll eingekup- 
peltem Motor noch Massenbeschleunigungen des Antriebes in Erscheinung treten. Somit gilt II, 
Ziff. 5, auch fiir den Kraftwagen; von den im Zahlenbeispiel untersuchten Gliedern andert sich 
bei voll eingekuppeltem Motor lediglich (R,);, und zwar nach der ersten Gleichung von (41) um 
das Zusatzglied —4sin y (H cos y + M). Hierfiir bekommt man mit den Zahlenwerten aus I, 
S. 348, fiir den Hinterrad-Antrieb 


— 5,1 sin 2 im direkten Gang, 

— 103 sin 2 y im ersten Gang, 
fiir den 4-Rad-Antrieb 

— 2,9 sin y (1 + cos y) im direkten Gang, 

— dlsiny(1-+ cosy) im ersten Gang. 


Kin Vergleich mit II, Ziff. 5, zeigt, daB die Massenbeschleunigungen des Antriebes im ersten Gang 
vor allem beim Hinterrad-Antrieb einen erheblich gréBeren Anteil an (R,); haben als diejenigen 


1 Siehe I, 2h), letzter Absatz, S. 346. 
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des ganzen iibrigen Fahrzeuges. Allerdings ist zu beachten, da8 man iiber die Abhangigkeit der 
erforderlichen Haftreibungskraft Ry vom Antrieb erst dann ein richtiges Bild erhalten kann, 
wenn man auch den auf den Antrieb entfallenden Anteil der eingepragten KraftgréBe K, kennt 


(siehe hierzu Ziff. 7). 


SchlieBlich sei noch als Ergainzung zu II, Ziff. 5, kurz auf den Einflu8 der Schwerpunktlage 
auf die vy-Glieder in Ry und Ry eingegangen. Diese spielen vor allem beim Hineinfahren in eine 
Kurve und beim Herausfahren aus ihr eine Rolle. Benutzt man fiir die betreffenden Koeffizienten 

‘in (41) die Naherungsformeln II (46), so ergeben sich die Diagramme der Abb. 2. Aus diesen ist 
zu ersehen, daf} — jedenfalls fiir das gewahlte Zahlenbeispiel — im praktisch interessierenden 
Bereich 1/4 <1,+/l1 < 3/4 die erforderliche Seitenkraft Ry sehr betrachtlich von der GroBe des vx- 
Gliedes abhangen kann; dieses wiirde beispielsweise fiir 1,«/1=1/2 und vy = 5 m/sec? etwa 350 kg 
ergeben. Der Antrieb tritt nur in (R,),; in Erscheinung, and auch dann nur beim Hinterrad- 
Antrieb und praktisch nur im ersten Gang bei gréBeren y-Werten. 


(Ri), a 
200 3 SF 
kgri"| sek? eet Nei kg 
aT AIT TIEL) 
== Lsx/li 
150 eS 1 
= ao | 
3 
¢ 
A. SS SSS eds 
¢ 
50 
————_— — 4 Gang 
| 
reel ; direkter Gang : x 
0 15° 30° 45° 0 
pee ee ee ee ee eee 
100 52 37m 


Abb. 2, 


7. Uber die K,, der wichtigsten eingepragten Krafte. Auch bei A + 0 gilt fiir die Schwerkraft 
II (40) und fiir den Luftwiderstand II (41), dagegen ergeben auf die Vorderrader wirkende Momente 
um deren Achsen an Stelle von II (43) wegen (10) folgende Naherungsformeln 


A 1 ; 
1 ~ Tr Cr: caer M,,)) cos 2 At t= a (M,,, —= M,,,) i oe (M,, a M,,1) sim Y ; 
nN : 
K, % — = (M,, — M,,)) cos x ey Mi en) 
M, M,, Mni oa Mnr 
K,=——, K,=——", IG aareeee Kees 
1 i 1 (4) 
K, ben (M,, = M,,) =e or Vine M,,) sin 234 71 = (M,,, = M,,) cos ¥ 


= a (M,,—M,)) sinx, 
Kg ~ < [M,, aa M,, 7) Ce ae M,,) sin 2 vA ° 
Fiir den Antrieb gilt zunachst nach I (51) 


DA, ae M,, (Dp = DPp) =F Mn (Dp, is Ds) y (47) 


Dasselbe bekommt man aber auch rein formal, wenn man in IT (42) M,, und M,, durch M,, und 
M,, und M,, durch M,,, ersetzt. Daher ergeben sich die Naherungsformeln des Antriebes einfach 
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aus (46) durch die gleiche Substitution: 


Ki 
Ky, 
Ky3 


K,; 
Ky 


Kgs 


fe ae Pitt aes ) 
a (Dc 

=K,.=—"**, 
ee 

x 2=* + 271" 08 x, 
a 2" eA sin 2 x. 


(Eingegangen am 15. Oktober 1956) 


(48) 


-Anschrift des Verfassers: Dr.-Ing. B. Stiickler, Berlin-Steglitz, Mariendorfer Str, 32b, 
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EinfluBflachen einer Dreiecksplatte mit Aufpunkt am freien Rand * 
Von A. Scherer 


1. Einleitung. Im folgenden sollen singulare Lésungen der biharmonischen Gleichung her- 
geleitet und das Verfahren zur Ermittlung von EinfluBflachen von F’. Schultz-Grunow! weiter ent- 
wickelt werden auf Platten mit teilweise freien und eingespannten Randern fiir eine am freien Rand 
wirkende Last. 

Es sei vorausgesetzt, daB der Elastizitatsmodul E und die Querkontraktionszahl » konstant 
seien, daB die Plattendicke A konstant und klein gegen die iibrigen Abmessungen sei, der Baustoff 
homogen und isotrop, die Durchbiegung klein gegeniiber der Plattendicke und ferner, da das 
Hookesche Gesetz giiltig sei. Die Differentialgleichung fiir die Durchbiegung w(x, y) einer homogenen 
isotropen Platte lautet in diesem Fall 


Adw(x,y) =< psy) - (1) 


mit dem Laplaceschen Operator A und der Plattenkonstanten 


. Eh? 
=P») 


sowie der Flachenbelastung p(x,y), abhangig von den kartesischen Koordinaten x, y. 


N 


2. Aufpunkt im Plattenbereich. Eine Liésung der homogenen Differentialgleichung (1) ist 
die Greensche Funktion, die die physikalische Bedeutung der Biegeflache einer in x, yy angreifen- 
den Einzelkraft vom Betrage P, hat. Fir die Last P, = 1 wird diese Biegeflache durch Anwen- 
dung des Maxwellschen Vertauschungssatzes die EinfluBflache der Durchbiegung fiir den Punkt 
(%p. Yo), den sogenannten Aufpunkt. Die Greensche Funktion G(%p,¥9; x,y) setzt sich zusammen 
aus einer Summe von Funktionen G, mit gewissen Singularitaten (im folgenden kurz Singularteil 
genannt) und einer Summe G, von reguldren Funk- 
tionen (kurz Regularteil): 


G(xo> ¥os XY) = Go(%o» Vos XY) + Gi(%or Yo % ¥)- — (2) 
Der Singularteil G, lautet? fiir einen Aufpunkt im 
Inneren des Bereiches (Abb. 1) 


r 
0 


Iz 
Gites Vos ts) = SNA rT Invr (3) 


(To 
Der Regularteil G,(%, yo; x, y) besteht aus einer Summe ee) 


von Potential- und Bipotentialfunktionen, die mit F,, 
bezeichnet sind: 


Ci (00 Yo3 %>¥) = 2! Cm Fig (4) 
Die freien Konstanten C,, bestimmen sich aus den 
Randbedingungen fiir Abb. 1. Aufpunkt und laufende Koordinaten in der 


Plattenebene. 


C(%o, Yos % Y) + 

Neben der BiegeeinfluBflache interessieren auch die EinfluBflachen der Schnittmomente und 
Schnittkrafte. Wie bereits von Pucher*® gezeigt, erhalt man die EinfluBflachen fiir die Schnitt- 
momente und Schnittkrafte, indem man die Ableitungen der BiegeeinfluBflache nach den Auf- 
punktkoordinaten xp, yy in die Ausdriicke fiir die SchnittgréBen einfiihrt. Die Randbedingungen 


* Auszug aus der von der Fakultat fur allgemeine Wissenschaften der Rheinisch-Westfalischen Technischen 
Hochschule Aachen genehmigten Dissertation. Referent Professor Dr. F’. Schultz-Grunow, Korreferent Pro- 
fessor Dr. F. Reutter. 

1 F, Schultz-Grunow, Z. angew. Math. Mech. 33 (1953) S. 227. 

2 Ph. Frank und R. v. Mises, Differential- und Integralgleichungen fiir Mechanik und Physik, Bd.I, 
S. 56 ff. New York 1943. 

3 4, Pucher, Die MomenteneinfluBfelder rechteckiger Platten, Berlin 1938. 
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fiir diese EinfluBflachen bleiben die gleichen wie bei der Biegeflache; sie werden spater naher 
beschrieben. 
Nae der Differentiation hat also der Regularteil wiederum die speziellen Randbedingungen 
zu erfiillen. Man erhalt dann fir die EinfluBflachen folgende Ausdricke: 
1) EinfluBflache der Durchbiegung: 


1 ' 
W(Xq> Vos XY) = Rae aN eet Vee a Cn Fin(%o> Yo3 % ¥)]- (5) 


2) EinfluBflachen der Kriimmung K,,, Ky, baw. der Biegemomente M, ,, Mg ,: 


G2) Gra lA Co ii 3 90), 
0?w 1 2 (x — x)? m eee ae 
OP 2. Gas Pgs Sop Vas oan) 
0?(wo) 1 2 (y= Yo)" m 
ia eae Te ‘| act anads aca ay’ Oy 
Me = Koy +» Ky, (ed 


3) EinfluBflache des Verdrillungsmomentes Moxy? 


oe Gril 2¥Yo> * ] 
Moe Ne (= »)] 2 @ 49) (y= yo) 2 (Xo> os Xs) i a 
Oxy OX, AY 8x r OX) OV 
4) EinfluBflachen der Querkrafte: 
0Ao Wo N _ 1 4 (x ae Xo) ' 0A, x Gr en cas Yoo xy) (8a) 
Qox vp mee Axo pm 8x r2 dx, ? 
BA, wy 1 [4(y—yq) , 240 Cm Fin(%o» ¥os ¥> ¥) io 
Dime tanta ee OE r? 5 2Yo : ( 
_ Pwo , Ow 
ace aaa 


Da in den Ausdriicken (6) bis (8) die Ableitungen von >) C,, F., (9,93 x.y) lediglich eine Summe 
m 
von Bipotential- und Potentialfunktionen bedeuten, kann man fiir alle EinfluBflachen anstatt 
dieser Ableitung die Funktionssummen selbst verwenden, wobei nur die Konstanten C,, andere 
Werte annehmen. 
Je nach der Lagerungsart der Platte treten nun folgende Randbedingungen auf: 
1) Eingespannter Rand: 


r) OG(Xo. ¥o3 Xs 
Wy Cla ) gn eet ae (9) 
2) Frei drehbar gelagerter Rand: 
PGC(Xo Vos X> O2G (Xo Vosees 
ty == G(x; 403% Y= Op AM ee Soleie (10) 


3) Freier Rand: 
M,, = 0°G(Xq. Vos XY) ae, eC Naea) | 


on? ar ? 
AGC (xo, Vos XY) CHE C5575 ca) (11) 
On = on? pa 9) 0s 


Darin bedeuten n die auBeren Normalen-, t die Tangentenrichtung am Rande, M), das Stiitz- 
moment und C,, die Stiitzkraft. 


3. Aufpunkt am freien Rand. Besitzt die Platte einen freien Rand, so sind die EinfluBflachen 
der Durchbiegung und der SchnittgréBen fiir einen Aufpunkt an dem freien Rand von Interesse. 
Kin analytischer Ausdruck fiir den Singularteil einer Biegeflache mit dem Aufpunkt auf dem freien 
geradlinigen Rand wird von A. Nadai? in Polarkoordinaten angegeben: 
lity 
i» 


P . -2 
Sl ea eters 3 rinr + (?p sin2.p —r? cos 2yInr) + 7, F (12) 


1 A, Nadai, Elastische Platten, S. 39, S. 206. Berlin 1925." 
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Es gelten hier folgende Beziehungen zu den kartesischen Koordinaten: 
2 a: i ee 
Ge Se > sae we tg. 


Die Einzelkraft liegt im Aufpunkt x, = 0, Yo = 9. Die y-Achse liegt hier in der freien Begrenzungs- 
geraden und die Platte in der Umgebung des Aufpunktes auf der Seite x > 0 (Abb. 2). In Fort- 
fiihrung der Arbeit von A. Pucher!, der die Herleitung der Singularteile fiir einen freien Platten- 
rand ankiindigte, wurde diese dann von W. Dworzak?2 angegeben. Er fiihrte in (12) kartesische 
~ Koordinaten ein und, um die Einflu®flachen der Schnittkrifte zu erhalten, die Aufpunktkoordi- 
nate Yo. 

Dabei wurde die Behauptung aufgestellt, da8 der Aufpunkt an den freien Rand gebunden sei, 
infolgedessen eine Aufpunktkoordinate Xj nicht existiere und auch nicht nach ihr differentiert 
werden kénne. Eine solche Auffassung ist abzulehnen, da selbstverstandlich die Last auch in 
«-Richtung, also in den Innenbereich der Platte geriickt werden kann. Es wire sonst ja auch nicht 
mdglich, die EinfluBflachen fiir die Schnittkrafte zu erhalten, bei denen eine Differentiation nach 

ry %) in jedem Fall erforderlich ist. Weiter wurde von Dwor- 
zak von vornherein in (12) die Querkonstraktionszahl » = 0 
gesetzt, welche Spezialisierung erst nach der Differentation 
gemacht werden diirfte. 


( Y, Y) 
treier Rand, 
der Platte & 
Y 
Abb. 2. Lage des freien Plattenrandes. Abb. 3. Veranschaulichung yon F, 


Es wird daher von der Uberlegung ausgegangen, daB die Potentialfunktion 
F, =r (y sin 2g —cos2qInr), (13) 
welche in (12) auftritt, die Aufgabe hat, die in den Randbedingungen des freien Randes im Auf- 
punkt auftretende Singularitat von r?Inr auszugleichen; (13) mu8 also im Aufpunkt ebenfalls 
Singularitaten besitzen; bei seinem Verschieben ins Platteninnere muf diese Funktion (13) aber 
im ganzen Bereich einschlieBlich des Randes einen reguléren Verlauf haben. 
Stellt man diesen Sachverhalt in Koordinaten dar, so ist statt F’, in (13) zu schreiben (Abb. 3) 


F, =r? (9 sin2 gp —cos2qgInr). (13a) 

In kartesischen Koordinaten lautet (13a) (Abb. 3) 
Fi, = (% + %)? — (y — Yo)" In Ve + 2%)? + (¥ — Yo)? —2 (% + %)(¥ —Yo) arctg a . (14) 
Der Ausdruck (14) ist in dem Punkt der Zahlenebene x = — xy, y = + Yo, der nur auf dem freien 


Rand oder, wenn der Aufpunkt vom freien Rande wegriickt, aber nur auferhalb des Platten- 
bereiches liegen kann, singular und hat dort die physikalische Bedeutung einer Scherkraft®. 
Der Singularteil G)(xp, yo; x,y) der Greenschen Funktion fiir einen Aufpunkt am freien Rand 
langs der Achse y = 0 laBt sich nun also schreiben 


Glo, os 29) = 2 ttn — 14 {Oe + 2g)*— (9 —90)*1 Ye EP FO IOP 
— 2 (x + %) (y —yo) arctg seed sare (15) 


Um den Faktor A zu berechnen, stellt man sich vor, daB an den freien Plattenrand eine symme- 
trische Platte angelegt wird, so daB der freie Rand eine symmetrische Trennlinie bildet, und auf 


1 4. Pucher, Ing.-Arch. 14 (1943/44) S. 258. 
2 W. Dworzak, Oesterr. Ing.-Arch. 1 (1947) S. 67. 
3 Siehe FuBnote 1 von S. 256. 
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dem Erganzungsteil ebenfalls eine Last P = 1 aufgebracht sei (Abb. 4). Dann mu die Querkraft 
in einem kleinen Kreisschnitt um den Aufpunkt zusammen mit der Scherkraft im Aufpunkt den 
Lasten P — 1 das Gleichgewicht halten. (Zu dieser Auffassung ist man berechtigt, da eine spezielle | 
Randbedingung noch nicht eingefithrt wurde.) 


Die Querkraft Q, des obigen Ausdruckes (15) iiber einen kleinen 


Kreis vom Halbmesser r betragt? 
¥ 
Platte Ae i 
Sper Q aia oat, Cork ate Bate (16) 
<—{e : pero ee 
=< Die Scherkraft betragt fiir (15) nach A. Nadai 
Erganzungsteil . RET, C= ») (M, ees M, ; ) 
te ait she a, 26, Las mee 
Abb. 4, Zur Berechnung von A. —A (1 — v) os hected ier —— OF oy ay —A-2 (1 ol. y) Tus ) 


wobei y = artg (y — yo)/(« + x) gesetzt ist. 
Es wird nun aus (16) und (17) 


25,+0,=82A42-2(1+7)7A4=2py 


18 
mit A= Po ( ) 
2(3+7)-2~ 

Mit (18) wird der Singularteil der Greenschen Funktion fiir den Aufpunkt am freien Rand 
= P Te 
© =reE Han [tint p=, {le +? —O — yl FH + O90) 
—— ¥ r2 
— 2 (x+4»)(y — yo) are tg 7— z + rae | ; (19} 


Sollen die Stiitzbedingungen von m Randern einer Platte beliebig sein, so erhalt man durch Super- 
position aus (19) 


Gy = [eine 


Mv 
i 


N = [ln Sei ed ey Vn + Xo)? (ce 9h 


— 2 (tt +9) * (Yq — Yo) are tg eal io | (20) 


m + Xo lt» 


Darin beschreiben u,, = («—a,,)-cosp,,—(y—b,,) * sing,, und v,,=(«—a,,) sing,, + (y—B,,) 
cosy, die Lage eines Randes. Der Ausdruck (20) stellt die allgemeinste Gestalt des Singu- 
larteils der Greenschen Funktion fiir eine Platte endlicher Ausdehnung mit beliebigen homo- 
genen Randbedingungen dar. Liegt namlich der Aufpunkt im Inneren, so kann die Summe der 
Funktionen in (20) weggelassen werden, da sie im Plattenbereich regular sein mu und infolge- 
dessen dem Regularteil zugeschlagen wird, so daf allein der friithere Ausdruck (3) tbrig bleibt 
und A = P,/8a wird. 

Liegt der Aufpunkt auf einem freien Rand, so wird nach (18) 

Te 1% 
— 234+ 0)/a * 
Dabei ist eine der Funktionen in dem Summenausdruck von (20) singular und beseitigt, da sie, 
wie oben bemerkt, die gleichen Singularitaten wie r?Inr besitzt, die Singularitat in den Rand- 
bedingungen am freien Rand. Mit anderen Worten: la8t man den Aufpunkt (xj, yo) vom Innen- 
bereich der Platte zum freien Rand wandern, so muB sich ein Teil des Regularteils der Greenschen 
Funktion zu einer unendlichen Reihe zusammenfassen lassen, die die Gestalt von (14) annimmt 
und die die am freien Rande in Biegemoment und Stiitzkraft auftretenden Singularitaten des 
Singularteils der Greenschen Funktion ausgleicht. 

Riickt der Aufpunkt an einen eingespannten Rand, wird wieder eine der Funktionen in dem 
Summenausdruck von (20) singular und beseitigt die Singularitaten in den Randbedingungen des 
eingespannten Randes, wenn man vor dem Summenausdruck das positive Vorzeichen einfihrt. 
Ks ist dabei fiir den eingespannten Rand selbst y = 0 zu setzen. Es wird 

ee P, 
A =— ee 


1 Siehe FuBnote 2 von S. 255. 
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Man erhalt die Singularteile der Greenschen Funktion fiir die EinfluBflache des Stiitzmomentes 
M,, und der Stiitzkraft C,, durch entsprechende Differentiation von (20). Fiir einen eingespannten 
‘Rand Jangs der y-Achse fiir x = 0 gilt fiir den Grenziibergang x,—> 0 


- 3, ®PGC (xo, ¥o3 XY) Py x 
Ma. x2 Sm In r? + =r) (21a) 
C= BE Yo *Y) __ y PCr 03 ¥Y) 
bie Ox? Oxy Oy? 


(21b) 


o Beh ee 

--2-(5+2) 
Die Funktionen In 7? und x/r? in den Ausdriicken (21a) und (21b) kénnen weggelassen werden, 
denn sie sind Potentialfunktionen, und es gilt der Satz, da® diese im Inneren des Regularitats- 
bereichs keinen Extremwert haben kénnen. Es wird also ihre Summe mit den Uberlagerungs- 


funktionen bei Erfiillung der ersten Randbedingung von (9) im ganzen Innenbereich Null. Man 
erhalt 


P i 

Me (21c) 
Dd eos 

Dy eee (21d) 


als gewissermaBen modifizierte Singularteile fiir die EinfluBflachen von Stiitzmoment und Stiitz- 

kraft, wie sie schon von J. H. Michell! und A. Pucher? auf anderem Wege gefunden wurden. 
Fiir den freien Rand miissen die EinfluBflachen fiir das Stiittzmoment M,, und fiir die Stiitz- 

kraft Cy, wegen der Randbedingung (11) Null sein und ebenfalls ihre Ableitungen nach yy: 


_ PC(xo os XY) PC (Xo Yos XY) _ 
LOG ax? ie dye ee me 
__ PGC (Xo ¥o3 XY) PC(X Vos ¥Y) __ 
OF ae a (2 — y) rr = () D) (22b) 
@Mox _ PC(xo yor XY) PC(XoYos XY) _ 1) 
ayy Ox§ Oo 4 dys a ee 
Aus (22a) ergibt sich 
Caan VN aay tees ea es #9) (22d) 
Ox? ays f 
aus (22c) 
BCX VoXY) __ P PC(Xo Yo XY). (22e) 
dys 0x5 Oo : 


Die Bedingungen (22d) und (22e) hat man bei der Bildung der Einflu®flachen fiir das Biegemoment 
M), und die Querkrafte Qo,, Qo, einzufiihren, wahrend sie in der EinfluBflache fir das Verdrillungs- 
moment M,,, keine Bedeutung haben, da in diesem Ausdruck nur die erste Ableitung nach x, 


und y, vorkommt. Man erhalt so die Singularteile der Einflu8flachen der Schnittkrafte Moy» 
Days oO Moxy fiir einen freien Rand langs der y-Achse fiir x) > 0, wie folgt: 


M,, 2 ae Ly an We (1 “Me 7?) OG, mi (1 a v) iE In r2 (1 = v) 2 ; (23) 


ve 8 2 G+naAN 
Qoy = oa | oa a NGF ”) fee PM pea ee) 
Bet tase) pa) awa ner tO a, em) 
Moxy = — (1—») oe = wis a : ee) | (1 +») aretg eo 3 (25) 


Die EinfluBflachen Q,, und M,, y sind nur wegen der Kirchhoffschen Annahmen vorhanden. 


1 J. H. Michell, Quart. J. Math., 32 (1901). 
2 A. Pucher, Ing.-Arch. 12 (1942) S. 76. 
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4. Die Bedeutung der konformen Abbildung und die Auswahl der Funktionen. Nach dem | 
Verfahren von F. Schultz-Grunow ist der Innenbereich einer polygonal umrandeten Platte konform 
auf das Innere des Einheitskreises abzubilden. Man denkt sich dazu die Platte mit einem Gitter- | 
netz iiberzogen (Abb. 8). Der Plattenrand geht in den Kreisrand tiber, die Abbildung der Um- 
gebung der Plattenecken erzeugt Haufungspunkte der Bildpunkte aquidistanter Gitternetzpunkte. 


Das bedeutet, da® eine im polygonalen Bereich erklarte Funktion wie z. B. G(%9, Yo; x, Y) durch | 
die konforme Abbildung an den Bildpunkten der Ecken singular wird. Die Ordnung dieser Sin- | 
gularitaten ist festgelegt durch den AuSenwinkel f der polygonalen Umrandung und betragt | 
(n B —2)/, wobei n jede ganze Zahl sein kann. Der Regularteil der Greenschen Funktion wird 
aus im Platteninneren regularen Potentialfunktionen aufgebaut, die im Kreis definiert sind. Diese 
Funktionen bleiben bei der konformen Abbildung Potentialfunktionen. Sie seien mit @,, und 
X,, bezeichnet, wobei ®, stets den zur x-Achse symmetrischen, X,, den zur gleichen Achse anti- 
metrischen Teil bezeichnen soll. Aus diesen Funktionen werden die Bipotentialfunktionen y,,, | 
V,, U,, gebildet, in dem die Potentialfunktionen ®, und Xx, mit x, y oder r multipliziert werden: | 

Y =xoD) Voy xe (Phe Ode (26) 
Es werden die Randbedingungen der Platte auf dem Kreisrand erfiillt, was sich als zweckmabig 
fiir die numerische Rechnung erwiesen hat. 

Diese Rechnung gliedert sich dann in zwei Abschnitte. Zuerst miissen die oben beschriebenen | 
Ecken-Singularitaten eliminiert werden, was mit einer speziellen Art von im Inneren des Kreis- | 
bereiches regularen Funktionen, die hier als ,,Eckenfunktionen“ bezeichnet sind, geschieht. Man 
entwickelt dazu diese Funktionen unter Beriicksichtigung von je zwei Randbedingungen von der | 
Ecke aus langs der zwei Rander in Potenzreihen und iiberlagert sie dem ebenfalls in eine Potenz- | 
reihe entwickelten Singularteil der Greenschen Funktion. Dabei werden die Glieder, die die Sin- | 
gularitat verursachen, identisch zum Verschwinden gebracht. Wiirde man alle Glieder der er-— 
wahnten Potenzreihe mittels der Eckenfunktionen zum Verschwinden bringen, so bedeutete dies 
die identische Darstellung des Singularteils der Greenschen Funktion und somit die Herstellung der 
trivialen Lésung. Deshalb sollen nur die Randbedingungen in den Ecken exakt erfiillt werden, — 
indem von den Potenzreihen der Eckenfunktionen und der Greenschen Funktion keine Glieder mit 
Exponenten > 1 beriicksichtigt werden. i 

Im zweiten Abschnitt der Rechnung werden dem so ausgeglichenen Singularteil Funktionen 
iiberlagert, die aus der analytischen Funktion €(,,6) in der Kreisebene hergeleitet werden. Mit 
diesen, hier ,,Mittelpunktsfunktionen“ genannten, Funktionen wird der Verlauf der mit den 
Eckenfunktionen tberlagerten Greenschen Funktion G(x, yo; x,y) den Randbedingungen an- 
gepaBt. Da die Eckensingularitaten bereits ausgeglichen sind, ist zu fordern, daB die mit C(j,o) 
konstruierten Bipotentialfunktionen in den Eckpunkten keine Singularitaten besitzen. Weil, | 
wie spiter gezeigt wird, die aus C(%,o) hergeleiteten Funktionen unter Umstanden mit den Ecken- _ 
funktionen identisch sein kénnen, mu8 zum anderen gefordert werden, daB sie nicht alle Potenzen 
des in Reihe entwickelten Singularteils bzw. nicht die gleichen linearen Abhangigkeiten enthalten 
diirfen, um nicht die triviale Lésung zu erzeugen. 

Zur Begriindung der Wahl gerade dieser beiden Funktionsgruppen in der Kreisebene ist anzu- 
fiihren: Die Eckenfunktionen erfiillen die Randbedingungen in der Ecke exakt. Da die Ecke in © 
der Kreisebene Haufungsstelle der Bilder aquidistanter Randpunkte der Gitternetzeinteilung ist, | 
werden von vornherein die Randbedingungen in guter Nahrung itiber einen relativ groBen Teil — 
des Polygonrandes erfiillt. Die Mittelpunktsfunktionen aus dem Kreis haben den Vorteil einer 
einfachen Bauart (langs des Kreisrandes ist r = konst.) und zum anderen eine bedeutend geringere 
Oszillation wie die entsprechenden Funktionen in der Polygonebene?. Hieraus ergibt sich die 
gute Konvergenz des Verfahrens. Weiter kommt noch hinzu, da8 die Funktionen @,, X,, V,, 
U,, voneinander linear unabhangig sind, wahrend bei einer Bildung von entsprechenden Potenz- — 
funktionen in der Polygonebene an Stelle von ®, und X,, zwei lineare Abhangigkeiten zwischen 
jeweils vier zu bildenden Funktionen auftreten wiirden, wie sich leicht aus dem Aufbau dieser | 
Potenzfunktionen zeigen 14Bt. Somit steht die doppelte Anzahl von Funktionen zur Verfiigung, 
was wegen der geringeren Oszillation eine weitere Konvergenzverbesserung des Verfahrens bedeutet. 
AuBerdem ist es bei diesem Verfahren auch méglich, auf einem Rand die Randbedingungen ver- | 
schiedener Lagerungsarten zu erfiillen. 


5. Berechnung von Einflufflichen einer Dreieckplatte. Im folgenden sollen die Biege- und 
die BiegemomenteneinfluBflache einer Dreiecksplatte berechnet werden. Das Dreieck sei gleich- 


1 Siehe FuBnote 1 von S. 255. 
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schenklig und rechtwinklig. Die beiden Katheten seien eingespannt, die Hypotenuse frei. Der Auf- 


punkt befinde sich in der Mitte des freien Randes. Die Kathetenlange sei a = V2, die Querkontrak- 
tionszahl wird mit y = 0,3 angenommen. Wegen der Symmetrie wird die Betrachtung nur fiir die 
eine Halfte des Dreiecks ausgefiihrt. Das Dreieck liege in der komplexen Ebene 2(x, ty) baw. in 
Polarkoordinaten ausgedriickt z(r, ¢), der Kreis in einer Ebene C(g, i) baw. C(o, @). Es ist dann 


Way ree Ce in Soe. 


Fiir die Abbildung des Innenbereiches eines Kreises der ¢-Ebene auf das Innere eines Vielecks 
in der 2-Ebene besteht die Schwarz-Christoffelsche Abbildungsformel: 


— Ae 


a 
ze fal—tk) = a +B. (27) 


Hier seien ¢x die Bildpunkte der Ecken in der Kreisbene: Ay ist die Differenz zwischen dem Aufen- 
winkel einer Ecke und z. B, C und ¢, stellen komplexe Konstante dar, die die Lage und GréBe 
des Vieleckes in der z-Ebene bestimmen (Abb. 5). Fiir ein gleichschenkliges, rechtwinkliges Dreieck 
ist Ay, = 2/2; Av, = Avs = 32/4. Die Abbildungsformel (25) geht damit in den Ausdruck 


c 


us LB (28) 


2=C 
| C=) C6) 4G) 


Abb. 5. Konforme Abbildung des Dreieckes auf den Einheitskreis. Abb. 6. Zur Ermittlung der Funktion €—¢,. 


iiber, wenn fortan z fiir z/a geschrieben wird. Da auf den Einheitskreis abgebildet wird, gilt fiir 
seinen Rand g = | und auberdem 


Cae, do =e dO. : (29) 


Fiir die Bildpunkte der Ecken soll 
2 4 


Cece ete A, eave: 9 eee (29a) 


gelten. : iat 
Es wird nun der Verlauf der Funktionen (—<,), (C—C,), (¢—¢;) auf dem Rand des Einheits- 


Ad 2 4 4 2 ; 
kreises in den Bereichen 0<O< 3 zt< O< 3% und 3% <O< z% also zwischen den 
Bildpunkten der Ecken untersucht. 
74, c 
Im Bereich 0<O0< 3% gilt 


sin O 3) 
lp —€; Sy rae on = | y; sin > 
ia (F—) 
Demnach ist (Abb. 6) 
C¢—= | 2 sin — 2 +3) 
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In Tabelle 1 sind die Randverlaufe der Funktionen zusammengestellt. 


Tabelle 1. Randverldufe der Funktionen (€ —¢,), (€ —l.), (€ — 93) 


| scocts | Sacocte Cer 
c—h | 2 sin -~ iG +9) | Bang iG +7) 2 sin iG+7) 
Se (cae aoc ce 
rom [faa (S + Se aan e095) nan (z 4 8) ted 


Wird eine Funktion bei der Annaiherung an eine Ecke = 0, so tritt beim Uberschreiten der Ecke, 


wie Tabelle 1 zeigt, ein Argumentsprung von —z auf. Der Argumentsprung um —z ergibt sich | 
aus der Bedingung, daB die Werte, die (¢—¢,) annehmen kann, auf einem bestimmten Blatt | 


der Riemannschen Flache-der Funktion (¢—C,)” liegen miissen. Dies sei durch Abb. 7 veran- 
schaulicht. S 

Driickt man die Randverlaufe der Funktionen (€—¢,), ((—,), (C—¢3) nach Tabelle 1 aus. 
so erhalt man in der Abbildungsformel (28) fiir die Randverlaufe t,, t,, t; (Abb. 5) folgende Aus- 
driicke: 


(2) 
: dO 
ere 7 1/2 iG O\ \3/4 rs @\ 3/4 | b, (30a) 
2) sin — ee Sith ee ee 
2 3 2 (3 4] 
0 
2/3 2 
; dO 
a iG 1/2 U Q\ |3/4 F QO\ \3/4 b, (30b) 
Tai ae cae 
(0) 
1} 
53 dO 
aos |. O12) (x O\/34|.. /x . 60\/3/4 +b, (30c) 
[meh =O) 
1/30 


die sich nach der Substitution von 


Abb. 7. Randverlauf der Funktion (€ —{,) [ein 

um —¢€, verschobener Kreis] sowie Argument- 

sprung um —z, wenn man den Bildpunkt der 

Ecke iiberschreitet, wobei die Pfeile den Vektor 
¢ —¢, fiir €¢—>€¢, darstellen. 


sein muB. Daraus ergibt sich Be 


("2 


cost g~ = 


(o) 
cos > 


1—(1 


S sin? ae 
3 2 


F durch die elliptischen Integrale erster Gattung mit dem 


Modul k = 1/y2 wie folgt darstellen: 


t; = c, F(p; k) + B, 
tp = €, F(p;k) + B, 


ts = €, F(p; k) + B, 


? 1% 

——, 1 
Pr 2 (3 a) 
|p ys 

=, b 
a 2 (SB) 
? 
_~ = Oe 
7 (31c) 


Fiir 1, gilt wegen der Symmetrie zur x-Achse (Abb. 5), daB fiir 


= hey ty een te 
9 >] Cire of 


C, 


1 
4 
rie 


= 0,53935260 ; 
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fiir t, und ¢, besteht die Bedingung, da8 ihr AbbildungsmaBstab auf der eingespannten Seite ein 
Minimum haben muB, und da 


P 7 ee 
2 Pr a . a y2 2 


also gleich der Seitenlange ist. Das Minimum tritt auf fiir cosy = — 1, p= +2 nach (30a) 
und (30b). Man erhalt dann fiir 


O= 0-1, =, O =i n> t,=0 


i cares 
9 => tat be |e = 125 
ferner . 


== = 0,38173988 , 


Die Randpunkte bilden sich dann, wie in Tabelle 2 aufgefiihrt, ab. 


Tabelle 2. Abbildung der Randpunkte 


ol] . 61] 2 


® 
= 
S 
ion 
> 
a 


A 0,0 67,5 0,48517668 120,0 0,0 

; 0,03828416 75,0 0,52426809 120,5 0,38742822 
5 0,05414223 82,5 0,56674445 121,0 0,34185739 
3 0,07657291 90,0 0,61456046 122.0 0,40665179 
aU: 0,10486917 aS 0,67103242 123,75 0,47610211 


0,14541084 105,0 0,74293306 127,5 0,56691652 
0,21047452 112,5 0,84943917 135,0 0,67641158 
0,25900145 116,25 | 0,93901050 142,5 0,75179931 
0,30109830 118,0 1,00797160 150,0 0,81221233 
0,33957706 119,0 1,07252839 157,5 0,86444182 
0,37630150 119,5 1,12695751 165,0 0,91187194 
0,41215515 120,0 1,41421356 172,5 0,95656355 
0,44814603 176,75 | 0,97835773 
178,0 0,98846699 
179,0 0,99423506 
179,5 0,99711761 
180,0 1,0 


’ 


SHUN SHMYUMAWNH SOS 
ounouounonnAnooune 


DNukwwnre 


ey 


Zur Abbildung innerer Punkte des Dreiecksbereichs auf das Kreisinnere wird nochmals an (28) 
angeknipft: 


fi St set a (Ge Py ta (St) fe a a (32) 


mit der speziellen Lage der Bildpunkte der Ecken nach (29a). Durch Substitution 


1/—4€—4 —5. 0) C—E) C—G) _ /_  — 0? 
= —— 5 = 33 
corp = J CBG — by G@+i+ 4 Sr 
wird das Integral wieder auf die Form eines elliptischen Integrals erster Gattung gebracht mit 
b= 1/2: 
ES Ke e 
z—=c F(g;k) be JCB (34) 


Die Konstanten ¢, b, y, errechnen sich aus den folgenden Anfangsbedingungen: 


1 & a5; 
I)o=0, ND as sae BE, z ) 


NE=1, cosp—0: ga=>, s=1=6 He;h) 


3)€=—1, cosop=1: yp=0, z=0=c F(g;k) 
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Daraus ergibt sich i 
pe z2—b ) 
#5 _ Fp;k, gp =am2 
C c 
ot elie one eg) US), 
cos p = cn =en( = 4iX}— @—t+) aa 
—b a—b —b 
14 2ent(7=") | on ( =) Y/12(1— ent? 
= c c c 
Aes ey Gall) 
G 
Es bedeutet hierin Ps 4 
a oe eee Oe dh ee 
z2—b c c ; c c c c 
cn —— = = —1 ; : 
1 sn? (2?) a aA 1 — sn? ES dn? 2) 
C c c C 


Die Trennung von Real- und Imaginarteil ergibt 


&=o0cos@O=Re[f], @=f/&2+7?, 


Horn? — incl te Or 


att 
E° 


In Tabelle 3 ist die Zuordnung der Punkte der Dreieck- und der Kreisebene angegeben (Abb. 8). 


Tabelle 3. Abbildung der Gitternetzpunkte des Dreiecks auf den Kreis 


00 | —1,0 0,0 20 | —0,85152883 | +0,52430785 | 50 | —0,56722327 | +-0,82356401 
01 | —0,73306117 0,0 21 | —0,66328894 | +0,43812174 | 51 | —0,51920015 | +0,78467926 
02 | —0,50949597 0,0 22 | —0,48377889 | +0,40000582 | 52 | —0,41245361 | +0,77150196 
03 | —0,30249597 0,0 23 | —0,30519426 | +0,39328676 | 53 | —0,37738982 | +0,87979429 
04 | —0,09632627 0,0 24 | —0,11814602 | +0,40926982 | 54 | —0,29335392 | +0,86253259 
05 | +0,11841868 0,0 25 | +0,08762038 | +0,44108809 | 55 | —0,19999997 | +0,97979593 
06 | +0,34461158 0,0 26 | +0,32290223 | +0,47766695 | 60 | —0,52803211 | +0,84922443 
07 | +0,57496061 0,0 27 | +0,59516350 | +0,49522055 | 61 | —0,50214191 | +0,82972098 
08 | --0,78690954 0,0 28 | +0,89378831 | +0,44848909 | 62 | —0,46422604 | +-0,83111801 
09 | +0,94217548 0,0 30 | +0,73287245 | +0,68036605 | 63 | —0,42193193 | +0,85979713 
010) +1,0 0,0 31 | —0,60387659 | +0,59709917 | 64 | —0,38189799 | +-0,92420447 
10 | —0,95666366 | +0,29119520 | 32 | —0,46445690 | +-0,56099438 | 70 | —0,50894114 | +0,86080133 
11 | —0,71324625 | +0,23384466 | 33 | —0,31611417 | +-0,56163934 | 71 | —0,49696756 | +-0,85368073 
12 | —0,50201517 | +0,20894721 | 34 | —0,15395387 | +0,59394916 | 72 | —0,47933107 | +.0,86020679 
13 | —0,30254831 | +0,20253746 | 35 | +0,03356481 | +0,65460699 | 73 | —0,46360522 | +0,88604186 
14 | —0,10111004 | +0,20823544 | 36 | +0,26523427 | +-0,73758210 | 80 | —0,50177129 | +0,86500031 
15 | +0,11188020 | +0,22077770 | 37 | +0,57306536 | +0,81950968 | 81 | —0,49769423 | +0,86381615 
16 | +0,34121249 | +0,23283384 | 40 | —0,63449126 | +-0,77293003 | 82 | —0,49288978 | +0,87009176 
17 | +0,58263607 | +0,23165010 | 41 | —0,55304070 | +-0,71016515 | 90 | —0,50011079 | +0,86596143 
18 | +0,81480089 | +-0,19863935 | 42 | —0,45309090 | +-0,68612330 | 91 | —0,49955678 | +.0,86628113 
19 | +0,99293983 | +0,11861931 | 43 | —0,34017569 | +0,69883658 |100 | —0,5 -++0,86602540 

44 | —0,21298158 | +0,74951277 

45 | —0,06263569 | +0,84363658 

46 | +0,13360166 | +0,99103513 


Um nun die Eckensingularitaten in der Kreisebene ausgleichen zu kénnen wird CoG@osVos 4) 
iiber den Bildpunkt von der Ecke aus auf dem Kreisrand in eine Reihe nach steigenden Potenzen 
von @ entwickelt. Zunachst wird die Entwicklung der Abbildungsfunktionen, die durch (31a) 


bis (3lc) gegeben sind, durchgefiihrt. Man erhalt folgende Anfangsglieder der Reihe bis zu Poten- 
zen S 7/4: 


dt, 1 13 

tte Ea 

dt, 1 2; — 3/4 1 ‘ 2 1/4 
rarest 0) eayea aaa 
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Die Integration mit den Werten aus (31 ff) ergibt 


2 13 2 
= es c (ove mag est mit oa c, = 0,40981935 , (36a) 
/3 /3 
22 2 ui tee(2 aM oy 
Ro dey (ee = 6) oe (Fx = 6) | mit 27? , —0,94029173,  (36b) 


— 3 3\3 va 
3 40 3 3 

= eee 2 1/4 1 2 5/4 ~ 

ty =3— (@ is | ae (0 7 os x) | mit 28 Co= 0,94029173 (36c) 
ya. ; 40/3 /3 


und damit an der Ecke | (Abb. 5) 


t t 
x = 1 ——+ =1 — 0,289786 OV?, y =— = + 0,289786 @!? 
V2 Ve (37a) 
r? = 1,0 — 0,579572 O'? + 0,167952 O!, 
an der Ecke 2 (eingespannter Rand) 
2 / i/4 
aoe 2 — 0,664887 (32 —6) a ae Re ee Fe 0,668887( 5 rt By 
2 1/4 fe 1/2 >) 
r? = 1,0 —1,329773 (3 cc 6) + 0,884149 (37 — 6) : 
an der Ecke 2 (freier Rand) 
me yee 7 = 1 0.940099 (0 he (37c) 


fp ly My Ie |e | 
2 \23_ (|b 


4 


~iy “in 
Abb. 8. Zuordnung der Gitterpunkte Dreieck—Hinheitskreis. 


Mit Einfithrung dieser Ausdriicke in (19) erhalt ‘man die Entwicklung von GoltGs Vos 4) an 
den Ecken in Potenzen von @. Ebenso lassen sich auch die Normalableitungen und hdheren 
Ableitungen mit den Formeln (37) an der Ecke entwickeln. Diese Ableitungen kénnen sowohl in 
der Polygonebene wie auch in der Kreisebene durchgefiihrt werden, da sich wegen der konformen 
Abbildung die Richtungen in beiden Bildern entsprechen. Im weiteren sollen fiir die Randbedin- 
gungen alle Ableitungen in der Polygonebene gebildet werden. Es ist dann auf der eingespannten 
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Seite, wenn n die Normalenrichtung in der Polygonebene bedeutet, 


ax ob ay 1. FF 
Lae =—, 5 =J2, (38a) 


én V2" an fa” 


auf der freien Seite 


scp Ia Sa pen 2 Ba (38h) 


én "On Gas 
Fiir eine im Kreis definierte Funktion TfL) lautet die Normalableitung in der Polygon- 


ebene, wie folgt: 


on 


Im (ee Goan (38c) 


R 
{Tmt SO ipl af(é) a% a 
Der Ausdruck d0/dn ist der Ma®stabsfaktor der konformen Abbildung. Man kann wegen (28) 


a _ || e|¢—z,h |¢—e, |¢ — Ca) (38a) 


00 
Oz 


bilden. 


Die Eckenfunktionen, die die Elimination der Potenzen von @ des Singularteils der Greenschen | 


Funktion bewirken sollen, und deshalb dieselben Potenzen in © wie jener enthalten miissen, haben 
die Gestalt 

R R R 

Rehan, fret —eome, [Rebeca (39) 


wie schon aus der Abbildungsfunktion zu vermuten war. Es sind Potentialfunktionen und durch 
Multiplikation mit x, y oder r? erhalt man die zugehérigen Bipotentialfunktionen. 

Bei den Funktionen (39) lassen sich aber bestimmte Beziehungen zwischen ihren Entwicklungen 
in © langs der beiden von einer Ecke ausgehenden Rander herleiten, wie dies an einem Beispiel 


gezeigt werden soll. Es sei die Funktion re (€ —€,)"? und ihre Normalableitung fir die 
Koordinatenbereiche t,, t, (Abb. 5) betrachtet: 


ie Fel (¢ —C,)¥2 = sin oat ia 


Im 


sin 
9 [Re _-\e 1. ft A (cos 5 3 O\ | ag 
an ITeat (2 c= 9 sin( 5 Z| el ae a) dz |? 
;. |Re| (= + vee fe OM leash ax 2) 
ia Tal ¢ @)) =e z) et eta) 


oO Re ro 12a LOW, 4 @\—1/2 cos 4 3 
in Im © ba) sin(3 z at naa el 
Fir O > 2/3 lauten die ersten Entwicklungsglieder 


w= Fae—e= ao) 


mn 

= 

B 
ee, 


(is | ; : (39a) 


SIBJe—or- Ho) 
fi {ref Etat? = (Fa — 0)" [enol (39) 
zfisje—cm—1(85—e)"Eals). om 


Man sieht, daB die Entwicklungsglieder der Potenzreihe in (39a) und (39d), sowie in (39b) und 
(39c) die gleichen Koeffizienten besitzen. 

Ahnliche Zusammenhange lassen sich fiir alle Glieder mit gleichen Exponenten m finden. 
Da sich aber alle verwendeten Potentialfunktionen durch eine Reihe der Form (¢C—C,)™* dar- 
stellen lassen, bedeutet das, daB bei Ausgleich einer Potenz von @ auf der einen Seite der Ecke, 
auf der anderen Seite ebenfalls eine bestimmte Potenz von O ausgeglichen wird. Die auf diese 
Weise gefundenen Bedingungen sind in Tabelle 4 zusammengestellt. 
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Tabelle 4. Zusammenhang der Potentialfunktionen 


Potenzen von Rand 1 Rand 2 
fir m = Funktion Normalabl. Funktion Normalabl. 
oO TFT —__ 
@, As.8n 1216 gaan ae eae ae 
fe -— > x 
276. 10, 14) 13 ar eT 
OP a ae eee, 
[x — x] <> x : 
1:35, 9813-17 te 8 [eee ee 
a7. Al, 16. 19 B Pel 2 
2's > ? [x il ee = - x 


Wegen dieser Beziehungen wurde zur Erleichterung der numerischen Rechnung fiir den Ecken- 
ausgleich der Singularteil (19) dargestellt durch 


G=*6,+G6n, 


wobei G; und Gy Potentialfunktionen bedeuten. Dadurch erreicht man, da Gy und Gy je fir 
sich mittels der Eckenfunktionen (39) ausgeglichen werden kénnen. Das geschieht in der Weise, 
daB die nach Potenzen entwickelten Glieder von G;, Gy und der Eckenfunktionen (39) in einem 
Gleichungssystem gleichgesetzt wurden. Die Rechnung in dieser Form hat den Vorteil, daB das 
hier vorliegende Gleichungssystem von 26 Unbekannten in mehrere kleinere Systeme mit héchstens 
5 Unbekannte zerfallt, da die Matrix desselben auf Grund obiger Beziehungen nur zur Halfte 
besetzt ist. Der in deren Ecken ausgeglichene Singulirteil besitzt die Form 


Cy =e > (Gi ¢,. ©,,) + > Gute, ©,,)> 


worin die Konstanten c,, und C,, sich aus dem Gleichungssystem ergeben. Fiir die BiegeeinfluB- 
flache gilt 


OP — xlnr—yaretg® + 775s, (40a) 
Cup G2 2) Inr— 2.2 y arctg Lee ae) 20) 
Bates ety y Dae Sama (1 29) y)> 
fiir die BiegemomenteneinfluBflache 
es Se (41a) 
2 V=Inr. (41b) 
Tabelle 5. Die verwendeten Eckenfunktionen 
@,, = Im (€ —¢,)4/? ®,,, = Im (¢ —¢,)3/? 
2, = rat [Rede —7,)t/4-P Im (¢ —¢5)¥4] P31 = ea [eG Crete (C5) | 
y2 V2 
P2. = a [Re (¢ —¢,)/2 — Re (¢ —¢,)}/?] P3.. = = ftin'( =<) tm (605) 8/* | 
1 é 1 ae 
P23 = 93/4 [he(C—-¢,)7/4 aa (C=, 3/4] P33 = 9a/4 flr? =£,)9/4—— Het? —-¢3)7/2| 
Po = —— [Re (6 —f,)5/4 + Im 6 —G,)5/4] Pay = —— [Im C —G,)8/4 + Re C—G)9/4] 
2/2 2/2 ; 
1 1 é ; 
Yi = 5 j2 [Re (C— ¢,)3/2 — Re ¢ — ¢,)4/2] D4 = 2/2 [Im (¢ —¢,)3/2 + Im (¢ C3)3/2] , 
1 1 
P27 = 2.95 [Re (¢ —¢,)7/4 — Im (¢ — ¢3)7/4] P= 5 ott (hash) ee (C——7,) 15] 
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Entsprechend der Symmetriebedingung wurden die in Tabelle 5 zusammengestellten Ecken- 
funktionen aufgestellt, wobei die erste Indexbezeichnung die Ecke, die zweite den Wert von m 
angibt. 


Fiir den Ausgleich der tibrigbleibenden ganzzahligen Exponenten von O wurden sogenannte 


Mittelpunktsfunktionen, das sind Funktionen der Form @,, =| 


Re 


o™ e'™® verwandt. Es sel 
Im 


bemerkt, daB die Eckenfunktionen ie (¢ —C,)™ fiir ganzzahliges m sich immer als Summen 


solcher Mittelpunktsfunktionen darstellen lassen, z. B. 


Re (¢ —C,)? = Rel? — 2 Re(C, + Rel 


Rel? — 2 Ref 41. 


Entsprechend den Symmetriebedingungen wurden folgende Mittelpunktsfunktionen aufgestellt : 


Tabelle 6. Die Konstanten des Singularitdtenausgleiches 


G, = Re(C!) = 1, 
Oy =Re(F) =p cos Go 
D, = Re) =e" cos 2 OZ 


cr Cm , ¢ @, = Re(C*) = 6 cos 3:07 
Biegung Moment Biegung Moment D, = Re(C*) es o! Cond 0 ‘ 

1,1 | + 0,725065 | —0,409819 | + 0,900702 | + 0,409819 Mit diesen Ecken- und Mittel- 
1,3 | —0,169795 | + 0,171295 | —0,225176 | —0,125398  punktsfunktionen wurde der 
2,1 | + 0,500266 | + 0,887129 | —5,477503 | +0,6807]8 FV enausgleich durchgefiihrt und 
3,1 | —2,597849 | —0,680719 | + 1,105323 | + 0,887129 Be Spat iee 
2,2 | —0,603885 | -- 1,207769 | -+ 1,892283 | —0,161811 dabei die Konstanten der Tabelle 
3,2 | —0,161810 | -+ 0323620 | + 3,778807 | + 0,603885  errechnet. 

2,3 | —0,089176 | + 0,535058 | + 0,615115 | —0,430581 Zu dem Singularitatenaus- 
3,3 | —0,215291 | + 1,291744 | —0,254790 | +-0,178352 eon 3 

25 | —0,463256 | — 1.801775 | —0,281622 | +0,227307 _ Sleich ist noch 2u een 2 Ae. 
3,5 | —0,386728 | —0,044754 | + 1,458046 | —0,415149 an den eingespannten Randern 
2.6 — 0,149436 | —0,991282 | + 2,090895 | + 0,476945 die Eckensingularitéten ausge- 
3,6 == Rene = Ane = een a eneGle glichen sind, wenn die Glieder der 
Mell == ()}, 82704 T +> armas) I ? + : 1 
37 | +0,048798 | — 27144711 | — 0.506483 | + 0.392581  F otenzreihe 3 Wee n Ce O, 
0 + 2,389813 | + 6,489717 | —0,399717 | —2,323259  dagegen am freien Rand erst, 
1 + 0,609183 | + 2.016635 | + 1,161052 | + 1,196505 wenn die Glieder bis 97/4 eliminiert 
2 — 0,053214 | —1,903354 | —0,038399 | + 0,421469 werden, da in den Randbedina 
3 + 0,013187 | + 0,882574 | —0,008250 | — 0,154710 : 
4 —:010187 | — 0.358393 | —0,001198 | +0.058880  Sumgen des freien Randes die 


zweiten und dritten Ableitungen 
vorkommen. 


6. Punktweise Erfiillung der Randbedingungen. Nachdem der Singularteil in den Ecken 
durch die Uberlagerung von Eckenfunktionen identisch dargestellt wurde, und somit in den Ecken 
die Randbedingungen exakt erfiillt sind, werden dem so ausgeglichenen Singularteil nun Mittel- 
punktsfunktionen in der Weise tiberlagert, daB die Randabweichungen méglichst klein werden. 
Die Mittelpunktsfunktionen als Potentialfunktionen sind mit ®, und X,, bezeichnet, wobei der 
Index m die Potenz von @ bzw. den Faktor von @ bedeutet: 


Da Re Gs realm (ahr 
Die zu tberlagernden Mittelpunktsfunktionen werden so gebildet, da® sie in den Bildpunkten 
der Ecken ebenfalls keine Glieder mit Potenzen von 9” mit m < 1 enthalten. Das wird erreicht, 
in dem mehrere Mittelpunktsfunktionen durch Summation zu einer einzigen zusammengefaBt 


werden. Dabei lassen sich alle méglichen linear von einander unabhangigen Kombinationsméglich- 
keiten ausschépfen. Man erhalt dann den Symmetriebedingungen entsprechend 


Ps =P, +29,—2,56,4+0,5D,, , 
Po =P, + DP —2@,, 

~, =O,+9,—29,, 

Ps =P, + D, —2G,, 

Py =P, +20,—3,, 

Po = Py +29,—3,, 

Pu = D,, +26,—3G,, 

Piz = Dyp + 3D) —49,. 


(42) 
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Durch Einsetzen von X statt ® erhalt man die antimetrische Funktion. Die zugehérigen Bi- 
potentialfunktionen werden wieder durch Multiplikationen mit x, y oder r? erhalten: 


Un 20), C= 2. Vee WeX (43) 


Diese Funktionen ¢,,, y,,, U,, V,, sind nun so dem Singularteil zu iiberlagern, daB die Rand- 
verlaufe folgenden Randbedingungen angenihert werden: 


eingespannter Rand 
Cox ci (Ci, Pm “lp Cm Ym 7 K, Us ie ie Van) = 90 ? 


aG d 
a Ss én 2 (C,, Pm = Cn Wm =e Ke U,, + Ss Va) a 0, 


freier Rand 


oe? 0? 
& Ret i) [Goy ie pa, (Ca Pm a Cm Ym mie biG us a Ks,)| ae 0, 


og 0. 
(Ga + C9) 5 5)[801 +S (Cm m+ En Yn + Ky Uy + Ky Ve] =O. 


Die erstrebte Annaherung an die Randbedingungen wird erreicht durch ihre exakte Erfillung 
in einer beschrankten Anzahl von Punkten. Um aber den héheren Rechenaufwand bei der Bildung 
der zweiten und dritten Ableitungen, die die Randbedingungen des freien Randes enthalten und 
die Randbedingungen fiir Stiitzmoment und Stiitzkraft formulieren, zu umgehen, werden nicht 
diese, sondern das Integral bzw. Doppelintegral iiber dieselben langs des freien Randes an ver- 
schiedenen Punkten zu Null gemacht; denn da zwischen zwei Nullpunkten eine regulare Funktion 
mindestens ein Extremum besitzt, so muB auch ihre erste bzw. zweite Ableitung zwischen bzw. 
in der Nahe zweier solcher Punkte mindestens einmal Null werden. Es werden also durch diese 
punktweise Erfiillung der Integralverlaufe auch Stiitzmoment und Stiitzkraft mindestens an etwa 
ebensoviel Punkten Null werden. 

Die Integration der Potentialfunktionen ®,, wird mit Hilfe der Laplaceschen und Cauchy- 
Riemannschen Differentialgleichungen durchgefiihrt. Die Bipotentialfunktionen y,, kénnen wegen 
«x = 0 am freien Rand integriert werden. Die Bipotentialfunktionen V,, und U,, werden fiir die 
Integration durch x O, + @, dargestellt. Es ergeben sich fiir die einzelnen Funktionen folgende 
Ausdriicke: 


*Re } 
| (2 +953) (eff (9) =—a—» ; R/O 
: 2M alee a Re ge) 
i | (as t=») = =) if 1) pes) 


le tad (= ee HO) (2 | res , 


(44) 
fleets 
Rel» 
| | & t (2—9) ea ("al tO)| dy? = (1—y)r “int Ee |F(0), 


ite VO) f 

e| es a. \~ ReJ 
HO] ay = (1+ 9) 29 [PF (Lao) 3 : an ; 
y 
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Die dem Stiitzmoment M,, entsprechenden Ableitungen des Singularteils der BiegeinfluBflache 
wie auch der SchnittkrafteinfluBflachen sind am freien Rande gleich Null, ebenso auch die Inte- 
grale langs des freien Randes iiber diese Ableitungen. Die der Stiitzkraft Cy, entsprechendea Ab- 
leitungen der Singularteile der EinfluGflachen sind ebenfalls am freien Rande gleich Null. Auch 
die Doppelintegrale langs des freien Randes ttber diese Ableitungen mit Ausnahme des Doppel- 
integrals tiber den Stiitzkraftausdruck der BiegeeinfluBflache sind Null. Dieses letztere jedoch 
nimmt folgende Gestalt langs des freien Randes an: 


1 [( ([ee 
7] |[ae+e 


Diese Funktion (45) hat an der Stelle y = 0 (Abb. 5), also im Aufpunkt eine Spitze. Da die Uber- 
lagerungsfunktionen aber alle dort regular sind, lieBe sich mit ihnen weder das Doppelintegral 
iiber die Stiitzkraft noch die Stiitzkraft selbst ausgleichen. Jedoch ist die zweite Ableitung nach y 
(lings des freien Randes) von (45) gleich Null und man kann deshalb (45) tiberhaupt weglassen, 
was soviel bedeutet, da man diese Funktion von y wieder in Abzug bringt. Fiir den vorher voll- 
zogenen Eckenausgleich bedeutet das aber, das nun am freien Rand in der Ecke wieder die Sin- 
gularitaten dieser Funktion von y (37c) auftreten. Von den Mittelpunktsfunktionen besitzt aber 
ebenfalls eine, namlich U; in (43) eine solche Singularitat, so daB mit dieser Funktion also der 
Ausgleich in der Ecke und dann mittels den iibrigen Funktionen die punktweise Erfillung aus- 


°C de 


0 
45 
v) Ox dy? (45) 


(3 4 vy) arctg 2 = (3+) Fy. 


i\x—> 0 


gefihrt werden kann. 
Es lassen sich aber durch Summation mehrerer Eckenfunktionen auch zwei Funktionen her- 
stellen, die diese zusatzlich auftretende Eckensingularitat ausgleichen. Sie lassen sich in der Form 


(C,, ®,, sie Cm Vin) 2 


ee 


m 


ye = > (C,, ®,, tr Ke is) 
schreiben. Eine der Funktionen U;, y,, V, muf also hinzugenommen werden zu der punktweisen 
Erfillung, da ohne sie eine auch nur naherungsweise Erfiillung der Randbedingung der Stiitzkraft 
fiir die BiegeeinfluBflache nicht méglich ist. Die Bedingungen fiir die punktweise Erfillung und 
die letztgenannte Bedingung wurden in einem Gleichungssystem zusammengestellt. Zusatzlich 
wurde noch die Erfillung von Stiitzmoment und Stiitzkraft selbst im Aufpunkt gefordert. Die 


Ausrechnung des Gleichungssystems ergab die Konstanten in Tabelle 7. 


Tabelle 7 
m om K,, m 
m Biegung Moment Biegung Moment Biegung Moment Biegung Moment 

0 — 50,711831 | +117,912629 | +149,013065 | —220,301587 | — 78,726474 | —100,674663 — —_— 

1 — 82,131071 | +210,369150 | +161,587416 | —335,131495 | —135,288946 | +425,712801 | + 32,615249 | + 48,038623 
2 — 36,826453 | + 89,347665 | + 99,861082 | —255,495318 | — 24,627343 | + 33,127950 | + 39,371866 | —183,789778 
3 +112,893538 | —279,484262 | —347,760343 | +500,167824 | +173,038082 | +315,286788 | +198,316904 | —764,028817 
4 +193,414886 | —500,024982 | —391,060086 | +841,046416 | +338,221528 | —1067,272255| — 65,918293 | —152,271840 
5 + 88,892888 | —228,058863 | —235,960225 | +623,666896 | + 60,148573 | —107,028053 | — 95,376382 | +432,816272 
6 — 71,028788 | +198,371341 | +242,730292 | —316,258489 | —107,144064 | —342,644595 | —161,742830 | +617,024285 
7 —140,917163 | +373,835454 | +297,391938 | —679,226065 | —268,082244 | +857,406107 | + 33,521375 | +156,925867 
8 — 66,777736 | +182,692497 | +172,299786 | —478,067292 | — 49,158491 | +114,672256 | + 72,120753 | —318,115351 
9 + 6,224286 | — 29,940412 | — 38,231815 | + 13,219854 | + 10,079778 | +142,127478 | + 46,442769 | —175,459384 
10 + 29,489167 | — 82,711740 | — 67,909064 | +172,552829 | + 65,897854 | —215,846653 | — 0,359170 | — 53,743233 
11 + 14,855490 | — 45,449181 | — 36,210848 | +110,654057 | + 12,889068 | — 40,772153 | — 16,257077 | + 68,038273 
12 + 2,622795 | — 6,859296 | — 5,751199 | + 23,172398 | + 2,752678 | — 14,095008 | — 3,539888 | + 14,089600 


Es hat sich im Verlauf der Rechnung gezeigt, daB es fiir eine gute Naherung erforderlich ist, 


auBer in den Mittelpunkten der eingespannten Seiten auch je Halbseite noch mindestens in drei 
Punkten, beim freien Rand je Halbseite in vier Punkten, die Randbedingungen zu erfiillen. Bei 
einer noch weitergehenden punktweisen Erfillung und damit einer noch weiterreichenden Er- 
fiillung der Randbedingungen anderten sich die KinfluBflachen nur noch geringfiigig, woran die 
Konvergenz der Naherung zu erkennen ist. Der maximale Randfehler, bezogen auf die maximale 
Durchbiegung unter dem Aufpunkt, betragt in der vorliegenden Rechnung bei der BiegeeinfluB- 
flache am eingespannten Rand 4,2%, am freien Rand etwa 4%. Eine Verringerung des maximalen 
Randfehlers bei der BiegeeinfluBflache am eingespannten Rand von 7%, auf 4,2° und am freien 
Rand von 7,5% auf etwa 4% ergab eine Veranderung des Maximalwertes der BiegeeinfluBflache 
von weniger als 2°. Somit ist also eine Naherung von 2%, Genauigkeit erreicht. 


A 
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Die Auswertung der EinfluBflachen, die in Abb.9 und 10 dargestellt sind, hatte folgendes 
Ergebnis: Fiir eine Platte in Form eines gleichschenklig-rechtwinkligen Dreiecks mit eingespannten 
Katheten der Lange a und freier Hypothenuse ist mit u = 0,3 


Abb. 9. BiegeeinfluBflache. 


46 7 2B 19 
100 = OS z 010 
vies = — 


SS >< 


Abb. 10. BiegemomenteneinfluBflache. 
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1) die Durchbiegung w der Platte in Hypothenusenmitte unter einer dort angreifenden Einzel- 
kraft P * 
Wo» = 90,0182 a N° 


2) die Durchbiegung wy, und das Biegemoment M),, der Platte in Hypothenusenmitte unter 
einer iiber die gesamte Platte gleichmaBig verteilten Trae q 


q 
Wy» = 0,00124 at F , 


M, 


ogy = 0,00685 a4 = 

7. Zusammenfassung. Nach einer Einfiihrung in die Singularitatenmethode und die Her- | 
leitung der Fundamentalintegrale fiir die SchnittkrafteinfluBflachen wird fir den am freien Rand | 
liegenden Aufpunkt das Fundamentalintegral der BiegeeinfluBflache aufgestellt. Dann wird daa 
Verfahren von F’. Schultz-Grunow zur Ermittlung von EinfluBflachen umfangsgelagerter Platten | 
bei von Rand zu Rand wechselnden Lagerungsbedingungen g gezeigt, und insbesondere numeriseh | 
die Biege- und MomenteneinfluBflache einer Platte von ies Fon eines gleichschenklig-recht- | 
winkligen Dreiecks mit eingespannten Katheten, freier Hypothenuse und dem Aufpunkt in der 
Mitte der letzteren ermittelt. Dazu mu die konforme Abbildung des Dreiecks auf den Einheits- | 
kreis durchgefiihrt werden. 

Ks Pie sich, daB man fiir Platten mit freien Randern bei der Ermittlung von EinfluBflachen | 
nach dem Verfahren von F’. Schultz-Grunow gute Ergebnisse erhalt, wenn auch der Rechenaufwand | 
naturgemaB héher ist als bei anderen Lagerungsarten der Platte. Weitere Untersuchungen in | 
Fortfiihrung dieser Arbeit, die demnachst veréffentlicht werden, haben jedoch ergeben, da8 sich 
der Rechenaufwand bei Beriicksichtigung folgender Punkte wesentlich vermindern laBt. 

1) Die Potentialanteile des Fundamentalintegrals werden nicht mehr als Funktionen des 
polygonalen Bereiches definiert, sondern im Einheitskreis gebildet. Die so gebildeten Funk- 
tionen sind in den Abbildungspunkten der Ecken nicht mehr singular, so daB der sogenannte 
Singularitatenausgleich entfallt. 

2) Die zur Verwendung gelangenden Bipotentialfunktionen werden gebildet als 


y, =x-d®D /dx 


und mit den Potentialfunktionen ®, so zusammengefaBt, daB die erste Randbedingung des 
freien Randes exakt erfillt ist. 

3) Die drei Ecken werden dicht nebeneinander im Einheitskreis abgebildet, so daB die Er- 
fillung der Randbedingung des eingespannten Randes nur noch etwa in einem Punkt nétig 
wird, um eine entsprechende Genauigkeit zu erreichen. 

4) Wegen 2) sind die verwendeten Bipotential- und Potentialfunktionen in der zweiten 
Randbedingung des freien Randes orthongonal zueinander, so daf nicht mehr die Lésung eines 
Gleichungssystemes fiir die punktweise Erfillung der Randbedingung erforderlich wird, sondern 
die Konstanten durch Bilden der Minimumsbedingung ermittelt werden kénnen. 

Dadurch wird das Rechnungsverfahren sehr stark abgekiirzt und so die Ermittlung von 
KinfluBflachen bei Platten mit freien Randern der praktischen Nutzung zuganglicher gemacht. 

Herrn Prof. Dr. F. Schultz-Grunow, der die Anregung zu der vorliegenden Arbeit und wert- 
vollste Hinweise bei ihrer Durchfiihrung gab, sei an dieser Stelle der verbindlichste Dank des 
Verfassers zum Ausdruck gebracht. 


Aus dem Institut fiir Mechanik der Techn. Hochschule Aachen 
(Eingegangen am 22. Oktober 1956.) 
Anschrift des Verfassers: Dr.-Ing. Alfons Scherer, Koblenz-Oberwerth, Goethestr. 10. 
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Grenzen des Tragfliigelauftriebs 
Von H. B. Helmbold 


1. Einleitung. Die Grundlagen der quantitativen Theorie des Tragfliigels endlicher Spannweite 
wurden 1918/19 von L. Prandtl gelegt. Da zu jener Zeit die Auftriebszahlen praktisch verwendeter 
Tragfliigelprofile den Wert c, = 2 nie iiberschritten, war die linearisierte Form der Theorie fiir alle 
Anwendungen vollig ausreichend; gleichzeitig wurde die Auswertung der Theorie durch die Linea- 
risierung wesentlich erleichtert. Dieser Stand der Dinge blieb lange Zeit unverdndert, um so ist es 
ganz erklarlich, daB tiber der alltaglichen Anwendung und den Erfolgen der linearisierten Theorie 
die Tatsache, da die Theorie in ihrer urspriinglichen Form hinsichtlich der zugelassenen Zirku- 
lationsstarke keinen Einschrankungen unterliegt, allmahlich aus den Augen verloren wurde. In- 
zwischen aber ist die Erzeugung bedeutend héherer Auftriebszahlen durch Absaugung oder durch 
Rotation des Fliigelprofils in den Bereich des praktisch Méglichen geriickt. Fiir solche Probleme 
reicht die linearisierte Theorie nicht mehr ganz aus, und man muf deshalb zur urspriinglichen Form 
der Theorie zuriickkehren. Man hat nur einige Einzelheiten der Theorie Prandtls, die in seiner 
Darstellung infolge Spezialisierung nicht in Erscheinung treten, wiederherzustellen, die seither hin- 
zugekommenen Kenntnisse iiber den Endzustand des aufgerollten Wirbelbandes weit hinter dem 
Fliigel mitzuberiicksichtigen und einige naheliegende Folgerungen zu ziehen, um zu der Einsicht zu 
gelangen, da der Steigerung des Zirkulationsauftriebs durch die Endlichkeit der Spannweite 
Grenzen gezogen sind. Die Erzeugung ,,unmittelbaren“ Auftriebs durch Reaktion eines Blas- 
strahls bleibt hier auBer Betracht. 


2. Allgemeine Beziehungen. Weit hinter dem Tragfliigel ist das System der freien Wirbel um 
den Winkel 6 gegen die Richtung der ungestérten Parallelstrémung (umgekehrte Flugrichtung) ab- 
warts geneigt. Absolut — d.h. in einem Koordinatensystem, in dem die Luft in unendlicher Ent- 
fernung vom Wirbelsystem in Ruhe bleibt, — bewegt sich das Wirbelsystem senkrecht zu seiner 
Langserstreckung mit der Eigengeschwindigkeit w’ abwarts. Bezeichnet V den Betrag der Flug- 
geschwindigkeit, so ist w’ = V sin 6 und der mit der Vorwartsbewegung des Tragfliigels verbundene 
Langenzuwachs des Wirbelsystems in der Zeiteinheit gleich Vcosd. Bekanntlich findet man in 
groBer Entfernung hinter dem Tragfliigel die beiden Halften der freien Wirbelflache in zylindrischen 
Raumen, den Wirbelkernen, zu Spiralen aufgerollt. Die Absolutstrémung auBerhalb der Wirbel- 
kerne ist eine ebene Strémung mit Abwartsgeschwindigkeiten w parallel zu w’ und Quergeschwin- 
digkeiten v in Spannweitenrichtung. Wegen der angenahert schraubenformigen Gestalt der Wirbel- 
inien im Kerninnern herrscht dort auch noch eine stromabwarts gerichtete Lingsgeschwindigkeit u. 

Von einem zweiten, mit dem Tragfliigel fest verbundenen Koordinatensystem aus betrachtet ist 
die Strémung stationar und besitzt weit hinter dem Tragfliigel die Geschwindigkeitskomponenten 
V cos 6 + u in Langsrichtung des Wirbelsystems, v in Spannweitenrichtung und V sin 6 — w senk- 
recht zur Langsrichtung des Wirbelsystems nach aufwarts; wie schon angegeben verschwindet u 
auBerhalb der Wirbelkerne. Damit wird das Quadrat der resultierenden Relativgeschwindigkeit 
weit hinter dem Tragfliigel 


V2, =(Vcos6 + u)? + vw? + (V sind — w)*. (1) 


Der Druckunterschied zwischen dem Unendlichen stromabwarts (V.., p) und dem Unendlichen 
stromaufwarts (V, py) wird 


2 2 2 
| Le erat ee ee, V (w sin 6 — u cos 6) —o~ Se ; (2) 


Fir 6 = 0 geht dies in die von Prandtl zwischen seinen Gleichungen (65) und (66) angeschriebene 
Gleichung iiber. Wird nun der Tragfliigel von einer Kontrollflache umschlossen, die aus zwei parallelen 
und zur Langserstreckung des freien Wirbelsystems senkrechten Ebenen weit vor und hinter 
dem Fliigel sowie seitlich ringsum aus diese Ebenen verbindenden Abschlufbflachen in groBer Ent- 
fernung vom Fliigel besteht!, dann ergibt der Impulssatz als Komponenten der Luftkraft am Trag- 


1 I, Prandtl, Nachr. Ges. Wiss. Gottingen, math.-phys. K]., 1918 (1. Mitteilung), 1919 (2. Mitteilung). 
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fliigel einen ,,Langsschub“ parallel zur Langsrichtung des Wirbelsystems stromaufwarts 
fo) Fx 
P, = [ (p—po) dF + [ @ (V cos 6 + u) udF (3) 
und eine Querkraft senkrecht zur Verbindungsebene der Kernachsen (Schwerlinien der Wirbelkerne) 


nach oben ee 
P,=[o(Vcosd +u) wdF = [0 VweosddF + f puwdF. (4) 


Hier bedeutet dF ein Flachenelement der Querebene, Fx die Querschnittsflache der Wirbelkerne, 
und 9 (V cos 6 + u) dF ist die in der Zeiteinheit durch dF hindurchtretende Masse. (DaB der »Langs- 
schub“ einen negativen Betrag erhalten wird, ist in diesem Stadium der Betrachtung noch nicht be- 
kannt.) Durch Einsetzen von (2) in (3) erhalt man 


a 4 2 
Pe [ e(Yw sin 3 ees fesar. (5) 


x 


Auftrieb und induzierter Widerstand, die Luftkraftkomponenten senkrecht und entgegengesetzt | 


zur Fluggeschwindigkeit V, ergeben sich aus P, und P, in folgender Weise: 
A = P,cos6 + P, sind, W, = P, sind — P, cos 6. (6) 
Durch Elimination von P, wird der Ausdruck fiir den Auftrieb auf die Form von Prandtls 


Gleichung (68), 


AES SMES (7) 


gebracht. (Der von Prandtl angegebene Ausdruck A = 9 VA dy Wate 6 folgt dann aus unseren 
spiteren Gleichungen (15) und (20).) Die Kraft P, laBt sich als Zuwachs des Querimpulses der 
Absolutbewegung in der Zeiteinheit auffassen; dann bedeutet P./V cosd den Querimpuls 
bezogen auf die Langeneinheit des Wirbelsystems. Einsetzen von (4) und (5) in die zweite 
Gleichung (6) ergibt den induzierten Widerstand 

ioe) ae 


w= [oe Z* cos ddr + | g (uw sind— 5 cos). (8) 


2 


Dies unterscheidet sich von Prandtls Gleichung (67) durch das erste Glied im zweiten Integral. Prandtl 

konnte dieses Glied fortlassen, insofern bei der von ihm gemachten Annahme starrer Rotation in den 
PK 

Wirbelkernen das Integral i; uw dF aus Symmetriegriinden verschwindet. Daf die Kraft W; wie 

stillschweigend vorausgesetzt tatsdchlich einen Widerstand darstellt, ergibt sich aus der Energie- 


bilanz der Absolutbewegung. Die kinetische Energie dieser Bewegung nimmt in der Zeiteinheit um 
a 2 2 2 
| (Y e088 ee 


Die Zunahme der Gesamtenergie je Zeiteinheit ist wegen (2) 


Bie 
dF zu, wahrend die Driicke die Arbeit i} (p — po) u dF leisten. 


Fe 


K=V [e ete cos 6 dF + V | 0 (uw sind — "3 cos) dr. (9) 


Der Vergleich mit (8) zeigt, daB E = W, V die Arbeitsleistung eines Widerstandes gegen die Vor- 
wartsbewegung des Tragfliigels darstellt. Die Gesamtenergie bezogenauf die Langeneinheit | 


W; 


des Wirbelsystems ist ae ee 


3. Ein Gedankenexperiment. Es ist bisher nicht gelungen, den Aufrollvorgang der freien Wirbel- | 


flache von Anfang bis zu Ende exakt zu verfolgen. Wir werden deshalb so vorgehen, daB wir den 
Sonderfall eines Wirbelbandes betrachten, dessen Eigengeschwindigkeit langs seiner Quererstreckung 
(Spannweite) von Anfang an konstant ist, und das an der Verformung seines geradlinigen Quer- 
schnitts zunachst verhindert wird. Man mag sich das im Bilde eines in der Wirbelflache liegenden, 
schwerelosen Bambusrollvorhangs vorstellen, dessen starre Stabchen die Querschnitte des Wirbel- 
bandes verkérpern. Nach Ablauf einer sehr langen Flugzeit reicht die Wirbelflache soweit strom- 
abwarts, da sie weit hinter dem Fliigel einen ebenen Streifen bildet. Jetzt wird die Wirbelflache 
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freigegeben, indem — um im Bilde zu bleiben — der Bambusvorhang verdampft wird. Die grund- 
legende Annahme ist nun die, daB die stabile Endgestalt, der die freigegebene Wirbelflache zu- 
strebt, mit der Endgestalt der wirklichen Wirbelflache, die schon am Fliigel selbst sich aufzurollen 
beginnt, identisch ist. Da der Aufrollvorgang der hypothetischen Wirbelflaiche weit hinter dem 
Fliigel ohne Einwirkung auBerer Krafte erfolgt, miissen der Querimpuls je Langeneinheit und die 
Gesamtenergie je Langeneinheit ebenso wie die Zirkulation um die beiden Halften des Wirbelbandes 
erhalten bleiben. Dies sind die Invarianten des Problems. Auftrieb und induzierter Widerstand 
_ aber werden sich bei festgehaltener Zirkulationsverteilung lings der Fliigelspannweite wahrend der 
hypothetischen Verformung andern, da dieser Vorgang im fliigelfesten Koordinatensystem instatio- 
nar verlauft (der wirkliche Aufrollvorgang ist auf den Fliigel bezogen stationar). Da ferner alle not- 
wendigen Angaben iiber das unverformte Wirbelsystem bekannt sind, braucht man nur noch die 
Eigengeschwindigkeit des fertig aufgerollten Wirbelsystems zu bestimmen, um die Luftkraft am 
Tragfliigel zu ermitteln. Diese Bestimmung der Eigengeschwindigkeit ist beim gegenwartigen Stand 
der Theorie zwar nur angenahert méglich, doch bleibt der Fehler voraussichtlich ganz unbedeutend. 
Eigengeschwindigkeit und Neigungswinkel des Wirbelsystems sind, wie sich ergeben wird, vor der 
Verformung gréBer als darnach [siehe Gleichung (19)]. 


4, Das unverformte Wirbelband. Um die im vorigen Abschnitt vorausgesetzte Konstanz der 
Eigengeschwindigkeit langs der Spannweite des Wirbelbandes zu erzielen, muB die Fliigelzirkulation 
langs der Spannweite elliptisch verteilt sein: 


2 
y ist hier die Spannweitenkoordinate, b die Spannweite und J das Maximum der Zirkulation 
in Spannweitenmitte y = 0. Dann ist die Eigengeschwindigkeit des Wirbelbandes weit hinter dem 


Fligel 
L 


i eas a (11) 
und der Abwindwinkel 6, des Wirbelbandes ist dort durch 
sin 6, = ~2 (12) 


V 


bestimmt. Der Langenzuwachs des Wirbelbandes betragt je Zeiteinheit dx /dt = V cos 0), und die 
scheinbare Masse des ebenen Wirbelbandes betragt je Langeneinheit 


dm I : 
tilde Nae Pa Sane & 
7 ae be (13) 
Demnach nimmt der Querimpuls je Zeiteinheit um 
dm dx, Ls 
PE — es Tr Wy = 0 va b2 V cos Os °Wo (14) 
zu und betragt je Langeneinheit 
Pvp ie: Bee i io) ke = 15 
ATE RAE Pak eGrty=ov {Tay. de) 


Da die Wirbellinien im unverformten System samtlich einander parallel sind, verschwindet die 
Langskomponente u der Absolutgeschwindigkeit, und die kinetische Energie wird zur Gesamtenergie. 
Sie nimmt in der Zeiteinheit um 


dm dx, w? oA w2 
0 ~ dx, oa’ | £8 Oo 16) 
zu und betragt je Langeneinheit 
| are | eee ae or 17 
Vist tesa 2 @-3 fo. 7) 


5. Das Wirbelsystem im Endzustand. Das fertig aufgerollte Wirbelband besteht aus zwel Wirbel- 
kernen von entgegengesetzt gleicher Zirkulation J, . Wie A. Betz’ gezeigt hat, bleibt der gegenseitige 
Abstand der Schwerpunkte beider Halften des Wirbelbandes wahrend des Aufrollvorgangs konstant 


1 4. Betz, Z. angew. Math. Mech. 12 (1932) S. 164 = NACA-TM 713, 1933. 
20 
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gleich 
+ b/2 
il Tt 
= =—b. 18 
l=7 | Ty) dy = (18) 
—b/2 


Nimmt man vereinfachend erst einmal wie Prandtl [Gleichung (64) und (72)] an, daB die Zirkulation 
in zwei Wirbellinien statt in Wirbelkernen endlichen Querschnitts konzentriert ist, dann ist 


7 oabiate 

te sat a mb ~ 70 2) 
die Eigengeschwindigkeit des Wirbelpaars und die Stromlinien der Absolutbewegung sind Apollo- 
nische Kreise. Nun hat Betz vorgeschlagen, als Querschnittsform der Wirbelkerne Kreise anzunehmen, 
und W. Kaufmann! hat auf Grund der von Betz aufgestellten Satze iiber Schwerpunkte und Tragheits- 
momente der beteiligten Wirbelgruppen die Wirbelverteilung in den Kernen naher bestimmt, indem 
er die Spiralwindungen ebenfalls durch Kreise ersetzte. Damit ergaben sich die Wirbel in unmittel- 
barer Umgebung des Schwerpunktes und aufdem K ernrand auf den vorerwahnten Apolloniusschen 
Kreisen angeordnet. Das bedeutet aber, da alle Stromlinien der Absolutbewegung auBerhalb der 


Kerne ebenfalls Apolloniussche Kreise sind und die Gleichung (19) im Rahmen der Betz-Kaufmann- . 


schen Kreishypothese streng giiltig bleibt. Kaufmann hat seine Theorie doppelt nachgepriift: in- 
dem er die berechnete Druckverteilung mit MeBergebnissen verglich, und indem er numerische Uber- 
einstimmung der berechneten Energie mit dem genauen Wert (17) nachwies?. Der bei der zweiten 
Kontrolle beobachtete Fehler von 1% wird darauf beruhen, da er die von ihm gefundene Wirbel- 
anordnung zum Zwecke einfacherer mathematischer Behandlung durchweg durch die nur wenig ab- 
weichende Anordnung auf Apollonischen Kreisen ersetzte. Trotz des Naherungscharakters der 
Kreishypothese wird man annehmen diirfen, daB (19) auch fiir die Eigengeschwindigkeit spiraliger 
Wirbelkerne eine geniigend genaue Bestimmung liefert. Die Annahme wird weiter gestiitzt durch 
die numerischen Ergebnisse von F’. L. Westwater®, der den zweidimensionalen Aufrollvorgang fiir 
20 Einzelwirbel untersuchte, die nach gleichen Inkrementen von I’(y) iiber die Spannweite des Ellip- 
senfliigels verteilt waren. Den sehr giinstig ausfallenden Vergleich dieser Ergebnisse mit der Aus- 
sage (19) veranschaulicht Abb. 6 von Spreiter und Sacks 4. 

Da nun wie oben ausgefiihrt (Abschnitt 3) Querimpuls und Gesamtenergie je Langeneinheit fiir 
das unverformte und das aufgerollte Wirbelband iibereinstimmen miissen, erhalten wir durch Ein- 
setzen von (19) in (15) und (17), indem wir beachten, daB w’ = V sin 6 ist, zunachst 


Pz bot} Px jv Ge 2172 ot 
cosd cosd, © 8 pee erie (20) 
und 
Le OS En es : 
Vcosd cosd V cos 09 = 0 395 v* sin® 6 (21) 


und schlieBlich 


5 
V = 055 b? V2 sin? 6 cos 6 , (22) 
sowie wegen (7) 
3 
A=07 B V2 sin d(1—7 sin?) (23) 


Werden die letzten beiden Gleichungen durch ae o V? b* dividiert, so ergibt sich die dimensionslose 


Parameterdarstellung des Zusammenhanges zwischen Auftrieb und Widerstand (Abb. 1): 


A . 7. ae 

a gel hanr sin* 3), (24) 
0. 0 

W; Cwi m0 * 
ie ee ND (25) 
eb 

2 


W. Kaufmann, Sitzungsber. Bayer. Akad. Wiss., math.-nat. Abt., 1946, S. 109. 
W. Kaufmann, Ing.-Arch. 17 (1949) S. 187. 

F. L. Westwater, Aer. Res. Comm., Rep. and Mem. 1692 (1935). 

J. R. Spreiter — A. H. Sacks, J. Aero Sci. 18 (1951) S. 21. 


1 
2 
3 
4 
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worin A das Seitenverhaltnis bezeichnet. Der induzierte Anstellwinkel ist durch W JA =c,,,/¢,=tg a; 
bestimmt. Die reduzierte Auftriebszahl ¢,/A besitzt ein Maximum 

Calim me 2 

= —— = 1,90 fur snd =—— = 0,368, 26 

A 3/3 a V3 ey) 


dem eine reduzierende Widerstandszahl 


Cin my/ A 
relhe eer pans TS I 
ii lh oo 2,40 (27) 


entspricht. Der Auftrieb verschwindet fiir sin 6 — zie = 0,637. Wenn diese Ergebnisse wegen 


der leichten Unsicherheit der Zahlenkonstanten von (19) auch noch nicht ganz als endgiiltig ange- 


Abb. 1. Ideelle Polare des Tragfliigels mit elliptischer Zirkulationsverteilung. 


sehen werden diirfen, so geniigen sie doch jedenfalls fiir den Nachweis einer Begrenzung des Auftriebs 
von Tragfliigeln endlicher Spannweite. Experimentelle Auftriebswerte von Rotoren ohne End- 
scheiben bleiben durchweg unter der durch (26) angegebenen Grenze. 


Fir kleine Neigungswinkel 6 gehen die Gleichungen (24) und (25) in 


3 : 5 
“ = = sind und + = ig sin? 0 (28) 


iiber, woraus die bekannte Formel der linearisierten Theorie folgt: 


co 
oi i (29) 
In dieser Naherung wird 
sin 6 = 2 tg a; = 0,405 tg a; . (30) 


(Eingegangen am 24. Oktober 1956.) 
Anschrift des Verfassers: Prof. H. B. Helmbold, R. D. 2, Greencastle, Pa. (USA). 
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Allgemeine Stabilititstheorie der Schalen 


Von G. Schwarze 


1. Einleitung. Stabilitatsuntersuchungen aus dem Gebiet der Schalentheorie sind bisher fiir be- 
stimmte Einzelprobleme durchgefiihrt worden. So wurde u. a. die Beulfestigkeit der Kugel- und | 
Zylinderschale, der Kegelschale, sowie des Rotationsparaboloids in grundlegenden Arbeiten be- 
handelt. Ziel der vorliegenden Arbeit ist es, eine allgemeine Stabilitatstheorie fiir Schalen aufzu- | 
stellen. Die allgemeinen Differentialgleichungen des Stabilitatsproblems und das zugehérige Varia- | 
tionsproblem werden mittels der auf dem Gebiet der Schalenstatik tiblichen Methoden entwickelt. | 
Allgemeine Theorien der Schalenstatik liegen bereits vor. Es wurde hier im wesentlichen auf Ar- 
beiten von W. Zerna zuriickgegriffen 1. 

Die konsequente Anwendung der Tensorrechnung ermdglicht nicht nur, allgemeine Formeln von 
eroBer Durchsichtigkeit zu entwickeln, sondern erleichtert auch die Behandlung spezieller Probleme, 
was von umso groBerer Bedeutung ist, als man sich bei Aufstellung der Gleichgewichtsbedingungen | 
am verformten System in erster Linie auf die Anschauung stiitzen muBte bei Nichtanwendung des 
Tensorkalkiils. Die Lésung von Problemen, denen schiefwinklige Koordinaten zugrunde liegen, 
diirfte ohne Benutzung der Tensorrechnung auf groBe Schwierigkeiten stoBen. 


2. Bezeichnungen und geometrische Grundbegriffe. Zur Beschreibung der Formanderungen wird | 
ein raumfestes rechtwinkliges Koordinatensystem zugrunde gelegt. Zur Benennung der Schalen- 
punkte werden die Koordinaten vor der Verformung benutzt, so daB 0,, O, O, als kérperfeste Ko- 
ordinaten (Substanzkoordinaten) aufzufassen sind. Das Zahlentripel 0,, O,, O, charakterisiert 
demnach einen ganz bestimmten Schalenpunkt bei jeder beliebigen Formanderung. Der Ortsvektor 


t= 1 (0, O,) 
beschreibt die Schalenmittelflache. ©, ist der Abstand eines Punktes von der Schalenmittelflache 
und lauft im Intervall 


t t 
ee <ee 
pee Ole a 


Durch die Einschrankung von @, werden zwei weitere Flachen definiert, die Laibungen der 
Schale, deren Abstand, die Schalendicke t, gegeniiber den iibrigen Abmessungen der Schale als klein | 
vorausgesetzt wird. Die Flachen werden dabei als doppelpunktfrei angenommen. 

M=r+ O,a 
bezeichnet den Ortsvektor des Schalenraumes, wobei a den Einheitsvektor in Richtung der Flachen- 
normale darstellt (Abb. 1). 
Q, =+ t/2 Nach der Definition geht der Ortsvektor t iiber in 
M t=r-+ Do > 
0, --t/2 wobei 0g den Verschiebungsvektor der Mittelflache 
darstellt mit 
Do = v* A, +A (t= 52 )e 
Die a, sind die kovarianten Grundvektoren der Mittel- 
flache 
Abb. 1. A, = se ae. 


Allgemein sollen griechische Indizes iiber 1, 2 laufen, wahrend lateinische Indizes iiber 1, 2, 3 laufen. 


3. Geometrie der deformierten Schale. Werden entsprechend ¢ die auf die verformte Schale be- 
zogenen GréBen durch einen Querstrich gekennzeichnet, so folgt aus 


T=r1+b, (1) 
Cn = (6, ie uv iy bs r) A ae (b.6 v8 ae V, 4) a ? (2) 


LWA Zerna, Ing.-Arch. 17 (1949) S. 149; .A. E. Green u. W. Zerna, Quart. J. Mech. appl. Math. 3, (1950 
S.9; A. E. Green u. W. Zerna, Theoretical Elasticity, Oxford 1954. : Sy as : 
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0 
wobei (---) |, die kovariante Ableitung beziiglich der unverformten Schale definiert. Es ist 


3 le iis eA 0, 
ON far ses 


Ebenso lassen sich alle weiteren GréBen auf die entsprechenden GréBen des unverformten Systems 
und den Verschiebungsvektor zuriickfiihren. Dabei wird jeweils bei den linearen Gliedern der Ver- 
schiebungen abgebrochen, die zunichst bei der Stabilitétsuntersuchung nach der Gleichgewichts- 
methode ausreichen. 


Als MaBtensor der Flache gilt 


= = = 0 0 
@,8 = A,° Ag = 4,5 + , |g + ¥g|,—2b, gr. (3) 
Weiter ist dann 
= = 0 
a =deta,g=a(1+2 “|, —2 bio), (4) 
oe De 0 0 
ee al 2 a 2 be a (5) 
Ae eta (OF BT B A B yA A 
a’ =a, a*? = (63 —v? |, + Bh v) a + (bv + a4 v,,) 0. (6) 
Als Flachennormalenvektor ergibt sich 
a =a(b,,0% + %4,) a. air) 
Die Kriimmungstensoren der Flache und die Drei-Indizes-Symbole zweiter Art sind 
=> a 0 0 0 0 
bi = bag —v by bag + 04], bap + 04 |p bay £4 bag |p +” lap (8) 
yo a ? a ° oe 9 , 
bg = bg + v|,ga*” —v bj bj + 04 |, ba — 0* |, bj + ve a7 bog |y + 2 v bg, b%”, (9) 
Te 3 4 a? (o4 oO O o 
Ly = Tp, + v* |p, —v 3 |, — v2 bog by + v2 bg, bp — v,, 63 —y 4b, + v,, 2% bg,. (10) 


Der Ortsvektor des Schalenraumes ist fiir die verformte Schale 
R=R+BV mit R=r+O,-a. 
Die geringe Schalenstarke erméglicht, folgende Annahmen als hinreichende Annaherung fiir den 


Verformungszustand einzufiihren: 


1. Die bei der Verformung eintretende Verschiebung eines beliebigen Punktes sei durch die Ver- 
schiebung des zugehérigen Punktes der Mittelflache derart bestimmt, daf die Normalen zur Mittel- 
flache auch nach der Verformung zur neuen Mittelflache normal sind. 


2. Die Spannungskompenente normal zur Mittelflache sei so klein im Vergleich zu den anderen, 
daB ihr EinfluB auf die Formanderungen vernachlassigt werden kann. 


Damit ergibt sich der Verschiebungsvektor %§ in der Form 
XY = vy + O, w 
mit 
by =v,0%+va und wW=w, 0%, 
wobei tv die Drehung der Normalen angibt und parallel zur unverformten Mittelflache verlauft, wenn 
die Komponente in Normalenrichtung gegeniiber den iibrigen Komponenten vernachlassigt wird. 
Die tv-Verschiebung laBt sich durch die Verschiebung der Mittelflache ausdriicken: 


Q 
W, = —V,, — b, 0. 


Damit wird 

B =v, a* + va— Oz (v,, + 55 v,) a. 
Der Ortsvektor des Schalenraumes nimmt damit die Form an 

R=1 + O,a + v)— 3 (rig + 88 v,) 0° 


oder 


R=T+0,4. (11) 
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Daraus lassen sich weitere GréBen ableiten: 


Es ia oie mae 
Be = oe = te — Op BE a, = (02 — Oe BE) 
= OR ee 
687 On =, a. 
Als MaBtensor des verformten Schalenraumes gilt 
Bap = 45 —2 Os bap + O3 bs et (12) 
Bas = 833 — 9 > 833 = 1, 
entsprechend ¥ ee, 
gh = ao? © 26, bo? 3 OF be bo (13) 
gre = gl = 0, g8B=—1, 
und es ist ae 
p=/£ =1-20,7+ GK > (14) 
mit 


= i! = = — ee = = 
H= 7+), K=KH— DE. 


4. Der Spannungszustand der verformten Schale. Die Spannung am verformten System wird 
definiert als das Verhaltnis der auf ein Flachenelement wirkenden Kraft zu der Flache des unver- 
zerrten Klements. 


Es sei t der Spannungsvektor der verformten Schale. Aus dem Kraftegleichgewicht an einem 
Tetraeder folgt dann unter Beriicksichtigung der Spannungsdefinition 
tds = t, ds, , 

wobei ds und ds; die Flachen eines Tetraeders vor der Verformung darstellen. 
Mit 

ds; = Vg n; ds , 
n; als kovariante Komponenten des Normalenvektors zu ds, wird 

t=Jt, Vein; , 
wobei t eine Invariante ist (Abb. 2). 
Daraus folgt, daB sich t, Vgit verhalten muB wie ein kontravarianter Tensor. 


Ein Vergleich der Transformationsformeln fiir den Spannungs- 
zustand am verformten System einerseits und am starren Ele- 
ment andererseits liefert eine vollige formale Ubereinstimmung 
der einzelnen Glieder bis auf bestimmte Koeffizienten der 
Spannungskomponenten 


il 0 1 0 
Av) =1— — (ry, lh — by, v) baw. f,(v) =1— — (v, le — bog) , 
ayy Ago 


die bei der Betrachtung am verformten System zusatzlich auf- 
treten!, 


Abb. 2. 


~ Daraus folgt zunachst, da bei Beriicksichtigung der Verformungen die Spannungskomponenten 
t*8 keine Tensoren mehr sind. 
Durch Einfiihrung neuer SpannungsgréBen T%° ergibt sich aber wieder ein kontravarianter 
Tensor. Die GréBen 72° sind nach den Transformationsformeln definiert durch 


T%} =fAe)™, T=fo)t™, BW=fe)T, Tf =f) (15) 
mit 


fo =1 _ (vo is boo) , (16) 


wobei tiber groBe griechische Indizes nicht summiert werden soll. 
Bei Vernachlassigung der Verformungen gehen die GréBen T%° und T*? wieder in 1°? iiber. 


* Vel. A. Pfliiger, Stabilitatsprobleme der Elastostatik, Berlin 1950. 
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Mit Veg i — g, ergibt sich als Momentengleichgewichtsbedingung 


4; x aM = () . 
Daraus folgt 


Der Tensor 7'* ist demnach symmetrisch. 


Die einander zugeordneten Schubspannungen 7'* und T*! dagegen sind nicht gleich. 


Damit ergibt sich dann 
t Vgi* = THO, 


oder 
ee ioe. 
i — 1 zis (@ — 0,58) a, +1 wa 
ee iF 
Mit 
E aie = (6 — 0,18) 7 
wir 


1 
t; = —— (aif a, + 73 a). 


FAN 


Infolge der geringen Schalenstarke laBt sich der Spannungszustand insofern vereinfachen, als 
man auf die genaue Kenntnis der Verteilung der Krafte iiber die Schalendicke verzichtet und sich 
auf ihre Resultierenden und auf die resultierenden Momente beziiglich einer in der Mittelfliche 
liegenden Achse beschrankt. Wird normal zur Mittelflache ein Schnitt durch die Schale gelegt, so 
hat man, um die in der ganzen Schnittflache ubertragene Kraft zu bestimmen, die Summe der 
Spannungsvektoren zu bilden. Macht man sich von der Lange des Schnittes frei, indem man die in 
dem Schnitt iibertragene Kraft auf die unverformte Langeneinheit der Schnittlinie mit der Mittel- 


flache bezieht, so verbleibt eine Integration der Spannungsvektoren iiber die Schalendicke t. 


Die Gesamtkraft, die an der Flache 0, = konst. angreift, ist 
4 1/2 


== i 
te cnt 0dr, 
epee | oa, + HD 
—1y2 
und mit 
= Ta 
i, = Ne, 
a a** 
wird 
4+ 1/2 + 1/2 
N, 0, Vg o%2 d@? +. q ia Vg 73 dQ? , 
— 1/2 — t/2 
See co 
Aw, =a [ |/£ aeao +a | //% x40" 
aneorets 
— 1/2 er! 
Mit 
eet 
nee — es o%2 dO 
a 
—y2 
als Langskrafttensor und 
+42 
q = yz 7x3 dQ 


als Querkraftvektor wird 


Durch 
yanve Fe) nxe — Nee 
jaq* — Yee qe == Q« 


(17) 


(18) 


(19) 


(20) 


282 Schwarze: Allgemeine Stabilitatstheorie der Schalen Ingenieur-Archiv 


=a ; 
mit f@) = 12 Sn NS dl il bi v (21) 
ergibt sich E a =e 

M, = Nea, = 0% 0 
Mit N, = ya axe 1, folgen die SchnittgréBen 


il, =—— (nea, + 7a). (22) 
ar 
Infolge der Verschiebung der Kraft tritt ein Moment auf mit dem Hebelarm 0, a 
+ 4/2 
iit,= ee [ Vea x Gs¢a, + 73° a) 0, der. 
aa i 
Mit ' 
m, = M,, 
a aso 
wird daraus 
+ t/2 


Mit ed Ee 
ag = Nye ; E41 = &g, = 9 
(Oi OX Ap = Eo, A bei ' aay =o 
12 21 
und 


mee = | |/£ oy 0,46" (23) 


als Momententensor wird 


wobei a 
Pay _ fa 7700 —— WV og 
yamre = 7 mee = Mve (24) 
ist. Damit wird dann 
Ti, = = meee, a. (25) 
ak % 


/\ 


5. Gleichgewicht der deformierten Schale. Infolge der Beschrankung des Stabilitatsproblems auf’ 


die Ermittlung der fiir die Bemessung eines Systems maBgeblichen kritischen Belastungen, die durch 
die Bedingung des indifferenten Gleichgewichts geliefert werden, kénnen die Verformungen des 
Grundzustands vernachlassigt werden. 

Im Falle des indifferenten Gleichgewichtszustandes mu8 sich auch das durch bestimmte vir- 
tuelle Verriickungen 1,, 2, 0; verformte System unter dem vorliegenden Spannungszustand im 
Gleichgewicht befinden. Dabei wird unter einer virtuellen Verriickung die erste Variation eines Ver- 
formungszustandes verstanden. 

Im Anschlu8 an die iibliche Bezeichnungsweise soll der durch bestimmte, unendlich kleine Ver- 
anderungen des Verformungszustands erreichte ,,Nachbarzustand“ durch den Index I gekenn- 
zeichnet werden. wahrend der auf die Art seines Gleichgewichts hin zu untersuchende ,,Grundzu- 
stand“ den Index 0 erhalt. Der Nachbarzustand ist dann gleich dem Grundzustand plus dessen 
erster Variation, z. B. fiir die Langskrafte 


B 2 ae F : 

ny =n F(é) + ne? = njl + ae, 
wobei (. ~ .) andeutet, daf es sich um eine virtuelle Verschiebung bzw. die erste Variation des Grund- 
zustands handelt. 


Damit laBt sich das Kriterium fiir das indifferente Gleichgewicht folgendermafen formulieren: 
Kin Gleichgewichtszustand ist indifferent, wenn es mindestens eine spezielle Variation des Verfor- 
mungszustands gibt, nach deren Ausfiihrung die Gleichgewichtsbedingungen auch noch fiir den 
Nachbarzustand ohne Anderung der Belastung des Systems erfiillt sind!. 


1 Vel. A. Pfliiger, a. a. O. 
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Dabei ist in (26) 


Fv) =fe=1 a (vp ie — bgp). (27) 


Der Faktor F(v) resultiert aus der Tatsache, daB die Schnittgréen n° des Grundzustands innerhalb 
von (22) auf die Grundvektoren ag; des Nachbarzustands bezogen sind. Nach (22) ist namlich fir 
den Nachbarzustand 

1 


== [are 
Daraus folgt zunachst, daB fiir die n}° bezogen auf ag, kein Tensor mehr existiert. Durch Einfiihrung 


neuer SchnittgréBen n*8 ergibt sich aber wieder ein kontravarianter Langskrafttensor. 
Es ist 


Nar (ue Qor “te qr ty) % 


ye Ee inee. (28) 


Durch die virtuellen Verriickungen ¥,, 0, v, entstehen also in der Schale zusatzliche SchnittgréBen 
zu den schon vorhandenen SchnittgréBen des Grundzustands. Um festzustellen, ob die verformte 
Schale im Gleichgewicht ist, miissen diese neuen GréBen mit in die Gleichgewichtsbedingungen ge- 
nommen werden. Da aber die zusatzlichen Spannungen beliebig klein sein kénnen, so sind die neu 
in die Gleichungen eintretenden Glieder klein gegen die schon vorhandenen, und es miissen, wenn es 
schon auf Glieder dieser GréBenordnung ankommt, auch die Produkte aus den %,, v,, v3 mit den 
SchnittgréBen des Grundzustands beriicksichtigt werden, die dadurch entstehen, daB diese GréBen 
jetzt nicht mehr an dem urspriinglichen Schalenelement angreifen, sondern an dem durch zusatz- 
liche Formanderungen verzerrten Element. Es treten also zweierlei Zusatzglieder in den Gleich- 
gewichtsbedingungen auf. Da die Gleichgewichtsbedingungen auch ohne diese Glieder schon be- 
friedigt waren, so miissen diese fiir sich gleich Null sein, und da sich die SchnittgréBen linear und 
homogen mit Hilfe des Elastizitatsgesetzes durch die Verschiebungen ausdriicken lassen, so ergeben 
sich lineare homogene Differentialgleichungen fiir %, Y5, 03. 

Bei Aufstellung der Gleichgewichtsbedingungen ist zu beachten, da zwar Glieder von der 
GréBenordnung der linearen Anteile der zusatzlichen Verschiebungen beriicksichtigt werden miissen, 
jedoch Glieder héherer Ordnung vernachlassigt werden kénnen. 


Die Gleichgewichtsbedingungen 
werden fiir den Nachbarzustand in /m, ,m,,26%)g@’ 
Richtung der verformten Achsen auf- 4 /” 
gestellt. Das Gleichgewicht der Krafte 
ist nach Abb. 3 


( Ror Nz 200)’ 


Mra + PrV¥a=9- (29) -x, ae’ (MR, +N, dO) ae? 
Das Gleichgewicht der Momente be- fier? 
zogen auf den Schwerpunkt des 
Flachenelementes ist -MN,dQ? 


Mer aga +O. Mr, =9- (30) 


Die Krafte aus der auBeren Belastung 


(X, +M7d O) 20? 
a G 


sind dabei ); Ja. , wenn die Spannung 
als Kraft pro unverzerrter Flache Abb. 3. 
definiert wird. 


Werden die Vektorgleichungen in Komponenten zerlegt, so gilt mit 
Nr. = Ny a,+Q;7G und py =P; Oper Pr a 
Nis, Op + NF Gra + Orn A+ OF G4 + Yo (P; + Py a) =O, 
(Nix + NT? Dep — OF be + Va Py) d+ (Ny? bya + OF5 + Yap; a= 0. 
Jede Komponente muB fiir sich verschwinden. Also ist 
Nie + NG? Peo —O; bs + Yop; =0, 
Nr? bon + Ot. + Yap; =0- 
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Mit (19) und (20) ergibt sich dann 


haere ere a 4 _9 
Fey” «—~ 60, [f@r ™ Sa “+z Pr é 


1 of(v) 1 Bs 
f(r) bo. ny © ob pay le we (foe qr 4 ee Py 0. (32) 


Mit 
Wii Oa, 
ergibt sich aus (30) 
Min pine My 0; = 06 
und mit (20) und (24) folgt 


Ls eB) OL) eat ae ee eee 33 
flv) nT | 20. [fOo)P mM, f(r) qd (33) 
Durch Elimination der Querkrafte laBt sich die Zahl der Gleichgewichtsbedingungen auf drei redu- 
zieren: 
Dt) AO. Mahe yl fey eye red eee i 
fe) "t «~~ 20, [f@) ” Be reo) pt b * 605 [fO)]? bs +2 Py , (34) 
as, Ge afte) V\) gla 90) 1a 
fio Pear + Feey™ |e — 70 (265 OL ™r —~ 305 TOP 
_ Fo) 1 pe), F).AO) 1 on | fe 3 
20, Tree ™ 6 + 20, 20, Troe Ve Pre (35) 


Der Grundzustand soll im folgenden noch etwas ausfihrlicher untersucht werden. Und zwar sei 
zunadchst der Grundzustand als Biegespannungszustand angenommen, wobei also auch die Ver- 
formung des Grundzustands beriicksichtigt werden mu. 

Als Anfangszustand hat man also jetzt nicht den unverzerrten Zustand anzusehen, sondern den- 
jenigen durch gewisse Werte der d4uSeren Belastung hervorgerufenen Spannungszustand 0, bei dem 
ein unendlich benachbarter Gleichgewichtszustand I méglich sein soll. Die kovariante Ableitung 


Ce JI soll sich hier also ausnahmsweise auf die im Grundzustand verformte Schale beziehen. 
Die Gleichgewichtsbedingungen fiir den Grundzustand sind nach (34) und (35) 


Laat Fo) 1 wap mi pt Flv) 1 Galo EA 9 we 
Fe)" '*— 86, [/eoP tome =F) int be 805 (fe) +)2 er 
1 5 nz? sev kel) tel | ha. Ove) tine hay 
Fey bea 2? + Fy" low Fe) | 00, Fe Fes) a™> — 26, [Feoe™ | 

flv») 1 fv) f(r) 1 a 


eee aN rab : Bo a ne 
On [f(»)] yaa b+ 00x 00g [f(eo)]* LE oF Ve Po 0. 


Werden die quadratischen Anteile der Verschiebungen des Grundzustands vernachlassigt, so gelten 
die Gleichungen 


aA «9 as 
ne], — Ui mee |, + VS ph =O, (36) 


ao Ba" a 3 
Byars? + mi" |p. + VE pp =0. (37) 


Der eigentlichen Stabiltatsuntersuchung geht die Bestimmung der Lésungen des Grundzustands 
voraus. 


Die Gleichgewichtsbedingungen fiir den Nachbarzustand sind dann 


1 nt of (v) 1 net pt m°* pt Fe) 1 Bo res! 
FO™ — 30, Tor — Baym t |p + os 26, Tree +4 Piel 


NICLOMe Reet he 10) he 
{295 70 |' ™ —~ 20, LAr! |s 


Oe OM ORES 
. g v) | Op Fv) 
_ Fe) pa), Fo) Fe) _1 ps x 
20, wre ™ |° + 26, 26, [For ™ +e Bia Oe (39) 


of(v) 1 
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Wird durch die Beziehung (26) der Grundzustand in die Gleichgewichtsbedingungen eingefiihrt, so 
ergibt sich bei Vernachlassigung aller eee as Anteile der Verschiebungen 


ey iS of (v ey v) 
ee ee ee ee 
Lr is Of (v 
Bow nee T F() bo a Tox = me eee @ mé e = lean — i Pe ae 
_ F) Ba of(v) 1 
tor Mb [roy po] ats + Va 2h ee 


Werden durch (26) die Verschiebungen des Grundzustands eingefiihrt, die im allgemeinen recht 
komplizierte Funktionen von @, sind, wahrend andererseits die Krafte und Momente des Grund- 
zustands durch das Elastizitatsgesetz mit den Verschiebungen des Grundzustands verbunden 
sind, so ist ersichtlich, daB in bezug auf die praktische Brauchbarkeit die Zugrundelegung des 
Biegespannungszustandes als Grundzustand zu einem unméglichen Arbeitsaufwand fiihrt. Deshalb 
wird der Grundzustand in erster Annaherung durch einen Membranspannungszustand wieder- 
gegeben. Diese Annahme setzt beziiglich der Randbedingungen voraus, daS bis zur Erreichung 
der Stabilitatsgrenze die Verschiebungen der Schale ungehindert vor sich gehen kénnen. Eine Be- 
hinderung der Verschiebungen am Rande bewirkt ein Herabsinken der Tragfahigkeit, was zu einem 
Biegungsbruch vor Erreichen der kritischen Last fithren kann. 


Wird der Membranspannungszustand zugrunde gelegt, d. h. sind die 
queer) 


und wird die meatcarc des Grundzustands vernachlassigt, so vereinfachen sich die Gleichungen zu 


i yee tt ‘a2 
ot), F 75 nee re che lox fi. 022 — b, me le + Ve P—P=), (42) 
0 ae 
bo, Ree + 7655 bon 868 — Big Nee a A pag 2 P;— P=, (43) 
wobei f, , = f ie bezeichnet. Dabei wurde bereits beriicksichtigt, daB die Gleichgewichtsbedin- 


gungen fiir den Grundzustand fiir sich gleich Null sind. 
Der Belastungsvektor, bezogen auf die Einheit der unverformten Mittelflache, sei gegeben durch 


p = p*a, + pra. (44) 
Bei Deformation der Schale ist dementsprechend 
p =p*a, + pea. (45) 
Mit 
0 =i 
= (6; — v4 |, +0 dj) a, — (vbyg + 0,4) a 
und 


a=a-+ (vb, + %,,) a 
ergibt sich 


0 
p= p* (i —v' tv b2) + p3 (v% by + at V» 4) | (46) 
p® = p® — p* (0° bys + 1.) - 
Die BelastungsgréBen fiir den Nachbarzustand sind dann 
~ 0 = a. ~ 
py = pt (6, — |, + 5 Bf) + pe (ob; + a d,,) 5) (47) 
= + Po— Po (bag + 4) - 


Nach (22) und (25) ergibt sich eine Beziehung zwischen den Tensorkomponenten und den zugehérigen 
physikalischen Komponenten der rie ie 


= PTs _ aX —_ axe 
nee — z,, Nee = == = aC, Ip bap v) Vee do p iN, 
Peo ao 0 
met — 2 .M,, = + ee pee (% |, — age aso (Fam, (48) 
ay ay ay 


pe a, ee) Gt |M ps ae 
La y= My» y=} Bass 22 v) a2 aah 2 
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6. Allgemeiner Forminderungstensor. Das Quadrat des Abstands zweier benachbarter Schalen- 
punkte ist gegeben durch das Bogenelement und betragt nach der Verformung 


ds? = g,, dO* dOF + (dO%)? . (49) 
Der urspriingliche Abstand der beiden Punkte ist 
ds? = g,, d@* dO# + (d0*)? . (50) 
Im Falle einer Verzerrung ist 
ds? — ds? = (8,3 — g,) dO* dO* . (51) 
Dabei ist 
Vag eee) (52) 


ein Ma® fur die Verzerrung. 
Mit (12) wird dann 
Mes — OP mins mi 
Yap = 5 (ap — Gap) — 3 (Bap — bap) +o (Be Bp, — Bo Bp) i es 


Vos = 36 = Ys = O- 


7. Allgemeines Elastizitatsgesetz. Die SchnittgréBen fiir den Nachbarzustand sind definiert nach 
(17) baw. (23) 


$2 +42 
o & a & o 
nse = | |/£ off aoe, ie | //& of2 0, 46" 
—5/P) =e 
mit 
a@ _ 90 0) a8 
orn = (0g ay 0; bp) Tr. 
und 


Vf =1-20,H+ KG}. 
a -: 


Solange die Untersuchung auf rein elastisches Material beschrankt bleibt, und eine Umsetzung von 
Formanderungsarbeit in andere Energieformen ausgeschlossen wird, existiert stets ein Potential der 
inneren Krafte 2; als Funktion des Verzerrungstensors y;,. 


Es gilt fiir elastisches Material 


zik Oni 
© Oa? 


Das Potential 7; wird als Potenzreihe angesetzt in der Form 
06 905 a Ey ee ei 
wobei @, die Dichte des unverformten Kérpers ist. 
hing, = 0 ast 7.— 0, also J) 0, 
Durch Differentiation wird mit E**'" 4. Elrik — Biklr 


Oni 2 ai ee F 
= qik — Ei Biklr ‘ee 
a) Vik g * te Vir [ft 


Kir yj\== Osta" alec oh ti* bei y+ Co lhe 
Werden hoéhere Glieder der Potenzreihe vernachlassigt, so ist 
te Tik — pik 4 Biklray 
Fir isotrope Kérper gilt 


iklr i r ir : i r 
Bes = G (sits! +8 iene 


mit G als Schubmodul und i als Poissonscher Zahl. 


Fiir den Nachbarzustand gilt dann das Elastizitatsgesetz, wobei die Verformungen des Grund- 
zustands vernachlassigt sind, 
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werden die SpannungsgréBen ue eingefiihrt, so gilt, wenn nur lineare Glieder der Verschiebungen 
beriicksichtigt werden, nach (15) 


11 os - 22 “if 
x SET) 730 Dee bere 40 gate 
ee 8 te yt BE his 65) 
12 ae 21 zl 2 
SS Ee ie eee ast 2 
Seca a! aera o. . 


In dieser Form ist das Elastizitatsgesetz fiir Stabilitatsuntersuchungen bereits geeignet, wenn auch 
die Bedingung der Isotropie verletzt wird!. 


Die Einfiihrung des Elastizitatsgesetzes bedeutet eine an sich willkiirliche Definition, fiir die nur 
zwei Gesichtspunkte maBgebend sind. Das Gestz soll zundichst das Verhalten der benutzten Werk- 
stoffe méglichst gut wiedergeben und andererseits eine méglichst einfache Rechnung liefern. 

Aus diesem Grunde soll der Ansatz (55) vereinfacht werden zu einem Elastizitatsgesetz, indem 
die Verschiebungsgréfen f,(¢) gegen 1 vernachlassigt werden: 

Tr = 4+ Burrs, , t = 72 + Berg), , | 


12 ‘ a 21 __ f,(v) . 
ie + Bilr Vir > ty = 2 une Je Bary), F 


(56) 
AG) 
wobei T;' zwangslaufig festgelegt ist durch die Momentengleichgewichtsbedingung. 


Eine entsprechende Vernachlassigung in 7;" ist nicht zulassig, da die Verletzung des Momenten- 
gleichgewichts im allgemeinen eine nicht mehr widerspruchsfreie Rechnung bedeutet und auch 
zahlenmabig zu einer véllig falschen Beullast fihren kann. 


Ebensogut hatte das Elastizitatsgesetz auch in die Form gebracht werden kénnen 


22 


ty = + BM f, eee (57) 
12 f,(v) oe 21 ee 
Tye Te Virs ty = Ty BUSTY, 


Fiir den Ansatz (55) liefert die Anwendung der Momentengleichgewichtsbedingung fiir die zu- 
geordneten SchubspannungsgréRen einander véllig entsprechende Ausdriicke, wahrend sich im 
Ansatz (56) und (57) ein Unterschied zeigt, was einen Nachteil bedeutet insofern, als er eine Ab- 
hangigkeit des Gesetzes vom Koordinatensystem ausdriickt. Damit wird die Isotropie aber noch- 
mals verletzt. Ein weiterer Nachteil der Ansatze (56) und (57) liegt darin, daB das Gesetz von der 
Erhaltung der Energie nicht mehr erfiillt ist und daB die Anwendung des Prinzips vom stationaren 
Wert der potentiellen Energie nicht mehr méglich ist?. Dies ist allerdings ein besonderer schwer- 
wiegender Nachteil. Der Vorteil von (56) und (57) liegt andererseits in der Einfachheit, was bei der 
praktischen Rechnung zu beachten ist. 

Aus (54), (55), (56) bzw. (54), (55), (57) 1aBt sich ein Elastizitatsgesetz aufbauen, in dem die Vor- 


teile der vorher aufgestellten Ansatze ausgenutzt sind: 


11 y 22 a g2lr 
Pat eiry, Pa + Brg, | 
_ (58 
io tae 12Ir x 21 __ To Biles 
19 = am B 16: ae -f VI 4 
Ea Oe ae r = #6) 


Damit ist das Gesetz von der Erhaltung der Energie erfiillt, es 148t sich also das Prinzip vom statio- 
naren Wert der potentiellen Energie anwenden, andererseits sind die Ausdriicke fiir die Schub- 


spannungsgroBen wieder symmetrisch. » 
Das Elastizitatsgesetz (58) soll den weiteren Untersuchungen zugrunde gelegt werden. Mit 


(1 fir A=B, 
1 


48 — J FG) fiir Al ee i == ll. (59) 


1 


aa ihe a —sle B=2 
PAC) eae 


1 Vgl. A. Pfliiger, a. a. O. 
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ist dann 
lle Se ie EAT oe (60) 
Aus (17) baw. (23) folgt mit (15) 
ete 
7S | 2 flo) 02° dO, (61) 
42 
42. 
mie = | YE le) of? a0". (62) 
iP 
Mit (60) ergibt sich 
eye 
mee = | Vf) 0$ — 0, 6) 64? a5F + BPHY5,,) d® (63) 
—1/2 
und 
‘esa 
Ba | VE Solé) (68 — O, 08) (04 ve" + Bon F,,) Oy dO. (64) 
=r 
Mit | 
a2 Le 
| VE 0,08) 8 a6" = nse, 
#2 
und 


$2 
[ VE § — 6,18) x30 0, ao" = mge = 0 
= 42 


wird nach (26) und (27) 
+ 4/2 


nee — (§4P _ 1) nee + f, | Ve (68 — O, b§) Bates , dO, (65) 
Pee! 7 — 1/2 
ae | ys (63 — 5 bp) Br?u y,, Os dOS . (66) 


2 


Werden die folgenden Abkiirzungen eingefiihrt : 
1 
Heevr — a [Ql —p) (ax ace + qx ae) a 2 Ub axe at” ‘ 


i 
Asn = — [(1 — ps) (4 bH ber +. 3 aH bev by 4 3 ae” bay By + 4 bo» bee 


+ 3 aX” bey bY + 3 ace bx? by — 2 BT be ane — 2 BE ag” aH — 2 BS aH BBY 
— 2 bf ax” bf» — 2 b; ax” bee — 2 b; ace bx” — 2 bp ax” bu 
— 2 bg aFu b*) + 2 uw (4 bre bu” + 3 axe bey by + 3 at” bY be 
— 2 by are buv — 2 db; av bxe —2 bs ath bury — 2 b; al” §x8)] 
++ (b; bs + bt b2 63 — 64 b? 63) Haber , 

Deeur — [(1 — p) (a ber + ae” be + a%” bee + ae be”) + 2 yy (ae bY” + aH” bee)] 
— (b; 63 + bg) He6e , 


x 1 ~ 4 =~ 

Ony = oy (v,, y + vy he a or v) ? 

os — (841, bye oh ben be ee 
py lu “av jy “Aw v OY pa te Au Olle 
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3 Tues poe eae ae ee Oe, one, 0 r 
Bur =z @ |e Oya by + 0” |, by, by —v by b,, by + 0? by y |g by +0 [va 

0 

| 


_0 ; 7 a0 0 
+e liv Ole pars by by Diack v , 5, VN oe ly Br bya “= 0° eM bi. 
+2» b,, b4y ) . 
so folgt schlieBlich 
4/8 
g 5 
| Ve (53 as O; b5) bee Vuv dO 
—y2 
E xQUY y aQUY ax B 
= eee [Hera ay t+ (Axeus Buy + Deew @uy + Heeovr Bu») 75 (67) 


+ 1/2 
a E 
| Vz (63 =a @3 b3) BaBbuy Nae QO, dQ == es ye (Deer Bae au. Heeur @ 


12 


wig: (68) 


8. Vereinfachtes Elastizititsgesetz. Im Interesse eines konsequenten Aufbaues der allgemeinen 
Theorie wird eine grundlegende Vereinfachung im Elastizitatsgesetz eingefiihrt, und zwar werden 
in den Schnittkraften alle GréBen vernachlassigt, die proportional zur Biegesteifigkeit sind, wahrend 
in den Momenten nur die jeweils héchsten Ableitungen der w-Verschiebung beibehalten werden. 
Diese Vernachlassigungen werden durchgefiihrt in Anlehnung an die tiblichen Vernachlassigungen 
in der von Zerna aufgestellten Schalenstatik!. 

Durch die Annahme vom Ebenbleiben derQuerschnitte, die in Hinblick auf eine praktische Brauch- 
barkeit notwendig ist, werden von vornherein Glieder einer bestimmten GréBenordnung unterdriickt. 
Die in dem Elastizitatsgesetz genannten Glieder entsprechen aber dieser GréBenordnung. Es liegt 
daher nahe, auch diese Glieder zu vernachlassigen. Damit wird aus (67) und (68) 


+ #/2 


| y= 0s — 18 G,) Beeu" 5, dO® =~ | > Heew |, +4, |, — 2 Byy 0) t 
50n 
[ VE 08 0) BPH, O40 = —  F Ho0H sp. 
42 
Damit me sich folgendes Elastizitatsgesetz: 
Beg) ese eine (69) 
me SS eo os i) a 2) Fach gs (70) 


Dabei ist allerdings das Momentengleichgewicht um die Schalennormale nicht mehr exakt erfiillt. 
Die Abweichung ist aber von einer GréBenordnung, die den vernachlassigten Gliedern im Elasti- 
zitatsgesetz entspricht. Damit kann die Abweichung in der Momentengleichgewichtsbedingung 
ebenfalls vernachlassigt werden. 

Fiir das Beispiel einer auf Torsion beanspruchten Zylinderschale wurde dies vereinfachte Elasti- 
zitatsgesetz benutzt und es zeigte sich, dafs das Gesetz dort zu einem brauchbaren Ergebnis fihrt. 
Ob die Vereinfachung im Elastizitatsgesetz allgemeingiiltig ist, bedarf allerdings noch einer weiteren 
Untersuchung. 


9. Anwendung der Energiemethode. .Mittels der Energiemethode werden die Gleichungen der 
Gleichgewichtsmethode kontrolliert und gleichzeitig wird ein Variationsproblem gewonnen, das den 
Ausgangspunkt fiir ein wichtiges Naherungsverfahren bildet. Die virtuelle Arbeit der inneren 


Krafte ist 
~ sik «. — 
At = — | [ [Bare ao dO? do? = — ec Vz dO dO? dO. (71) 


Fe me 
: Q 
1 Vgl. W. Zerna, Ing.-Arch. 17 (1949) S. 149. 


Mit 
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wird dann 


At =— | Jf a* 0 /g dO dO? do. (72) 


Die GriBe %* wird als vollstandiges Differential der Funktion 7; vorausgesetzt; 7, soll nur ab- 
hangen vom Verzerrungszustand, vom MaBtensor des verformten Schalenraumes, sowie von Gré8en, 
die die physikalischen Eigenschaften der Schale charakterisieren, im Wesentlichen alsoz, = 7; (7;;)- 


Es ist 


Diese Bedingung stellt die Definition eines ,,elastischen Kérpers“ dar; 7; ist das elastische Potential 
bezogen auf die Volumeneinheit, 7; ist die erste Variation des elastischen Potentials. 


Unter dieser Bedingung ist 
At = A, =—TT, =— ||] %,0 Wg dO do? do* = — [f[n,0 Vg dO1 dO? d@® ~— (73) 


Uv 


mit //, als elastischem Potential der gesamten Schale. 


Mit 0 Vz == (0 Je (Satz von der Erhaltung der Masse) wird 
IT, = | [ {7:0 yg dO! dO? dO? . (74) 
Ein Potenzreihenansatz fir 7, liefert 
1 ep 
09%; = = tO" Veo 15 tet Vas Be 
Fir den Nachbarzustand gilt entsprechend 
I ~ eee: Boot 
QM; pgs VB a +17" fol?) Oy oe OS 


und durch Einfiihrung von (60) ergibt sich 


1 e 1 en 
Qo Mri =aao (lL + fg°0*”) TH? Yap Sap Oe Yur V«B ° (76) 
Mit | 
+2 
Ss aoe ey OGRA 
| Via Vas ae nee Vap 2 (77) 
Hp: 


wobei y,3 den Verzerrungstensor der Mittelflache darstellt, folgt fiir das elastische Potential be- 
zogen auf die Flacheneinheit 


+ e2 = + 4/2 
ed@3 — V2 ! ete ra ee: 
ae Yodo =F + fad) Ra mge +1 | VE BHM Fei dO. (18) 
ey — t/2 


Das Potential der inneren Krafte ist nach (74) 


+ 4/2 
i ‘ me 
LT; ; =F { {a + fp 04) ene? + | fe Ber y Vay A> 
— 1/2 


a dO' d@?. (79) 


Die praktische Berechnung von J/,,; wird noch erheblich vereinfacht, da die exakte Beriicksichtigung 

der Verformung des Grundzustands fortfallt, und infolgedessen Verformungsgréfen in jedem Fall | 

nur noch bis auf quadratische Glieder berechnet zu werden brauchen. | 
Mit den Abkiirzungen 


1 
a 2 (a jt) (ax af” + a2” abe) + 2 uy a%P an) , 
Buy 1 ‘ y 
HOP rece — f) (3 ax" bFY b, + 3 a® be” bY + 4 bee BBY + 3 as BBY bg 


+ 3 abu bor bY + 4 be bP) 4.2 (3 a ber BY 43 amr bor BE + A ba 7) 
Capuv = (L—p) (ae 58” + ab” baw + ax? bPe 4 abu bor) 4 2 ph (ax? bv + guy Boh) , 


SaBuy al Xx B uy ? 
NoBuv = Xp Ony ++ Cen Ove = D9H- Keep 6 oe ; 


PaBuy = Xy pg Bue ae Bin» Pag -+- Ong Vy —2 I &yp Oyy — 2 HS,; Onp ++ K &yp Sis : 
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geht das elastische Potential iiber in 


~ 


1 £ E 
Hy; age [I yrs’ Yop (1 + fe os) 7 [7 Eeoue Hobo ¢ 
= (ene enue oF Nat ae Gabur Si Pa Bur Hise?) a — er (80) 


Damit ist das Potential der inneren Krafte 7,; bestimmt, und ,, enthalt den Anteil, den die 


bei der Verformung neu auftretenden SchnittgréBen leisten, sowie den Beitrag der GréBen n*? des 
Grundzustands. 


Wird mit mit a* die virtuelle Arbeit der auBeren Krifte pro Einheit der Mittelflache bezeichnet, 
so gilt bei einer Variation des Nachbarzustands 


Qa — pr Ov, ail Pi OV, ae Pi Ov, . (81) 
Insgesamt wird die virtuelle Arbeit der auBeren Krafte 
Aj, = || at, Jad do, (82) 


wobei pj und p; die Tensorkomponenten der auBeren Belastung am verformten System sind und 
zwar bezogen auf die unverformten Koordinaten. 


Unter der Voraussetzung, daf es eine Funktion gibt mit 
d hg ata p) 


existiert fiir die 4uBeren Krafte ein Potential mit 


IT,, = | { =, ¥adO' dO@2. (83) 
Damit ]a8t sich das Gesamtpotential aufstellen: 
a Used 
HM, =; + Wy, =1h + 611 “rege thot. (84) 


wobei der eine Querstrich tiber dem Variationszeichen andeutet, daB es sich um eine spezielle Varia- 
tion handelt. 

Die erste Variation der potentiellen Energie des Grundzustands muB fiir jeden Gleichgewichts- 
zustand verschwinden. 


Fir den Fall des indifferenten Gleichgewichts ist somit 
i Lag a 
Oo = = (0?) a0). (85) 


Es muB demnach die erste Variation von mindestens einer speziellen zweiten Variation des Gesamt- 
potentials bei ungeadnderter Beslastung verschwinden. 

Die Eulerschen Gleichungen dieses Variationsproblems liefern dann drei Differentialgleichungen 
fiir die VerschiebungsgréBen 0, vy, 0, Andererseits bildet (85) den Ausgangspunkt fiir das Ritzsche 
Verfahren zur Lésung eines speziellen Stabilitatsproblems. 

Es wurden nach der oben entwickelten Theorie verschiedene Probleme durchgerechnet, wie 
Zylinderschale unter konstantem Druck sowie unter Torsionsbeanspruchung, ferner Kegel-, Kugel- 
und Rotationsparaboloidschale unter Eigengewicht. Es zeigte sich dabei an allen Fallen véllige 
Ubereinstimmung mit den nach bisher tiblichen Methoden abgeleiteten Gleichungen. Das Beispiel 
eines hyperbolichen Paraboloids zeigte dariiber hinaus, daB obiges Verfahren prinzipiell auch auf 
Probleme mit schiefwinkligen Koordinaten anwendbar ist. 


(Eingegangen am 31. Oktober 1956.) 
Anschrift des Verfassers: Dr.-Ing. G. Schwarze, Oberhausen-Sterkrade, Kirchhellener Str. 8 A. 
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Nichtlineare Schwingungen bei zwei Freiheitsgraden 
I. Mitteilung 


Von C. Schmieden 


1. Finleitung. In zwei Noten im Ingenieur-Archiv hat Th. Péschl versucht, im Falle eines Doppel- 
pendels mit endlichen Ausschlagen fiir die periodischen Lésungen — von ihm Hauptschwingungen 
genannt — die Beziehung zwischen Frequenz und maximaler Amplitude zu bestimmen?. Leider 
enthielt sein Ansatz einen prinzipiellen Fehler, wodurch seine Ergebnisse hinfallig wurden *. Das 
Problem blieb also offen; denn auch die Ergebnisse, die G. Bradastilow® in seiner Dissertation erzielte, 
zeigten nur, daB die periodischen Lésungen der linearisierten Differentialgleichungen — die Normal- 
schwingungen der Theorie der kleinen Schwingungen — fiir hinreichend kleine Amplituden auch fiir 
das nichtlinearisierte System existieren. Zudem wurde dabei der eigentlich interessante Fall aus- 
geschlossen, da die Frequenzen von Grund- und Ober-Schwingung in ganzzahligem Verhaltnis 
stehen. 

Da das System der Bewegungsgleichungen des Doppelpendels bei endlichen Ausschlagen auBer- 
ordentlich kompliziert ist, soll hier ein anderes, mathematisch einfacher zu diskutierendes System 
untersucht werden; es ist so aufgebaut, daB es fiir kleine Ausschlage in das System des Doppelpendels 
gleicher Fadenlangen, aber beliebigen Massenverhaltnisses iitbergeht, auBerdem aber auch die Ver- 
allgemeinerung der bekannten Differentialgleichung von Duffing auf zwei Freiheitsgrade darstellt. 
Fiir dies System werden wir die periodischen Lésungen fiir endliche Ausschlage aufstellen, und zwar 
auch in den Fallen, wo die Frequenzen fiir Grund- und Oberschwingung in ganzzahligem Verhdltnis 
stehen. 


2. Die periodischen Lésungen im allgemeinen Fall. Wir legen das folgende System von Differen- 
tialgleichungen unseren Untersuchungen zugrunde: 


I+ep +o29 + por =0, 
P+ $t0ip+ pore =0. 


Setzt man hier § = 0 und gleichzeitig w? = g/l, so erhalt man die Gleichungen fiir die kleinen Schwin- 
gungen des mathematischen Doppelpendels gleicher Fadenlangen 1, mit der oberen Masse m, (Aus- 
schlag ?) und der unteren Masse m, (Ausschlag y), wobei 


(1) 


Dietpone a am 2 
ao nts (O25 721) e. (2) 
Fir die Grundschwingung erhalt man in diesem Fall 
2 
Oh = p49 Gg =A sino, t, “0 =]= aA sn@,, ft. (3a) 
ebenso fiir die Oberschwingung 
2 
Oe gp=Asino,yt, P=—aA sin Wat. (3b) 


Wir wollen nun die periodischen Lésungen des Systems (1) durch Potenzreihen nach dem als klein 
vorausgesetzten Parameter f darstellen. Dazu bringen wir das System (1) zunachst in eine andere, 
einfachere Form, indem wir setzen 
Ss db , d y 

T= Oy) Woe: FU T=¢- (4) 
Damit geht (1) iiber in 

(+1) (0" +0°9") +9488 =0, 

: +7) 0" + 9") +e+ Bg =0. 


* Th. Péschl, Ing.-Arch. 20 (1952) S. 189 und 21 (1953), S. 396. 
* Th. Péschl, Berichtigung. Ing.-Arch. 22 (1954), S. 294, 
° G. Bradastilow, Math. Ann. 116 (1939) S. 181. 
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Da fiir kleine Schwingungen Grund- und Oberschwingung durch Vertauschung von a mit — x inein- 
ander tibergehen und der Parameter a nur quadratisch in (5) eingeht, geniigt es, die Grundschwingung 
weiter zu verfolgen; die Oberschwingung erhalt man durch Vertauschung von a mit — q. Wir 
machen nun fiir ?, y und x die folgenden Ansatze: 


1 
Lteapigtak tmp te, 
yg=Asint+@,sn37r+-.--, (6) 
O = @,sint + O,sin37-+---, 


wobei wir die ®, und @, in Form von Potenzreihen nach § ansetzen: 


O,=aA+BpInt+Pay+---, 
7 
OP go POP Uey D, =Pvo + Br + Boos + °°: (y 2 3). ©) 
Ist das Verhaltnis von @ 9 ZU @g9 nicht ganzzahlig, so beginnen die Entwicklungen von ®,, 1 und 
Oon+1 mit Bp” (s. u.). 


Wir berechnen jetzt zunachst #* und gy soweit, wie wir es auch fiir die Ausnahmefalle gebrauchen 
werden, und erhalten 


49° — Q, snt + Q,sin37+Q; sin571 +Q,sin 71 + Q, sin 97 (8) 


mit 
Q, = 3 (Of —O7 O; + 2 O, OF — 2 O, O; O; + OF O; — OF O,) , | 
0, = — 63 + 3 (2 G2 O, + O8— G20, + 20, O, O, —2 O, O, O, — O20) , 
0, =3 (— G20, + O, 02 + 2 C20, — 62 0, +2 0, O, 0, —2 0, 0,0, +2.020,),$ (9) 
Q, = 3 (— O, 3 — O2 6, + 2 6? O, — O? O, + 20,0, O, + 20, O; O, + 2.030,), 
Q, = — G3 + 3( G20, + 2020, 20, 0,0; +20, 0,0, +2 020,) . 


Die entsprechenden Formeln fiir g* erhalt man, wenn man in (9) 0, durch A und die iibrigen 0, 
durch @, ersetzt. 


Setzen wir nun zunachst voraus, daB 0, und ®, von der Ordnung f ausfallen, so erkennt man aus 
(9) sofort, daB dann die ®2,,,, und @2,,4; von der Ordnung f” ausfallen. Setzen wir demgemaB 


Oon 44 = lise Oon 41 > Don 41 = BS P2n+1 ) (10) 


so reduzieren sich die Q,,, wenn man Terme von der Ordnung f° vernachlassigt, erheblich, und man 
erhalt durch Koeffizientenvergleich in bezug auf die sin 7, sin 3 t usw. die folgenden Gleichungen: 


— (1 +7) (0; +0? A] + 0; + 7B [Of —B O29, + 26? O, 8] = 0, 


—(l+%) [0,44] +4+ 4B [4 42g, + 26? Agi] =0; 


—9 (1 +) [B, +0°g.] + 9s +4 [— OF + 68 OF 8,—3f* O28] =0, 


—9 (1 +x) [9s +96] +¢+ ql 4? + 6B 42g, — 3B A?gs] = 05 


— 25 (1 +) [9 + 029] + 9; +7 [— 7B, + BO, 82 +26 OF8,] =0, 


—25 (1+ %) [95+ 96] + +71 4295+ PAG +26 42y] =0; 


— 49 (1 + x) (8, + 0? @,] +o +7 [—O, B— OF 8] =0, 
3 
—49 (1+) [+o] +e + 74g — 4*95] =. 


alle 
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Durch Koeffizientenvergleich der niedrigsten Potenzen in B folgen daraus zunachst die Gleichungen | 
1 
(8 — a) Dos + 9 0? Gs + Go? (1 + «) A8 = 0, 


9 Bq + (B—2) 9m +4 (1 +2) 4° = 0; 


(12) 
Dy —A (1 +a)? +702 (1 +0) 42 =0, 
By +A (1 +a)?x,—+ (1 +a) 4? =0 
mit den Lésungen 
= tits zO—8a+0%), P= 59 (8 —a —9 04), a 
=> A? ee Dy => A8 (1 +a) (1—o%). 
In der n&chsten Ordnung entstehen die Gleichungen 
a yy — Aa (1 +0)? %— (1+) 4 Oy +7 APa? (1 +) [3 Pn — Onl= 0, 
Dig +A (1 +0)? + (1 +0) On + 42 (1 +0) 91 = 95 
(8—2) Dyg +902 p59-+ + A202 (1 +a) Oy —> A202 (I +0) Dg + 9(1 +01) 2 (On +0 Pat) =O, a 


9 Bop + (8—a) ? —> 4? (1 + &) Ga +9 (1 + 4) % (Ig, + Pa1) = 9; 


(24 — a) Bsp + 25 0259 + + A? (1 +0) Dy, = 0, 
3 
25 Bs ++ (24 — a) P59 + 7 A? (1 + 4) qs, = 0. 


Um unnétige Formelschreiberei zu vermeiden, geben wir hier nur die Werte von ),, und x, explizit — 
an: 


Liawe!) (ei 28 1— a) (8 — 9 
Dre = 395 4° r = = eon +a) [ Wha re ioe (15) 
9 l1—at , 3 (1 ++ «?)? (8 — a) I ee 9 
= Se a te ae se 


wahrend wir fiir die iibrigen vier GréBen nur bemerken wollen, daB bei qs. und Os jeweils ein Term 
auftritt, der dem Faktor (4 — 5a) bereits quadratisch im Nenner enthalt, wogegen beig,, und Os» 
die Faktoren (4 — 5 a) und (12 — 13 «) jeweils im Nenner auftreten. Schreibt man die Gleichungen 
fiir x, und 0,3 an, so sieht man, da in ihnen bereits 75. und qs. auftreten, also x, den Faktor (4 — 5 a) 
ebenfalls quadratisch im Nenner enthalt. Bei x,, in dessen Bestimmungsgleichung auch gy; und ?55 
eingehen, tritt erstmalig der Faktor (12 — 13 «) im Nenner auf, und so geht es fort. Dies Verhalten 
zeigt jedenfalls, daf fiir ein x in der Nahe dieser kritischen Werte die Konvergenz unserer Entwick- 
lung sowohl fiir die Amplituden wie fiir die Frequenz aufhéren mu8. Dies gilt allerdings nur fiir die 
Grundschwingung, da ja fiir die Oberschwingung « durch —a zu ersetzen ist und damit an den kriti- 
schen Stellen keine Singularitaten auftreten. Fiir hinreichend kleine f ist also fiir die Oberschwin- 
gung die Konvergenz unserer Entwicklungen gesichert, es sei denn, a sei sehr klein. 

Denn beide Entwicklungen versagen auBerdem noch fiir sehr kleine Werte von a, da bereits in x, 
und #1, Terme mit « | auftreten und, wie man leicht sieht, daraus in den Termen héherer Ordnung 
Glieder mit héheren Potenzen von a+ resultieren. Wir haben also zunachst folgendes Ergebnis: 

Nur im Bereich B <a <4/5 —f ist die Entwicklung fiir die Grundschwingung fir hin- 
reichend kleine f sicher konyergent. Die Ungleichung fiir « ist dabei als Ungleichung hinsichtlich 
der GréBenordnung in f zu verstehen. Bei der Oberschwingung gilt nur die untere Schranke. 

Wir wollen jetzt die physikalischen Ursachen fiir das Aufhéren der Konvergenz der Entwick- 
lungen aufzeigen und untersuchen zu dem Zweck zunachst die Singularitat fiir « = 4/5 naher. 
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3. Der Fall a = 4/5. Setzt man den Wert a = 4/5 in die Frequenzgleichungen der kleinen Schwin- 
gungen (3) ein, so sieht man sofort, daB sich diese Frequenzen jetzt wie 1:3 verhalten, also das zu 
sin 3t proportionale Glied der Grundschwingung in Resonanz zu dem Anfangsgliede der Ober- 
schwingung steht. Bei den kleinen Schwingungen kann man also Grund -und Oberschwingung mit 
beliebigen Amplituden addieren und erhalt wieder eine Normalschwingung. Fiir das nicht-lineare 


_ System entfallt natiirlich diese beliebige Uberlagerung; aber dafiir treten, wie wir zeigen wollen, an- 


stelle der einen periodischen Lésung, die fiir « < 4/5 die Fortsetzung der Grundschwingung bei 


- B=0 darstellt, in der Umgebung des kritischen Wertes von a drei verschiedene periodische Lésungen 


auf, bei denen die Amplitude des zu sin 3 t gehdrigen Termes von der gleichen GréSenordnung wird 
wie die Amplitude des Grundtermes. Es liegt aber keine Verzweigung der Lésung vor! Um die 
periodischen Lisungen nahe bei « = 4/5 zu untersuchen, setzen wir 


6 
ate, ltu=stmp 4 
; (16) 
0,=~A Ogee: sas BO Van ct P Va he = Cay ee Pay os 


Da jetzt ®, und OQ Terme endlicher Ordnung enthalten, werden daraufhin ©, bis O, von der Ord- 
nung f, so dafs nunmehr in den Q, der Gleichungen (9) alle hingeschriebenen Ausdriicke beriicksich- 
tigt werden miissen, wenn man die Lésung bis auf Glieder in f? einschl. erhalten will. Um aber das 
charakteristische Verhalten der Lésungen zu erkennen, geniigt bereits die Kenntnis aller zu B pro- 
tionalen Glieder, so da wir uns auf diese Naherung beschranken kénnen. Aus (11) erhalten wir jetzt 
als Gleichungen nullter Ordnung in f die folgenden: 


ae bg) 4 +5 4=0, is F1jA+4=0, 


= 16 
—5 (Ban + 35 Poo) + Pu = 0, — 5 (P39 + P30) + Pao = 9- 


(17) 


Die beiden ersten sind natiirlich, wie es sein muB, identisch erfiillt, wahrend die beiden letzten iiber- 
einstimmen, und liefern 


4 
3 Dy + 4 P59 =0, Dsq = — = Psp - (18) 


Wiirde man jetzt A = 0 setzen, so erhielte man bei willkiirlichen g,,) natiirlich gerade die Ober- 
schwingung. Es muf also A = 0 beibehalten werden. 

In der nachsten Ordnung erhalt man fiir x,, 3,1, 031, Ps, die folgenden Gleichungen, wobei bereits 
Os durch ~z) ausgedriickt ist: 


4 36 5 3/43 

5 Ou — xm A—ZtA r(5] (A? + A? qs + 2 AG) = 9, 
5 9 3 

9 Ou 541A 4 7 (— 4? + A? Pa — 2 A Yip) = 0 5 


(19) 
9 25 4 , 
5 Bey + 4 Qs, == %1Ps0— | § Pao 5 (4 + 6 A? qq + 3 GBp) » 


9 il 
5 Ps, + 4G) = ad hae aa (A® — 6 A? day — 3 Yip) - 


Sollen die beiden letzten Gleichungen keinen Widerspruch enthalten, so miissen ihre rechten Seiten 
iibereinstimmen. Das liefert aber eine Gleichung fiir das noch unbekannte 5) und x, mit & als Para- 
meter. Wir haben also jetzt drei Gleichungen fiir die drei unbekannten o,, 0; und Qo. Eliminiert 
man aus ihnen 7#,, und x,, so erhalt man fiir gz) als Funktion von & mit 3) = A x 
ae 

‘1 x3 — 9 x? — ——=- * S> at WN 20 

Al x 9x Al x oe gets (20) 
Diese Gleichung hat fiir positive & drei reelle Wurzeln, eine negative und zwei positive. Fir negative 
& wandern die beiden positiven Wurzeln wegen des fiir Amplituden der GréSenordnung Eins sehr 
groBen Zahlenfaktors von & in dieser Gleichung rasch ins Komplexe ab, wahrend fiir groBe Betrage 
von & die negative Wurzel in erster Naherung gegeben ist durch 


325s I 
x= 2.58 |g] +a]: 
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po SPE, gy ee 


Andererseits liefert unsere Grundlésung nach (13), die ja fiir groBe Werte von Ig| noch brauchbar ist, 
wenn man in sie 


34 AB 1 
He gilt also hier rg a [a (21) 


16 4 25 
=e, a=Z(lt+ggbet--) 


eintragt, also setzt | 
226 3 4. 34 | 
Ama Ci as ieee bey oh —— ig = een ene | 


2 


fiir ®, nach (13) 
3 
Ge zal 


64 G52) 


Nery poe 
26 [ee eee 


Somit gilt fiir é << 0 und grofe Betrage von & | 
4 3 
are ee 

Beide Entwicklungen stimmen also im Leitterm iiberein. Da die Grundlésung bis zu lE,| | 
= 0,13224 A? konvergiert [vgl. w. u. Gl. (27)], miissen fiir groBe |E| beide Entwicklungen tiberhaupt 
iibereinstimmen, also nur Umordnungen voneinander sein. Das ist fiir das Glied ~ 1/? auch noch 
leicht zu zeigen, wenn man den expliziten Ausdruck fiir 3. nach Potenzen von 1 /é entwickelt. 
Pzq enthalt ja einen Term der Form f?/(4 — 5 «)®, der: offensichtlich ~ 1/& ausfallt. Aber dieser 
explizite Nachweis ist iiberfliissig. da ja nach (27) die Konvergenzgrenze der Grundlésung bekannt 
ist. Entsprechendes gilt fiir den Bereich rechts von dem kritischen Wert von a. 

Fiir groBe positive € hat unsere Gleichung (20) eine kleine positive Wurzel, die in erster Naherung 
genau wie oben gegeben ist durch 

se Zk 


Ps0 = 9.55 pce fe (23) 


Diese Lésung schlieBt nun wiederum glatt an die Grundlésung in dem Bereich 4/5 < a < 12/13 
an, wie aus (22) unmittelbar folgt. Denn die Grundlésung ist ja in diesem Bereich fiir hinreichend 
kleine 8 wiederum konvergent, wenn nur gilt (als Ungleichung hinsichtlich der GréBenordnung in f) 


Bee (24) 


Die beiden anderen Wurzeln werden im Falle groBer positiver & ebenfalls betragsmaBig sehr groB; 
man erhilt fiir sie in erster Naherung 


oe ee VAS 1 


2 1 il 
(P0)3 = = i aE = A a . 
Die negative Wurzel (79); entspricht der Fortsetzung der Grundlisung aus dem Bereich a < 4/5 
in dem Bereich « ~ 4/5 und zeigt, daB dort beliebig groBe Amplituden des Gliedes sin 37 auf- 
treten, dieser Zweig der Grundlésung also gegen die Oberschwingung konvergiert, die dem Grenzfall 
A —> 0 entspricht. 
Um das Verhalten der Grundschwingung an dieser Stelle voll iberblicken zu kénnen, bestimmen 
wir noch den kritischen Wert von &, fiir den die beiden positiven Wurzeln in eine Doppelwurzel zu- 
sammenfallen. Mit den iiblichen Methoden der Algebra findet man 


E) = — 0,13224 42, X19 = 0,37288 , X%, = — 0,52630 . (26) 


Wenn man |pz9| itiber a auftragt, erhalt man Abb. 1, wobei der MaBstab auf der a-Achse so gewahlt 
ist, das £8, das ja von der Ordnung f ist, als endliche Strecke erscheint. Aus Abb. 1 entnimmt man 
zunachst auf Grund bekannter Satze der Reihenlehre, da® die Konvergenzgrenze unserer Grund- 
lésung genau an der Stelle 


(25) 


4 5 4 
a= = — a 0,13224 f A? = — —e (27) 
liegt, sie also fiir x = 4/5 + ¢ wiederum konvergent wird, falls man der zweiten Singularitat bei 
= 12/13 nicht zu nahe kommt. 
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Weiter entnehmen wir aus Abb. 1, daB keine Verzweigung der Lésung vorliegt, denn auf dem 
von links herkommenden Ast ist P30 < 0, wahrend auf dem anderen, bei a = 4/5 — ¢ einsetzenden 
Ast 39 > 0 ist. Der Schnittpunkt kommt nur dadurch zustande, daB wir den Betrag von Qs) auf- 
getragen haben. Vielmehr gibt es in diesem Bereich von a zwei periodische Lésungen mit der Fre- 
quenz der Grundschwingung, die einen so groBen Anteil der Oberschwingung enthalten, daB sie sehr 
rasch nicht mehr von dieser unterschieden werden kénnen; nur bei der dritten Liésung klingt dieser 
Anteil der Oberschwingung sehr rasch ab. Aber diese ,,echten‘“‘ Grundschwingungen kénnen beim 
-Durchgang durch die kritische «-Stelle nicht stérungsfrei ineinander iibergehen, da das Vorzeichen 
VON 9 bei diesem Durchgang wechseln muB. In 
dem Bereich hart links von a = 4/5 enthalten 
alle drei Lésungen ungefahr gleich starke Anteile 
der Oberschwingung, und es ist a priori nicht klar, 
ob und welche dieser drei Lésungen instabil sind. 
Genauen Aufschlu8 dariiber kénnte natiirlich nur 
eine Stabilitatsrechnung bringen. Vermutlich wird 
wohl die Grundschwingung, wenn man von links 
her mit « gegen die kritische Stelle geht, in die 
ja sicher stabile Oberschwingung iibergehen, also 
kein Sprung auf den unteren Ast der Kurve ¢) > 0 
stattfinden. Kommt man jedoch von rechts her 
an die kritische Stelle heran, so mu die Lésung 
auf dem unteren Ast irgendwann instabil werden, ead 
spatestens natiirlich fiir « = 4/5 —e und in die pate ee ee Beane Preeue 
Lésung des Astes ~) < 0 umschlagen. 

Den Verlauf der Schwingung iiber die erste Viertelperiode, wie er sich aus der Uberlagerung der 
beiden endlichen Bestandteile der Lésungen ergibt, zeigen fiir « = 4/5 —«, also fiir den kriti- 
schen Wert von « Abb. 2 und 3, und zwar Abb. 2 fiir die positive Doppelwurzel, Abb. 3 fiir die zu- 


gehérige negative Wurzel gemaB (27). Beide Bilder zeigen deutlich, daB sowohl 8B wie ( zweimal ver- 
schwinden, die Abb. 3 deutlicher, weil hier der Anteil der Oberschwingung gréfer ist. 


Ap 


Abb, 2. Verlauf von @ und @ fiir die Doppelwurzel bei Abb. 3. Verlauf von # und 9 fiir die dritte, negative Wurzel 
a = 4/5 —« wabrend der ersten Viertelperiode. bei « = 4/5 —¢ wahrend der ersten Viertelperiode. 


Es ist hier der Ort, kurz auf die Verhaltnisse beim Doppelpendel gleicher Massen und gleicher 
Fadenlangen einzugehen. Wenn auch unsere Gleichungen (1) nicht die Bewegungsgleichungen des 
Doppelpendels bei Beriicksichtigung aller Glieder dritter Ordnung darstellen, sondern nur fiir 


B =—1/6 den Anteil der Schwere in dieser Naherung richtig wiedergeben, so kann man doch 
erwarten, daB der Vorgang wenigstens qualitativ richtig erfa8t wird. Und da fiir gleiche Massen gilt 
a = = w 0,107 = + — 0,093, (28) 


V2 
so sieht man sofort, daB fiir Amplituden A nach (27) mit 


ZAPZ0T+16, AZY672 x 26 (28b) 


bereits der kritische Bereich erreicht wird. Man wird also vermuten diirfen, daB bei diesem 
Doppelpendel etwa fiir Ausschlage der GréBenordnung /2 in der Schwingungsform der Grund- 
schwingung sich bereits ein merkbarer Anteil der Oberschwingung zeigen wird. 
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Bevor wir uns zur Auflésung der Divergenz unserer Grundlésung fiir sehr kleine « wenden, 
wollen wir uns noch kurz mit dem Falle « = 3/5 beschaftigen, bei dem das Verhaltnis der F're- 
quenzen (49:9 gleich 1:2 ist. Hier lage also, wenn die Grundschwingung ein Glied mit sin 2 t 
enthalten kann, wiederum Resonanz mit dem Grundterm der Oberschwingung vor. Wir wollen 
zeigen, daB in der Umgebung dieses «-Wertes in der Tat neben unserer Grundlésung auch Lésungen 
mit der Frequenz der Grundschwingung existieren, die auch die geraden Sinus-Glieder enthalten; 
aber im Gegensatz zu dem eben untersuchten Falle werden wir zeigen, daB hier eine echte Ver- 
zweigung der Liésung vorliegt. 


4. Der Fall « = 3/5. Da die Rechnung hier prinzipiell genau so verlauft, wie im zweiten Ab- 
schnitt, wollen wir uns méglichst kurz fassen. Wir setzen also # und jetzt als Fourier-Reihen mit 
allen Sinusgliedern an, wobei neben ©, = A und @, jetzt auch ®, und @, jeweils einen endlichen 
Term enthalten sollen, waihrend alle iibrigen ®, und O, mindestens von der Ordnung f sind. Es gilt 
also 


at = +E, 1+x=2 +64 ; 
ned @,=2ALPOut..., (29) 
©, =Yo tPPat--- >» 0, = Boy +B Oy +.--> > Oo) eee 

Damit erhalten wir fiir die sechs Koeffizienten %,, 941, Yoo « - » 95, nach der gleichen SchluBweise 


wie Ziff. 3 die folgenden Gleichungen, in der die Glieder der Ordnung f voll beriicksichtigt sind: — 


} 


A+B On +4) +354] +54 +B Ont 4 Bg) +3 4 Ob] =O, 


—(¢+%8)(F4 +8 On) +4 + GBA +2 49%) = 05 


4{3 9 9 ibe 3 A\2 GO 
ss (7 +8) i + 35 Fao + B(Ega0+ Pas + 35n)| + Pn + BOs) 4 4 BM + 2/ 5 )n| =0. 


5 ; 
—4(5 +B) [Pn + Poo + B (Oar + Pe1)] + Poo + B Gar 4 +B (P% +24? p20) =0. 


Die Glieder nullter Ordnung liefern 
5 


24. 
7 ote a =O, —A+A=0, 
5 9 5 (31) 
5) (20 x50) + Bon = 0, ras) (P29 + Poo) 4 Po = 9, 


wobei die ersten beiden Gleichungen wieder identisch erfiillt sind, wahrend die beiden anderen iiber- 


einstimmen, und liefern 


Oog + is 0), eae — = an : (31a) 
Die Glieder erster Ordnung liefern, wenn wir Igy = —> gap gleich einsetzen 
ee ice 2(5) (4 +2 Agi) =0, 
By + oem A 2(5) (4? + 2 Ag}) = 90; 
vat =n a 3é P20 7574 P20 + ; (5 y [P20 + 2 A? Poo] = 0, _ 
Pat son + Py a = [P20 + 2 A? Goo] = 0. 


J 
wobei die beiden letzten Gleichungen nur wieder vertraglich sind, wenn 


dy aN 3:[, ah Al eeanean | 
95 “1 P20 =z § Pao +i5| = (5 | (P20 + 2 A? Peo) « (33) | 


| 


{ 
\ 
i 
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Eliminieren wir aus dieser Gleichung und den beiden ersten Gleichungen (32) x, und 3, und setzen 
Po = A x, so erhalten wir fiir x 


25 6 34 : 
pee ee pet (ia) 0, (34) 
also entweder 
x= 0 
oder 
155 é 
Pc ee 
Ce esos ak G2) 
Somit werden die beiden anderen Wurzeln nur reell und dann entgegengesetzt gleich, wenn 
RPA 
Fiir € = £, werden auch die beiden nicht-trivialen Wurzeln Null; also tritt fiir 
9 9 36-17 
Met eon ey (37) 


eine echte Verzweigung der Liésung auf; die dem Verzweigungsast (Abb. 4) zugeordnete periodische 
Lésung enthalt einen fiir positive € sehr rasch wachsenden Anteil der Oberschwingung, ist also sehr . 
bald nicht mehr von der Oberschwingung zu unterscheiden. Schon fiir = 0 haben die beiden end- 


¥20 
a 
J 
eo) 
Abb. 4. Die Amitude ¢., iiber « in der Abb. 5. Verlauf von # und g wahrend der ersten Halbperiode 
Umgebung der Verzweigung bei fiir den dem Werte « = 3/5 zugeordneten Wert von 2) auf dem 
a = 3/5. Verzweigungsast. 


lichen Bestandteile der Fourierentwicklung fiir ? und @ gleiche Amplituden. Aus Abb. 5, in der fiir 
diesen Fall # und ¢ fiir eine Halbperiode aufgetragen sind, ersieht man, daB wahrend einer Halb- 


periode sowohl & wie @ zweimal verschwinden. Nach Verlauf einer Viertelperiode liegt jetzt kein 
-Umkehrpunkt der Bewegung mehr vor, sondern beide Massen haben gleichen Ausschlag, aber ent- 
gegengesetzte Bewegungsrichtung. Diese Verzweigung diirfte beim Doppelpendel gleicher Massen 
nicht auftreten, da eine ahnliche Rechnung wie oben fiir die dazu erforderliche Amplitude liefert 
A> A.2 also A > 1. : 

Genau wie eben kann man in allen anderen Fallen vorgehen, wo die Frequenzen der Grund- und 
Oberschwingung der kleinen Schwingungen in ganzzahligem Verhaltnis stehen. Wir wollen dabei drei 
Falle unterscheiden. 

a) Das Verhiiltnis Wo: @49 ist eine ungerade Zah]. Hier wird qualitativ immer dasselbe auftreten, 
wie im Falle n = 3, den wir oben ausfihrlich besprochen haben: es gibt keine Verzweigung der 
Lésung, sondern es treten jeweils drei periodische Lésungen auf, von denen zwei einen so starken 
Anteil der Oberschwingung enthalten, daB sie, von einem schmalen Streifen der Ordnung f um den 
kritischen Wert von « herum abgesehen, nicht von der Oberschwingung zu unterscheiden sind. 


300 Schmieden: Nichtlineare Schwingungen bei zwei Freiheitsgraden Ingenieur-Archiv 


b) Das Verhaltnis gp :@,p ist eine gerade ganze Zahl. Hier gibt es immer eine echte Verzweigung der 
Lésung von der Art, wie wir es eben fiir n = 2 besprochen haben. Eine Lésung enthalt keinen Anteil 
der Oberschwingung, die beiden anderen enthalten die Oberschwingung, und zwar auBerhalb des 
B-Streifens so stark, daB sie wie im ersten Falle nicht von der Oberschwingung zu unter- 
scheiden sind. 


c) Das Verhiltnis 59:49 ist rational, aber nicht ganz. Wiederum ergibt sich eine echte Ver- 
zweigung, mit den gleichen Higenschaften wie im zweiten Fall. Wir haben den einfachsten Fall 
(9? @19 = 3:2 durchgerechnet, wollen aber auf die Wiedergabe dieser Rechnung verzichten, da sie 
prinzipiell nichts Neues bringt, Nur soviel sei erwahnt, daB man, um die Verzweigung der Lésung zu 
erhalten, fiir ® undgy Fourierreihen ansetzen muB, die nach Vielfachen von 1/2 fortschreiten, also 


2 : 5.8 
& = Oj)/2 sin > + Q, sint + O3/2 sin aces 
(38) 


as : eet: 
CD = Pi )2 sin > + A sint + Os)2 S185, es a 


Dabei werden auf dem Verzweigungsast, also fiir die Lésungen, die einen Anteil der Oberschwin- 
gung enthalten, die Glieder O,, 3/2, D3/2 endlich, wahrend alle anderen @, und ®, mindestens von der 
Ordnung f ausfallen. 


5. Die periodischen Lésungen im Grenzfalle sehr kleiner a. Es bleibt noch die Singularitat un- 
serer Grundlésung fiir sehr kleine « zu klaren. Wir fragen dazu zunachst, welche Aussagen die Theo- 
rie der kleinen Schwingungen in diesem Grenzfall zu machen gestattet. Lat man in (3) «> 0 
gehen, so fallen Grund- und Oberschwingung zusammen, also 1) = Wo9 = @) undy = A sin Wot 5 
= 0. Es gibt also nur noch eine Normalschwingung, da ja die Bewegungsgleichungen jetzt lauten 


i +020=0, 


eee (39) 

O+o+earp=0. 
Damit ist ? entkoppelt und kann sich unabhangig von gy bewegen. Gibt man aber ? einen auch noch 
so kleinen Ausschlag, so schwingt es mit der Frequenz w,, steht also mit dem unteren Pendel, das 
ja ebenfalls wy als Kigenfrequenz besitzt, in Resonanz und liefert demnach fiir y die Lésung 


: Boot : ; 
yg =A sin Wo t— —J*- COS Wo t mit 7 = Bicin ot, (40) 


deren Amplitude im Laufe der Zeit beliebig gro®B wird. Also mu8 in der Tat B = 0 sein, also & 
identisch verschwinden, wenn man in diesem Fall eine Normalschwingung erhalten will. Wir wollen 
untersuchen, wie sich die Verhaltnisse im nicht-linearen Fall gestalten, vor allem feststellen, ob 
sich beim Vollzug des Grenziiberganges « —> 0 in der Lésung selbst auBer der trivialen ? = 0, die 
natiirlich auch im nichtlinearen Fall auftritt, noch eine zweite periodische Lésung ergibt, mit einer 
von W, verschiedenen Frequenz. 


Dazu entwickeln wir die Loésungen der Gleichungen (11) nach Potenzen von « und behalten je- 
weils nur die Glieder héchster Ordnung in a bei. Um das Bildungsgesetz dieser Glieder zu erkennen, 
muff man in der Berechnung bis zu den Gliedern in f* gehen. Wir wollen hier nur die Ergebnisse 
dieser an sich einfachen, aber recht zeitraubenden Rechnung angeben; man erhalt 


3 3 


Mi eae 8) ae Oy =ZAR(L+a)+--, 
P32 = 2 (8 Oa) eos Bye = gag BAB, 
“= — ig tet te aed (41) | 
ea ee 


XXYV, Band 1957 Schmieden: Nichtlineare Schwingungen bei zwei Freiheitseraden 301 


Es gilt also fiir den Term jeweils héchster Ordnung in x"! die Relation 
Ait Ua Oi m % —Ax,, fir i= 23 AN (42) 


Beschranken wir uns also auf diese Terme héchster Ordnung in a4, so erhalten wir dafiir aus den 
Gleichungen (11) die folgenden Gleichungen: 


a Di, m— % A Hy, — (4 Cee + Oss rhea ote gay yg) OS 
th = A um = (4 aa = Xo mee i 7 os Hm—1 P41) x 0 : 


Denn wenn auch in den von f 3° bzw. By herrithrenden Gliedern allmablich Terme mit Potenzen 
von a + auftreten, so beeinflussen diese Terme doch nicht die Glieder der héchsten Ordnung in a 1 
fiir 3, ,, und x,,. Nehmen wir die Relation (42) bis m — 1 als bewiesen an, so zerstéren sich in der 
Klammer das erste und letzte Glied, und der Vergleich beider Gleichungen zeigt, daB Y, ,, % —A Xins 
also (42) auch fiir m gilt. Setzt man also O1,m © —Ax,, in die erste Gleichung ein, so sicht man, daB 
fiir m = 4 fiir die x,, die Funktionalgleichung — als GréBenordnungsrelation — besteht 


(43) 


20 Mm We Hine + Hg %mg °° * t+ Mm 9%, (m4). (44) 
Setzen wir nun 
Y = He B? -+ Hy BF +--+, (45) 
so wird offenbar 
ay? = 15 B* + (Hg %y + Hg %e) B® +> , (46) 
so daB die Funktionalgleichung (44) gelést wird durch 
2 a (y — x2 B® — x5 B®) = y* . (47) 
also 
y=a— ort 2 pat —% ps ae (48) 


(das Plus-Zeichen der Wurzel scheidet fiir « > 0 aus). 
Damit erhalten wir 
eal 
1l+«e 
O,=aA+2pA(L+0)—cA+ Al = pA +o0)+4y/ +[04. 


ee 


5 +B A?(L—a) t}a—f =14+2649(1—a)—/ + [04], 
(49) 


Ist «a <0, so muB bei der Wurzel das andere Vorzeichen gewahlt werden. Diese Teilsummation 
ist zweifellos gestattet, solange die Wurzel nach fallenden Potenzen von « entwickelt werden kann, 
also solange |x| > 3/8 8 A® ist. Da aber (48) auch einen Sinn behilt, wenn diese Grenze fiir  unter- 
schritten wird, so wollen wir versuchsweise den Grenziibergang « —> 0 vollziehen. Dabei diirfen wir 
wegen der iibrigen Vernachlassigungen natiirlich die Losung, wenn iiberhaupt, nur einschlieBlich 
der Terme in f benutzen. Wegen des doppelten Vorzeichens der Wurzel erhalten wir so zwei Lésun- 
gen, eine langsame Schwingung 


ye seed VN C= AGI te, 9 =f APsint 


und eine rasche 
L+e=1fop4, g=Asnr, $=0. 


Die zweite ist die triviale, mit # =0 und der Duffingschen Differentialgleichung fiir p, deren 
Frequenz auch den richtigen B-Term liefert. Die langsame liefert die Frequenz wo, aber fiir 0 eine 
Amplitude der Ordnung 6. Wir wollen jetzt zeigen, daB diese Lésung in der Tat fiir das nicht- 
lineare System im Grenzfall « — 0 als rein periodische Lésung existiert, indem wir ansetzen 


O,=BPI,+ PO. +:-::, 1l+x=1-+ », p+, pt4+---:, 
Os = B25 +--+ D, = 3 Gag ++: 
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Die Gleichungen (11) liefern dann, mit a? = 0, in der tiblichen Weise durch Koeffizientenvergleich 


3 ; 3 AB 03 
Ong Oe ie A ee a lem Dos = — 3p ° 2 
3 A2\3 1/3 A2y4 (60) 


Damit ist gezeigt, daB auch fiir diesen Grenzfall zwei rein-periodische Lésungen auftreten, deren 
Fourier-Entwicklungen fiir hinreichend kleine f konvergieren. 


(Eingegangen am 3. November 1956.) 
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triebwerke. Kolbentriebwerke und, ihre Kennlinien. Aufbau der Kolbenmotoren und ihre Einzelteile. 
Strémungstriebwerke und ihre Kennlinien. Aufbau der Gasturbinentriebwerke und ihre Einzelteile. — 
VI. Flugleistungen. Waagrechtflug. Steigflug. Reichweite. Messung der Flugleistungen. — VII. Steuerung 
und Stabilitit. Flossen, Ruder, Steuerungen. (Wirkung der Ruder. Lings- und Querstabilitit. Auto- 


matische Flugzeugsteuerung. — VIII. Abflug und Landung. Hilfsmittel fiir Abflug und Landung. Der ne 


Abflugvorgang. Der Landevorgang. — IX. Das Drehfliigelflugzeug. Allgemeines iiber Schraubenflugzeuge. 
Der vertikale Flug. Der Vorwartsflug. Blatthewegungen. Hubschraubersteuerung, Stabilitat, Schwin- 
gungen. — X. Kursbestimmung und Flughbeférderung. Einfluf des Windes auf Kurs und Flugdauer. 
Funknavigation. Die direkten Kosten der Flugbeférderung. Einige neuere Flugzeugmuster. — Zusammen- 
stellung der wichtigsten Bezeichnungen. — Namen- und Sachverzeichnis. — Zahlentafeln. 


Die sechste neubearbeitete Auflage der bekannten v. Misesschen Fluglehre nimmt ebenso wie ihre -Vor- 
ginger eine Sonderstellung in der Fachliteratur ein. Das Buch vermeidet einerseits die Oberflachlichkeit 


popularer Darstellungen und weicht den begrifflichen Schwierigkeiten, die in der Natur der Sache liegen, — 

nicht aus, erfordert jedoch andererseits yon dem Leser nur ein Minimum an Vorkenntnissen. ‘Die physi- | 
kalischen Prinzipien des Flugwesens werden von Grund auf entwickelt, ohne dabei die formalen Methoden 

' der Differential- und Integralrechnung zu verwenden. Fiir Tragflachen- und Luftschraubenberechnung — 


werden elementare aerodynamische Theorien zusammen mit verallgemeinerten Versuchsergebnissen heran- 
gezogen. Die Darstellung idealisierter thermodynamischer Prozesse der Kolbentriebwerke, der Gas- 
turbinen und der Diisentriebwerke wird erginzt durch eine, Beschreibung der. wesentlichen Einzelteile 
und Funktionen dieser Triebwerke. Die Flugleistungen werden sowohl vom Standpunkt ihrer Voraus- 
berechnung wie auch ihrer Ermittlung im Fluge behandelt. Die Ausfithrungen liber Steverung und Stabili- 


tat umfassen die Grundlagen der automatischen Flugzeugsteuerungen. Die aerodynamische Theorie der: 


Drehfliigelflugzeuge schlieBt den vertikalen Flug und den Vorwartsflug ein, An Sonderproblemen werden 


erértert: Abflug und Landung, Kursbestimmung (Funknavigation) und Flugbeférderungskosten. Da trotz 


der Beschrinkung auf die elementarsten mathematischen Hilfsmittel wo immer méglich bis zur quanti- 
tativen Behandlung der Fragen vorgedrungen wird, kann das Buch auch dem angehenden Tngenieur, der 
sich der Flugtechnik erst zuwenden will, als erste Einfiihrung dienen. 
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